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ABSTRACT 


That  a  self— gravitating  equilibrium  perfect  fluid  in  empty 
space  has  a  spherical  configuration  if  it  is  nonrotating  or  has  a 
unique  configuration  characterized  by  some  physical  parameters 
if  it  is  slowly  rotating  is  considered  physically  evident  but  has 
not  been  rigorously  derived  from  Einstein’s  field  equations 
together  with  suitable  asymptotic  conditions.  We  therefore 
consider  the  equations  for  a  general  relativistic  stationary 

spacetime  which  is  asymptotically  flat  (and  dif  f  eomorphic  to  iR  'j 
and  consists  of  an  exterior  vacuum  solution  and  an  interior 
perfect  fluid  in  rigid  motion  with  a  given  equation  of  state  p(p) 
satisfying  O^p^p,  given  total  mass  m  and  given  surface 
temperature  Tt  .  Global  analysis  techniques  are  used  to  show 
that  there  is  no  family  of  such  static  (nonrotating)  spacetimes 
which  depends  dif ferentia bly  on  a  parameter  and  contains  the 
spherically  symmetric  solution  except  for  families  containing 
only  solutions  dif f eomorphic  to  the  spherically  symmetric  one 
In  the  general  stationary  (rotating)  case,  if  we  require  that  the 
solutions  be  close  to  the  static  spherically  symmetric  ones  (in 
the  sense  of  a  suitable  topology  on  the  set  of  stationary 
spacetime  metrics)  we  can  show  that  there  is  at  most  a 
one  — dimensional  differentiable  family  of  global  axially 
symmetric  solutions  parameterized  by  the  value  of  the  angular 
momentum  J  except  for  those  obtained  by  coordinate 


l  v 


transformations.  A  recent  result  of  Lindblom  (1976)  which 
shows  that  all  such  spacetimes  are  axisymmetric  makes  the 
latter  result  quite  general. 
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NOMENCLATURE 


Throughout  this  dissertation  we  will  use  geometrized  units 
in  which  the  speed  of  light  c=l  and  Newton’s  constant  G  =  1  . 
Greek  indices  will  run  from  0  to  3  except  when  they  are  used  as 
n  — tuple  indices.  The  usage  will  make  this  clear.  Small  Latin 
indices  will  run  from  1  to  3  while  capital  Latin  indices  will  run 
from  1  to  2.  Tensors  will  be  written  sometimes  with  and 

sometimes  without  indices  when  no  danger  of  confusion  exists. 
The  set  of  reals  will  be  denoted  by  iR  and  the  integers  and 
complex  numbers  by  K  and  (C,  respectively. 
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CHAPTER  I 


EQUILIBRIUM  RELATIVISTIC  FLUIDS 


1.1  Introduction 

Consider  a  star.  The  word  conjures  up  an  image  of  a 
distant,  hot  radiating  ball  of  matter.  A  closer  look  reveals 
many  different  types  of  stars  with  radically  different 

compositions  ranging  from  ordinary  matter  stars  (which 
include  the  planets)  to  dense  neutron  stars.  Among  these  there 
are  many  different  populations  with  wide  variances  in 
luminosity,  temperature,  mass,  angular  momentum  and  density 
corresponding  to  differences  in  their  internal  structure  such  as 
differences  in  composition,  nuclear  and  chemical  reactions, 
convective  motion  and  differential  rotation.  But  all  these 
quantities  are  determined  by  the  star’s  history  of  evolution, 
starting  perhaps  with  a  gravitationally  collapsing  cloud  of  low 
density  matter  whose  gravitational  energy  is  converted  into 
thermal  energy,  heating  up  regions  in  which  different  reactions 
which  produce  more  heat  can  take  place.  The  higher  pressures 
thereby  produced  can  slow  or  perhaps  stop  the  collapse.  But 
how  can  one  determine  a  star’s  structure  this  way?  Disney  (1976) 
has  suggested  that  the  structure  of  an  evolving  protostar  is 
very  dependent  on  the  initial  conditions  and  numerical  analysis 
of  the  gravitational  collapse  of  clouds  lends  credence  to  this  (c.f. 
Buff  et  al  (1979)).  Yet  even  though  there  is  in  some  sense  an 
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infinite  number  of  such  initial  conditions,  when  we  look  at  the 
stars  they  seem  to  be  very  similar  in  their  configurations, 
namely,  they  are  almost  spherical.  For  example,  the  sun,  with 
an  equatorial  rotation  period  of  about  25  days  and  the  interior 
rotating  slightly  faster,  has  polar  and  equatorial  radii  that 

differ  by  only  about  1  part  in  105  (c.f.  Hill  and  Stebbins  (1975)) 

Physically,  the  reason  is  very  intuitive.  All  nuclear  and 
chemical  reactions  proceed  towards  their  endpoints  as  the  star 
evolves,  with  electromagnetic  radiation  carrying  off  any  excess 
heat.  Viscosity  damps  differential  rotation  Perturbational 
analysis  of  a  spherical  star  indicates  that  gravitational 
radiation  carries  away  the  energy  of  all  but  the  1=0  and  1=1 
spherical  harmonic  components  of  the  perturbation  (c.f.  Thorne 
(1967)).  Thus  a  star  gradually  evolves  towards  some  equilibrium 
non  — radiating  body  and  can  be  described  adequately  by  a 
perturbation  of  this  equilibrium  model  for  some  latter  portion 
of  its  life,  even  though  its  evolution  is  very  complex. 

But  what  constitution  of  materials  and  what  spatial 
configurations  can  such  an  endproduct  of  stellar  evolution  have 
It  is  well  known  that  there  are  three  main  classes,  the  white 
dwarfs  supported  by  electron  degeneracy  pressure,  the  neutron 
stars  supported  by  neutron  degeneracy  pressure  and  black  holes 
the  end  result  of  total  gravitational  collapse. 

The  black  hole  equilibrium  configurations  are  well 
understood.  There  is  no  internal  structure  to  worry  about  and 

the  no  hair  uniqueness  theorems  of  Israel  (1967)  and  Carter 
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(1971)  show  that  their  gravitational  field  is  uniquely  given  by 
their  mass  m  and,  if  rotating,  their  angular  momentum  J. 

(Similar  no  hair  theorems  hold  if  the  black  hole  has  a  charge 
(c.f.  Israel  (1968)).  However,  even  though  most  stars  probably 
have  magnetic  fields  and  are  neutral,  dealing  with  charged  stars 
or  those  with  magnetic  fields  is  not  a  trivial  extension  of  the 
present  work  so  throughout  this  dissertation  we  will  consider 
only  stars  with  no  charge  and  no  magnetic  fields.) 

We  are  here  interested  in  the  non  — singular  equilibrium 
stellar  models.  There  are  no  such  configurations  for  large 
masses  (generally  for  masses  greater  than  about  two  solar 
masses)  but  whether  a  star  collapses  totally  depends  on  its 
equation  of  state  and  its  history  as  well  as  its  mass.  Some  large 
mass  stars  eject  some  mass  before  getting  to  the  endpoint  of 
evolution,  thereby  avoiding  total  collapse.  We  will  therefore 

assume  that  the  mass  of  a  configuration  is  always  in  [0,mcrit) 
where  mcrit  corresponds  to  the  critical  mass  where  total 
gravitational  collapse  occurs.  We  can  again  ask  the  question, 
what  are  the  equilibrium  configurations  of  such  stellar  models? 
To  determine  all  equilibrium  solutions  is  quite  hopeless.  The 
purpose  of  this  dissertation  is  to  determine  the  general 
properties  which  any  solution  corresponding  to  a  nonrotating  or 
slowly  rotating  isolated  equilibrium  stellar  model  must  possess. 

Fortunately,  as  long  as  the  density  of  matter  is  much  less 

than  10  gm/cm  ,  which  occurs  for  all  the  proposed  realistic 
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equations  of  state,  the  physical  interactions  in  a  star  occur  on 
two  widely  differing  distance  scales:  the  short  scale  nuclear  and 
chemical  interactions  and  the  long  scale  gravitational 
interaction  (c.f.  Thorne  (1967)).  (This  is  not  true  for  magnetic 
fields  but  as  we  noted  before,  we  will  not  be  dealing  with  this 
case.) 

This  enables  the  study  of  equilibrium  stellar  models  to  be 
divided  into  two  parts:  one  studying  the  inter-particle 
interactions  on  a  small  scale  to  determine  the  constitution  of 

cold  (non  — radiating)  catalyzed  matter,  neglecting  gravitation, 
and  the  second  studying  the  large  scale  structure  of  the  matter 
and  gravitational  field,  treating  the  matter  as  a  smoothed  out 
fluid  distribution.  A  review  of  small  scale  properties  of  matter 
for  the  first  study  is  given  in  Zeldovitch  and  Novikov  (1971).  In 
the  following  we  take  up  the  second  study  for  non  — rotating  or 
slowly  rotating  isolated  stars.  More  specifically,  we  will 
investigate  the  global  symmetries  and  the  uniqueness  of  such 
configurations.  Unfortunately  the  problem  of  existence  of  such 
configurations  seems  to  be  even  more  difficult.  This  will  be 
mentioned  occasionally  with  a  fuller  discussion  of  the 
difficulties  being  given  at  the  end  of  the  last  chapter. 

We  start  in  the  next  section,  then,  with  an  investigation  of 
thermodynamic  equilibrium  stellar  models.  It  will  be  seen  that 
global  stationarity  will  be  a  reasonable  assumption  to  make  for 
such  stellar  models  so  we  will  then  present  some  formalisms  for 
stationary  spacetimes  which  will  be  very  useful  for  our  later 
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work.  (Stationarity  rules  out  any  radiative  processes.)  The 
appropriate  boundary  conditions  resulting  from  the  assumption 
that  we  are  dealing  with  an  isolated  system  will  also  be  given. 

Although  both  non  —  rotating  Newtonian  stellar  models  (c.f. 
Lindblom  (1978))  and  non-rotating  black  holes  (Israel  (1967))  are 
spherically  symmetric,  i.e.  they  admit  the  group  S0(3)  as  a 
group  of  isometries  with  the  group  orbits  being  spacelike 
2  — surfaces  (which  are  then  necessarily  of  constant  positive 
curvature),  this  has  not  been  proved  for  non-rotating 
equilibrium  relativistic  stellar  models.  It  is  intuitively  so 
physically  obvious  though,  especially  because  of  the  Newtonian 
and  black  hole  results,  that  it  has  been  assumed  and  used  in 

many  studies  of  stellar  structure  (c.f.  Thorne  and  Zytkow  (1977)) 
as  well  as  in  studies  of  gravitational  collapse  (for  example,  the 
Zeldovitch  sequence,  c.f.  Harrison  et  al.  (1965)).  In  the  slowly 
rotating  case  the  corresponding  uniqueness  theorems  must  be 
weaker  in  the  sense  that  one  expects  uniqueness  (with  some 

OC  0  o 

physical  parameters  fixed)  onl}^  for  configurations  close  to 
spherical  symmetry  with  small  angular  momentum.  Thus  they 

00  90  00  00 

will  be  local  not  global  results. 

In  chapter  2,  then,  we  take  a  look  at  these  historical 

uniqueness  theorems  and  see  what  sort  of  directions  one  might 

try  to  take  for  the  relativistic  stellar  model  case.  For  reasons 

00  0  0 

which  will  appear  clearer  later  on,  our  results  are  local  for 
both  the  rotating  and  non— rotating  cases.  They  are  thus  weaker 
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than  the  historical  uniqueness  theorems  for  the  nonrotating 
case  but  are  more  or  less  equivalent  to  those  for  the  slowly 
rotating  case.  In  fact  our  results  are  valid  for  very  general 
equations  of  state  but  cannot  be  used  to  increase  the  generality 
of  the  Newtonian  slowly  rotating  case  because  the  Newtonian 

case  is  more  degenerate  than  the  relativistic  case.  This  will 
become  clearer  later. 

We  will  then  outline  the  approach  which  we  will  take 
towards  proving  that  for  a  fixed  equation  of  state,  a  fixed 
surface  temperature  and  a  fixed  mass  the  corresponding  isolated 
equilibrium  solutions  of  the  gravitational  field  equations  are 
unique  (up  to  coordinate  transformations)  provided  the  solution 

9  9  9  9 

is  close  in  a  differentiable  sense  to  a  spherically  symmetric 
solution.  The  reasons  for  the  latter  condition  will  become  more 

apparent  later.  In  particular,  the  solution  will  be  spherically 
symmetric  if  it  is  not  rotating  and,  if  slowly  rotating,  will  be 
uniquely  determined  by  the  small  angular  momentum.  Since 

99  99 

we  must  be  close  to  spherical  symmetry  we  then  take  a  look 
at  spherical  spacetimes.  This  is  followed  by  an  investigation  of 
the  conditions  on  the  equations  of  state  wThich  we  must  impose. 

It  is  seen  that  they  are  not  very  restrictive. 

Chapter  3  consists  of  an  introduction  to  the  mathematical 
techniques  wThich  we  use  to  obtain  our  results.  Chapters  4  and  5 
present  the  statements  and  proofs  of  the  uniqueness  theorems 
for  the  non  — rotating  and  slowly  rotating  isolated  equilibrium 
stellar  models,  respectively. 


a  •  i  »  t-  t-  **  •••••. b 
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1.2  Thermodynamic  equilibrium 

In  a  fluid  in  which  the  reaction  rates  are  too  slow  to  be 
important  on  the  time  scales  of  interest  the  chemical 
composition  of  the  fluid  can  be  assumed  to  be  fixed  uniquely  by 
two  thermodynamic  variables  such  as  the  number  density  of 
baryons  n  and  the  entropy  per  baryon  s.  This  is  obviously  the 
case  for  a  fluid  in  equilibrium  so  let  us  assume  that  we  have 
such  a  simple  relativistic  fluid  with  viscosity  and  heat 
conduction.  The  fluid  mechanics  which  we  present  here  is 
standard  and  can  be  found,  for  example,  in  Misner,  Thorne  and 


Wheeler  (1972)  §22  and  Weinberg  (1972)  p.53ff. 


The  stress  energy  tensor  for  such  a  fluid  with  unit 


four— velocity  ua  (  u°ua=-l  )  is 


(1.1) 


wThere  p(n,s)  is  the  total  mass  energy  density  in  the  rest  frame 
of  the  fluid,  p(n,s)  is  the  isotropic  pressure  in  the  rest  frame, 
£(n,s)>0,  77(n,s)>0  are  the  bulk  and  shear  viscosities 

respectively, 


O 


(1-2) 


(13) 


are  the  expansion  and  shear  deformation  rates  respectively 
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l  u  v „  )>  8  g  +u  u  is  the  projection  operator  orthogonal 

to  u  and  q  is  the  heat  flux.  Round  brackets  around  tensor 
indices  will  refer  to  symmetrization  while  square  brackets  will 

refer  to  antisymmetrization.  The  angular  velocity  tensor  cj  is 

given  by  — P  ^P^  Vj-^u^j  We  assume  that  qa  obeys  the  general 
relativistic  version  of  the  Fourier  law  of  heat  conduction,  first 
proposed  by  Eckart, 


q“=-Kp“'J0(sT+Tu/()  (1.4) 

where  /c(n,s)>0  is  the  coefficient  of  heat  conduction 

Furthermore,  we  take  over  the  standard  thermodynamical 
and  particle  conservation  laws.  Bar  yon  conservation  implies 

Va(nua)=0-  (1-5) 

while  the  first  law  of  thermodynamics  is  given  by 

dp=n  1(yo+p)dn+nTds  (1.6) 

where  the  temperature  T(n,s):=n  (dp/d s)  .  As  well  as  the  third 

law,  namely  if  T=0  then  s=0,  we  require  that  all 
thermodynamic  quantities  are  positive.  Then  the  Bianchi 
identity  together  with  the  Einstein  equations 


(1.7) 
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yield  the  conservation  laws  VaTw^=0  which,  when  projected 
parallel  and  orthogonal  to  ua,  yield 

Va(pua+qa)=-p@+2  7]  vap<Jap+(d2-uaq° 

(p  +  p)ua=-Pa/?a/Sp  +  P“v^(.faP^+  Tj  (T^-2q^u^). 

(The  convention  we  use  for  defining  the  Riemann  curvature 
tensor  is  that  used  in  Hawking  and  Ellis  (1973)  and  is  such  that 

vaV'M-VaH=Rh(3x''.) 

These  equations  correspond  to  the  Newtonian  energy 

conservation  law  and  Euler’s  equation  respecti vely.  The  first 
can  be  rewritten  as 

TVa(nsu<x+T  \a)  =  ^62+7.T1G<xpGa^qaqa/KrY  (1.8) 

Defining  the  entropy  four-vector  to  be  sa.  =  nsua-fT  1qa  the 
second  law  of  thermodynamics  takes  the  form  VasQ^0  with 
equality  iff  the  fluid  is  in  thermal  equilibrium.  (In  the  fluid’s 

rest  frame  Vas  is  the  rate  at  which  entropy  is  being  generated 
per  unit  volume.) 

As  Lindblom  (1976)  has  noted,  thermal  equilibrium  then 
implies  the  following  properties.  Since  the  fluid  is  shear  free 
and  expansion  free,  by  (1.8),  it  is  rigid  in  the  sense  that 
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P^V^-gV^g^+u^O  (19) 

This  notion  of  rigidity  is  equivalent  to  the  definition  given  by 

$  0 

Born,  Herglotz  and  Noether:  A  body  is  called  rigid  if  the 

distance  between  every  neighboring  pair  of  particles,  measured 

with  respect  to  the  world  line  of  either  of  them,  remains 

0  0 

constant  along  the  world  line.  (c.f.  Trautman  (1965))  If  in 
addition  there  is  no  rotation,  i.e.  cJa;S=  0,  then  Vau^— uaii£=0  so 
that  ua  is  hypersurface  orthogonal, 

u[«V„r° 

Furthermore,  qa  vanishes  (also  by  (1.8))  so  the  fluid  is 
perfect,  ua=— 9alog  T  and  all  thermodynamic  variables  are 
constant  along  ua, 


T=n=s=0.  (1.10) 

Euler’s  equation  then  becomes 

dlog  T=(p+p)  !dp,  (1.11 ) 

the  equation  of  hydrodynamic  equilibrium.  Integrating  this 
shows  that  there  is  an  arbitrary  integration  parameter  Tb,  the 
value  of  the  temperature  at  the  surface  of  the  star  given  by 
p=0,  which  does  nothing  but  change  the  temperature  of  the  star 
by  a  constant  amount;  so,  when  discussing  the  uniqueness  of 


i 


' 
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configurations  we  will  always  take  Tb  to  be  fixed.  Using  (l .11) 
in  the  exterior  derivative  of  (1.6)  it  is  easily  seen  that  dpvdp^U 

so  the  fluid  is  barotropic,  p=p( p).  Defining  ©a=T-1ua  it  is 

easily  verified  that  the  rigidity  condition  and  (1.10)  imply  that 

is  an  infinitesimal  symmetry  of  the  spacetime  structure  and 
all  thermodynamic  variables  in  the  region  occupied  by  matter 
(but  is  not  defined  in  vacuum). 

An  analogous  result  has  been  obtained  by  Stewart  (1971) 
from  a  statistical  mechanics  approach.  Using  the  relativistic 
Boltzman  equation  with  the  underlying  physical  assumption  that 
the  fluid  is  not  too  cold  or  dense,  so  that  particles  which  are 
about  to  collide  have  uncorrelated  momenta,  he  shows  that  in  a 
collision  dominated  equilibrium,  if  one  component  of  the  fluid  is 
not  massless,  the  fluid  must  be  rigid  and  the  spacetime  is 

locally  stationary  under  8°. 

It  should  be  noted  that  if  T  =  0  so  that  the  fluid  is 

trivially  in  thermodynamic  equilibrium,  results  analogous  to 
those  of  Lindblom’s  above  still  hold.  Namely,  the  fluid  is  rigid, 
barotropic  with  an  equation  of  hydrodynamical  equilibrium 

dp=(  p  +p)d(log  n)  (1.12) 

and  (  p  +p)  nu^  is  an  infinitesimal  symmetry  of  spacetime  and 

all  thermodynamic  variables.  Defining  Te=n~  (p+p)  we  then  see 

that  dp  =  (  p  +  p)d(log  Te)  and  that  ©“=Te  u°  is  an  infinitesimal 
s y m metry  of  spacetime  so  all  our  results  could  be  extended  to 


1 


■ 


12 


cold  matter  at  T=0  by  substituting  Te  and  6e  for  T  and  © 
respectively.  Although  cold  matter  is  often  assumed  in 
astrophysical  investigations,  particularly  of  equations  of  state 
(c.f.  Zeldovitch  and  Novikov  (1971)),  this  is  not  an  expected 
physical  situation,  especially  in  light  of  some  formulations  of 
the  third  law  of  thermodynamics  (c.f.  Pippard  (1957)  p.5l).  In 
the  following  we  will  assume  T>0. 

A  stationary  spacetime  is  defined  as  one  which  admits  a 
globally  timelike  Killing  vector  field,  along  which  all  physical 
fields  are  Lie  transported.  Such  a  non— singular  spacetime  which 
admits  a  Cauchy  surface  and  is  asymptotically  Euclidean  (which 
we  will  define  precisely  later)  can  be  shown  to  be  in  a  state  of 
thermodynamic  equilibrium  (Lindblom  (1976)).  Note  however 
that  this  is  not  quite  the  converse  of  the  preceding  results  since 
there  we  have  only  local  stationarity  in  the  region  occupied  by 
matter. 

It  would  seem  reasonable  to  conjecture  that  for  an 
isolated  equilibrium  fluid  in  empty  space  the  spacetime  region 
outside  the  matter  should  be  stationary  as  well  as  the  interior 
region.  This  will  clearly  depend  on  the  asymptotic  conditions 
imposed,  but  is  not  easy  to  prove  rigorously  even  if  one  demands 
pseudo— stationarity,  i.e.  existence  of  a  global  Killing  vector  field 
which  is  asymptotically  timelike.  However,  that  the  conjecture 
is  physically  reasonable  is  enforced  by  a  result  of  Friedman  and 
Schutz  (1975)  which  shows  that  stars  with  an  ergosphere  (a 
region  in  which  the  asymptotically  timelike  Killing  vector  field 
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becomes  spacelike)  are  unstable  to  a  non— axisy mmet ric  mode. 

On  the  other  hand  it  is  a  stronger  conjecture  than  one  can 
justify  in  Newtonian  theory  (where  stationarity  implies  all 
variables  describing  the  fluid,  including  the  velocity,  are 
independent  of  time)  since  the  non— axisy mmetric  rigid  rotators 
such  as  the  Jacobi  ellipsoids  are  examples  of  non— stationary 
thermal  equilibrium  fluids.  But  in  general  relativity  such 
objects  would  radiate  gravitational  waves  and  hence  would  not 
be  expected  to  be  in  thermal  equilibrium. 

For  solar  sized  black  holes  the  time  scales  for  approaching 
stationarity  are  of  the  order  of  milliseconds  (Carter  (1972))  so 
stationarity  is  a  good  assumption.  Also,  Thorne  (1969)  has 
shown  that  quadrupole  vibrations  in  neutron  stars  are  damped 
by  gravitational  radiation  in  the  order  of  one  second  while 
Langer  and  Cameron  (1969)  have  argued  that  many  other 
vibrations  are  damped  via  nuclear  reactions,  which  create 

thermal  heating,  on  the  same  sort  of  time  scale.  However  even 
after  a  star  has  reached  its  final  state  as  a  white  dwarf  or 
neutron  star,  the  rate  at  which  chemical  and  thermal 
equilibrium  is  reached  and  at  which  differential  rotation  is 
damped  out  is  very  slow.  For  example,  in  rapidly  rotating  white 
dwrarf  stars  Kippenhahn  and  Mollenhoff  (1974)  have  shown  that 

5 

differential  rotation  is  damped  out  in  a  time  scale  of  10  years. 

This  is  slightly  faster  than  the  time  scale  of  10  years  in  which 
thermal  equilibrium  is  approached.  Having  made  the 
assumption  of  thermodynamic  equilibrium  however,  we  see  that 


■ 
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little  is  lost  in  assuming  stationarity  with  £a ,  say,  being  the 
global  timelike  Killing  vector  field. 

The  vector  field  ©”  must  either  be  proportional  to  or 
©°  and  (Q  must  be  non  — parallel  commuting  Killing  vector 
fields.  In  the  former  case  we  have  Uj-a^]=0  which  in  the  case  of 
a  spacetime  containing  a  perfect  fluid  (and  obeying  Einstein’s 
equations)  is  equivalent  to  the  material  staticity  condition 


(1.13) 


Lichnerowicz  (1955)  has  shown  that  if  the  spacetime  is 
asymptotically  flat.,  asymptotically  source  free  and  topologically 
Minkowskian  (which  our  stellar  models  are  as  we  will  see  when 
we  discuss  boundary  conditions  and  spherical  solutions)  then 
(1.13)  is  equivalent  to  the  metric  staticity  condition 

f[aVvr°  (114> 

which  means  and  hence  ©°  are  everywhere  orthogonal  to  a 
family  of  spacelike  hypersurfaces  by  Frobenius’  theorem 

In  the  latter  case  the  interior  of  the  star  is  axisymmetric. 
Lindblom  (1976)  has  shown  that  by  using  the  Cauchy — Kowaleski 

theorem  and  assuming  that  |a  is  C4,  the  spacetime  is  C&  and 

the  boundary  of  the  star  is  smooth  enough  ©”  can  be  extended 
analytically  into  the  exterior  region.  This  extension  is  linearly 

independent  from  and  commutes  with  and  an  asymptotic 


I,  I 
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argument  indicates  that  it  must  be  a  combination  of  a  time 
translation  and  a  rotation,  so  the  spacetime  is  axisymmetric. 
Unfortunately  it  is  difficult  to  rigorously  match  the  above 

procedure  with  our  boundary  conditions  so  we  will  take  the 
above  result  as  an  indication  that  there  is  little  loss  in 

generality  when  we  assume  a  priori  that  our  stationary, 
non  — static,  stellar  models  are  axisymmetric,  i.e.  there  exists  a 
global  one  —  pa rameter  isometry  group  whose  orbits  are  closed 
spacelike  curves  and  which  commutes  with  the  timelike 
isometry  group.  A  result  of  Carter  (1970)  shows  that  there  is  no 
loss  in  generality  in  assuming  this  commutativity.  We  will  let 

77°  be  the  Killing  vector  field  generator  of  the  closed  spacelike 
orbits. 

From  now  on  we  will  consider  our  stellar  models  to  be 
isolated  stationary  thermal  equilibrium  solutions  of  Einstein’s 
equations  wThich  are  either  static  or  axisymmetric.  We  have 
seen  that  they  must  then  be  perfect,  rigid,  barotropic  fluids.  In 
the  next  section  we  will  present  a  convenient  formalism  for 
■writing  the  equations  governing  such  fluids  and  present  the 
boundary  conditions  one  imposes  to  ensure  that  the  fluid  is 
isolated. 

1.3  Stationary  spacetimes 

Let  (M,g)  be  a  stationary  spacetime.  For  simplicity  let  us 

OO 

assume  that  M  is  C  .  The  differentiability  of  g  across  the 
star  boundary  will  depend  on  the  form  of  the  equation  of  state 
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P(p)  as  p  >0  and  its  differentiability  in  the  interior  will 
depend  of  course  on  the  differentiability  of  p(p)  with  respect  to 
p  for  p>0  .  We  will  see  that  we  can  assume  the  Lichnerowicz 

(1955)  junction  conditions,  namely  that  g  is  at  least  C1, 

3 

piecewise  C  (the  second  and  third  derivatives  are  continuous 
except  at  a  finite  number  of  hypersurfaces  where  they  have 

finite  limits  on  both  sides).  The  metric  will  actually  be  C3  in 
the  exterior  region  so  by  a  result  of  Muller  zurri  Hagen  (1970)  it 
can  be  taken  to  be  analytic  in  the  exterior.  Lindblom  (1978)  has 
shown  that  if  p{ p)  is  analytic  so  is  g  in  the  interior  (for  rigid 

rotation).  We  also  assume  that  (  is  C  piecewise  C  .  With  the 

1  3 

Killing  equations  this  implies  that  (  is  C,  piecewise  C  . 

The  form  for  the  metric  g  which  we  will  derive  below  was 
originally  given  by  Ehlers  (1958)  but  the  geometric  derivation 
presented  here  first  appeared  in  Klinzle  and  Savage  (1980a). 

Since  M  is  stationary  the  flow  of  the  vector  field  £  corresponds 
to  the  orbits  of  a  one-dimensional  isometry  group  and  these 

2 

orbits  are  timelike  C  submanifolds  diffeomorphic  to  &  since 
there  are  no  closed  timelike  lines  (causality  condition).  Let 

e2U:=-g(£>£)  (115) 

1  3 

so  that  £(U)=0  and  U  is  C  ,  piecewise  C  . 

We  assume  that  the  strong  causality  condition  (Hawking 

and  Ellis  (1973),  p.192),  which  has  been  shown  to  be  valid  for 
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physically  realistic  solutions,  holds  for  our  stellar  models. 
Namely,  for  all  x0GM  there  exists  a  neighborhood  V  of  x0  that 
is  intersected  by  each  timelike  orbit  in  at  most  one  connected 
segment.  This  is  true  for  the  space  of  orbits  £  of  £  iff  £  has 
a  manifold  structure  such  that  the  canonical  projection  tt.  M 

■ — >£  is  C2  (Palais  (1957)  p.20).  Thus  M  is  a  principal  fibre 
bundle  over  £  with  structure  group  (&,+).  (See  for  example 
Choquet— Bruhat  et  al.  (1977b)  p.128  for  this  and  p.288  for  the 
theory  of  connections  used  below.) 

We  orthogonally  decompose  TXM  into  a  1-dimensional 
vertical  subspace  Vx  parallel  to  £  and  a  complementary 

horizontal  subspace  Hx  by  writing  ^£TXM  as  =  +  with  <f°= 

_2U 

e  g(<^,£)  •  This  decomposition  defines  a  connection  since  it  is 

invariant  under  the  group  action  of  (&,+)  .  It  can  be 

equivalently  characterized  by  a  connection  form,  i.e.,  in  this  case 
a  real  valued  1-form  co  that  satisfies 

(i)  1.^=0  (ii)  gj(Xx)=0<==>Xx2Hx  (ill)  u(£)=l  (1.16) 

These  equations  define  co  uniquely  if  the  horizontal  subspace  is 
given.  The  curvature  form  H  associated  to  the  connection  form 
co  is  then 


H=dco+coAco  =  dco 


and  H=7i  H  for  a  unique  2-form  H  on  £,  since  £JH-0  and 
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so  that  H  is  baselike.  Similarly  we  have  U  =  7r*fj. 
We  define  a  Riemannian  metric  7  on  E  by 


7(x>t  )|x-e2Ug(x',^)]x  (1-17) 

for  any  xCtt  \x),  where  x>t^TxM  are  the  unique  horizontal 
lifts  of  respectively.  Then  we  find 


g--7T  (e2^ )co®co+u  (e  ^7).  (1.18) 

If  (M,g)  is  also  static,  then  £Ad£=0  where  ?=g(^)  by  (1.14). 

^  ^ 

Using  (I.16iii)  it  is  seen  that  ^=-rr(e  )co  so  cjaH  =  0  which 

rsj 

implies  H=0.  The  converse  is  also  true. 

We  now  drop  the  tilde  and  consider  only  the 

three-geometries  (E,7,cj,U)  assumed  globally  defined,  C  in  the 

2 

exterior  region,  C  except  at  the  star  boundary  with  7  a 
positive  definite  Riemannian  metric,  co  a  connection  form  on  an 

tR-principal  bundle  over  E  and  UGC^E.tR)  .  Note  that  there  is 
no  reason  for  the  closed  2  — form  H  to  be  exact 

On  a  local  neighborhood  of  x£M  there  exist  coordinates 
(t,x1)=xa  such  that  $=dt,  i.e.  £^g=0<=>dtg/Jlf-0  .  Letting 


2U  2U  -2U  2U 

goo=~e  gok=-e  ak  gki=e  7ki“e  aka! 


(1.19) 


we  have  that 
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«=dt+akdxk  7  =7kldxk®dx]  Hi.:=29^a.j  (1.20) 

We  now  use  such  local  adapted  charts  to  derive  the 
three-geometry  formulas.  First  we  observe  that  for  any  vector 

field  x  such  that  =  X°=g(£,x)  and  X*^  are  a  scalar  and 

vector  field,  respectively,  on  E.  By  writing  all  the  tensorial 
equations  with  the  0  index  lowered  and  the  Latin  indices  raised, 
the  ak  do  not  appear  except  through  We  use  to  lower 

and  raise  all  Latin  indices  (which  will  always  be  considered  to 
run  from  1  to  3). 

Einstein’s  equations  for  a  perfect  fluid, 

R^-TRg^n^P+prXVPSW  (121) 

then  become  in  the  3-dimensional  formulation 

8ij=27irarU7]SasU+^e4Uh1hJ+(/)+p)e_4UT2010J 

-(2pe_2U+(p+p)e~4UT2a2)7ij,  (1.22) 

AU:=7ijViajU=M-^e4Uh2+(p+p)e"4UT202,  (1.23) 

£ijkVj(e4Uhk)=2(p+p)e"2UuT0i,  (1.24) 

where  h1:=£1^kHjk,  h  i^y^h  hJ,  9  :  =  ©  ,  6  \=y ^6  9\  M:  =  ^(p+3p)e 

2  2U  22  6 

and  u:=-uo>0  so  that  v  =e  +  T  6  .  On  the  other  hand,  ©  is  a 
symmetry  of  the  spacetime  metric  and  all  thermodynamic 
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quantities  iff 

^eU=0  ^e7  =  0  £ep  =  lep=I6T=0  (1.25) 

The  rigidity  conditions  together  with  Euler’s  equation  take  the 
form 


0JH+d(uT  'e  sv)=0 


(1.26) 


which  implies  that 


c£6H  =  d(0JH)  +  @JdH  =  O  whence  £6h  =  0. 


(1.27) 


Since  H  is  closed  we  also  have 


Vrh  =0. 


(1.28) 


If  M  is  static  H  =  0  and  equation  (1.24)  implies  8  =0  . 
Einstein’s  equations  then  reduce  to 


Rii=-2aiuaiu-2p7i, 

1J  1  J  1  '  Ij 


(1.29) 


AU—  M 


(1.30) 


__  2U 

where  p:  =  pe  while  (1.26)  becomes 

T=TceUc'U  (1.31) 

in  the  interior,  where  the  subscript  c  will  always  refer  to  the 
value  at  the  center.  We  will  see  later  that  for  all  the  cases  we 
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are  interested  in  there  will  be  a  uniquely  determined  center. 
The  equation  of  hydrostatic  equilibrium  (l.ll)  then  becomes 

dp+(p+p)dU=0.  (1.32) 

We  can  now  state  the  boundary  conditions  which  are 
usually  imposed  to  restrict  to  solutions  which  are  isolated, 
namely,  we  demand  that  spacetime  is  asymptotically  flat  at 
spacelike  infinity  in  the  following  sense  (c.f.  Lichnerowic-z  1955). 
In  terms  of  the  three  geometry, 

(i)  there  exists  a  compact  KCE  and  a  dif feomorphism  ip:E\K — : 
3 

iR  \B  where  B  is  a  closed  ball  centered  at  the  origin; 

3 

(ii)  with  respect  to  the  standard  coordinate  system  on  iR 


T'^ij+Oflxr1)  dk7ij=0(lxr2) 

u=o(|xi_1)  aku=o(lx|"2) 

hi=0(|xr2)  (1.33) 


where  lx|2=  E  (x1)2. 

i=l 

These  conditions  are  implied  by  more  sophisticated 

definitions  (c.f.  Geroch  (1972)  for  the  special  case  of  E^tR3  which 
we  will  see  later  will  hold  for  the  spacetimes  we  are  interested 
in.  (We  do  not  require  any  conditions  at  null  infinity  since 
stationarity  has  already  ruled  out  any  dynamics.)  With  such 
asymptotic  conditions  Lichnerowicz  ((1955)  p.126  et  seq.)  has 
shown  that  in  the  static  case  (130)  together  with  p^ 0,  p>0, 


a 
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implies  that  U^O  on  £  and  if  there  is  a  compact  region  DC£ 
such  that  p=p  =  0  in  V=£\D  then  U  has  no  maximum  or 
minimum  in  V  and  no  critical  point  in  some  neighborhood  of 
00  unless  the  space  is  flat.  Furthermore  Dc:=5x€!£|U(x)^c|  is 

compact  for  all  c<0  and  Sc=dDc=U  *(c)  is  dif feomorphic  to  the 
Euclidean  2  sphere  for  c  sufficiently  close  to  0,  i.e.  £  is 
asymptotically  spherical. 

If  £  is  homeomorphic  to  tR3  (we  will  see  later  that  £  is 

3 

actually  dif f eomorphic  to  tR  for  all  the  situations  we  are 

interested  in)  the  Poincare  lemma  (c.f.  Choquet-Br uhat  et  al. 

(1977b)  p.213)  implies  that  the  closed  form  H  is  exact.  Thus, 

even  though  the  ai  are  not  necessarily  the  components  of  a 
1-form  on  £  we  can  find  a  1-form  a  on  £  such  that  H  =  da  . 

This  a  is  determined  uniquely  up  to  a  closed  1  — form  which  by 
the  Poincare  lemma  must,  be  exact,  i.e.  a  is  determined  up  to 
the  differential  of  a  function  f  on  £  .  Imposing  a  gauge 

condition  would  give  us  a  unique  a  since  then  Af=0 

which  implies  f=const.  However,  we  will  not  impose  this  as  an 
extra  condition  (the  reason  for  this  will  be  apparent  later). 
Rather,  the  gauge  freedom  will  be  treated  in  a  manner 
analogous  to  the  coordinate  freedom. 

We  now  want  to  define  the  mass  and  angular  momentum 
of  our  isolated  stellar  model.  There  is  a  difficulty  in  general 
relativity  in  introducing  such  physical  quantities  which  describe 
the  structure  of  a  gravitating  system  as  a  whole.  This  results 
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from  the  fact  that  in  flat  space  such  quantities  are  associated 
with  the  action  of  the  Poincare  group  but  in  a  curved  space  the 

corresponding  translation  subgroup  is  much  larger  than  the 
constant  vector  fields  one  has  in  flat  space.  This  is  due  to  the 
fact  that  one  must  admit  as  asymptotic  Killing  vector  fields  all 
those  vector  fields  which  are  linear  in  position  with  coefficients 

that  form  a  skew  tensor  (  rotation  )  so  the  translation 
subgroup  includes  all  vector  fields  which  added  to  this  do  not 
disturb  its  asymptotic  behaviour.  The  faster  the  spacetime 
metric  approaches  the  Minkowski  metric  the  closer  the 

asymptotic  symmetry  group  will  be  to  the  Poincare  group,  but 
the  smaller  the  amount  of  information  available  from  the 
asymptotic  behaviour  of  the  gravitational  field  will  be.  The 
general  situation  is  therefore  very  difficult  but  in  our 
particular  case  (stationary,  axisymmet ric,  asymptotically  flat 

and  E  homeomorphic  to  iR  )  we  can  define  quantities  which  in 
an  asymptotic  expansion  of  the  metric  are  seen  to  correspond  to 
the  physical  mass  and  angular  momentum  (c.f.  Misner,  Thorne 
and  Wheeler  (1972)  §19). 

Ehlers  (1965)  has  shown  that  the  total  gravitational  mass 
of  an  asymptotically  flat  stationary  system  in  which  S  is 

3 

homeomorphic  to  iR  can  be  defined  by 

m  =  (l/47T)/aD(V1U  +  Jre4',Hlkak)dS,  (134) 

where  dS,  is  the  normal  surface  volume  element  of  dD  and 


■ 
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where  D  is  any  domain  of  E  containing  all  the  matter.  It  is 
easily  verified  from  Einstein’s  equations  (i. 22— 1.24),  the  rigidity 
conditions  (1.26)  and  the  asymptotic  conditions  that  the  value  of 
the  integral  is  independent  of  the  choice  of  D  as  well  as  a.  In 
the  static  case  this  reduces  to  the  usual  integral 

m  =  (l/ 4:Tt) j  •  pMdV  obtained  from  integrating  the  relativistic 
Poisson  equation  (1.30),  where  dV  is  the  volume  element  of  E 
(c.f.  Ehlers  et  al  (1962)  p.68).  A  calculation  shows  that  thi  s  mass 
agrees  with  that  found  from  the  ADM  four  — momentum  for  an 
asymptotically  flat  stationary  spacetime  (c.f.  Jang  (1979))  and 
also  agrees  with  the  coefficient  of  — l/ix|  in  an  asymptotic 
expansion  of  U,  i.e. 


U=— m/lx|+0(|xl  2)  (1.35) 

at  infinity. 

A  general  definition  for  the  angular  momentum  in  the 
general  case  is  somewhat  trickier.  (See  for  example  Geroch  (1972) 
or  Ashtekar  and  Streubel  (1979).)  However,  in  the  axisymmetric 
case  a  formula  corresponding  to  (1.34)  can  be  found  Recall  that 
r)  is  the  axisymmetric  vector  field  on  E,  so 

7  =0,  £v  U  =0,  £va=0  whence  t£7?H  =  0  (1.36) 

and  in  the  asymptotic  region  we  can  take 


77=x1d2-x2d1+0(lxl  !) 


(1.37) 


. 
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in  terms  of  the  cartesian  coordinates.  By  fixing  77  in  this 
manner  we  are  eliminating  some  of  the  coordinate  freedom  so 
we  will  expect  the  slowly  rotating  isolated  equilibrium 
stationary  stellar  models  to  form  a  1  —  parameter  family  rather 

than  a  3— parameter  family.  Now,  T)a—{ 0,77*)  is  a  4— Killing  vector 
field  that  commutes  with  £=dt  and  is  everywhere  spacelike. 
Defining  the  amgular  momentum  J  by 

167TJ:=JSDV-|g|V°;p1(511J2k3dxJAdxk  =  JDV -|glR°^dx1Adx2Adx3(1.38) 

wThere  |gi=det(g^y),  it  is  clear  that  J  is  independent  of  D  as 
long  as  D  contains  all  the  matter.  This  can  also  be  written  as 

J=(l/16i7  )J 3DS-^Hlk(?7k+J-e4Jaka1?)I)  +  ak(2?7k7,e?J'7kl)aIujdS1 

(1.39) 

which  can  be  verified  to  be  independent  of  the  choice  of  a. 

1.4  (2+1)— dimensional  formalism 

The  formalism  we  present  here,  valid  in  the  case  where  U 
has  only  one  (nondegenerate)  critical  point,  was  used  by  Kilnzle 

(1971)  in  an  analysis  of  the  linearized  static  Einstein  equations 
because  it  allowed  a  result  for  2-dimensional  Riemannian 

manifolds  of  positive  curvature  to  be  used.  (This  result,  with  a 
corrected  proof,  is  given  in  appendix  1  )  We  need  this  formalism 
so  that  Klinzle’s  analysis  can  be  used  as  well  as  to  aid  in  the 
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analysis  of  the  linearized  stationary  Einstein  equations.  When 

we  look  at  the  spherically  symmetric  solutions  we  will  see  that 

U  then  has  only  one  critical  point,  which  is  nondegenerate  and 

defines  a  unique  center  xc  such  that  U(xc)=Uc,  diU(xc)=0  and 

<9i9jU(xc)^0.  For  technical  reasons  we  will  be  forced  to  consider 

$  $  ©  © 

only  spacetimes  -which  are  dif ferentiably  close  to  the  spherical 

solutions  so  that  this  will  hold  true  for  the  situations  we  are 

interested  in.  We  are  then  able  to  re-write  our  equations  in 

terms  of  the  two-dimensional  geometry  of  the  equipotential 

surfaces  which  will  all  be  Riemannian  2-spheres. 

2  2 
Let  S  be  an  abstract  2-sphere  and  ic:  S  — >E  be  the 

imbedding  map  of  S2  into  E  such  that  ic(S2)=U  1(c)=:Sc  for  any 


c£(Ue,0)  .  Then  any  £  1-form  is  characterized  by  cj:  =  i  o  and 


u°:  =  i*(VUJcj). 


(1  40) 


This  decompositon  extends  to  any  E— tensor  field  since  E  is 
Riemannian.  An  intrinsically  defined  normal  derivative  to  the 
Sc  surfaces  is  given  by 


D«=i*(|VU!  2Vvuco) 


(1.41) 


while  the  second  fundamental  form  of  Sc  is 


n=i*(v(|vul  Vu)) 


(1.42) 


(c.f.  Kobayashi  and  Nomizu  (1963)) 
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To  find  the  local  coordinate  expression  let  (x  ,A=1,2)  be  a 

chart  of  S2  and  ic:  x  A^(U=c,xA=  x  A(c))  (using  (U.xA)  as  a  chart 
of  E  ).  Then  for  a  1-form 

6JA=a)1<9U/ dx  A+cjbc)xB/3x  A-coA, 
w°=i"(VUJW)=i"(coi71JajU)=co1.  (1.43) 

In  particular 


7ab=~ab  700=7ii31U3JU=:W-2  7°  =  rlrSrU?iA=aAU  =  0  (1.44) 

and  is  used  to  raise  and  lower  all  indices  of  S  — tensors.  If 

A  _ AB 

we  also  -write  y1A=PA  and  P  :=y  PB  we  obtain 


/ 

/  w 


(7ij)= 


"  +  P  h 

f  A 


(1.45) 


Although  P^  does  not  transform  as  a  tensor  on  S  (if  xA — »  x  A= 


x  (U,x)),  by  treating  PAaA  as  a  vector  field  on  S  we  find  that 
for  any  S  -tensor  K  (with  possible  U  -dependence), 


DK=auK-ypK  (1.46) 

is  invariantly  defined  since  DK— DK  can  be  expressed  in  terms 
of  S  —tensors.  For  example,  if  U  is  a  1-form,  DwA-DuA= 

W“3cj°aAW-W  1fiABaJB  .  In  particular, 
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D7AB=2W‘nAB.  (1.47) 

Note  though  that  when  U  has  only  one  critical  point,  VU 
is  a  vector  field  that  vanishes  nowhere  except  at  the  center  and 
is  orthogonal  to  the  surfaces  Sc.  It  can  thus  be  used  to 

construct  a  global  polar  type  coordinate  system  (U,xA)  on  £ 
such  that  PA=0  so  D=du  .  To  simplify  matters  we  will  always 
consider  this  to  be  done  when  dealing  with  this 
(2  +  1)  — dimensional  formalism,  but  it  is  not  necessary.  By 
writing  the  1— indices  up  and  the  A  — indices  down  one  gets  the 
same  equations  which  we  derive  below  and  Kiinzle's  analysis, 
done  with  PA+0,  is  merely  simplified  a  bit. 

After  some  work  it  is  found  that  the  Einstein  equations 
take  the  form 


2 

DW=-0+MlC-±e4u(W-3h°  +  W~‘h  ) 

+  (p+ pJWV+V  (1  48) 

8[AhB]=0  (149) 

DhA=aA(h°Wa)-4hA+2(p  +  p)W,e_6+T-aBeB  (1.50) 

DQAB=2w“1oACQg-W1QnAB-2W3aAvaBw 

+  w  2Va3bW  +  4w  R7ab-Tw  e  hAhB 

-(p+p)w+'4V(eAeB-e27AB)+2pw“1e":2U7AB  d-5i) 


I  |  |  I  II  I  III 
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1e4Uh°hA 


(1.52) 


R-n2+nAB0AB+2W2-^e4U(h2-W  2h°  )=— 2p  (1.53) 


T-S  _ABy-  y-  ytS  — AByy  -tt 


where  h“:  =  r“hAhB,  6  ‘  =7',ue  AS  B,  n:  =  y4BnAB,  TAB:=W£1AB 
=  v  y  5  ab>  Rab="vR7ab  anc^  V  is  the  covariant  derivative  with 


respect  to  the  connection  defined  by  ~y  AB.  Using  a°=a1diU  =  0, 
equation  (1.25)  becomes 


£^¥=0  (1.54) 

D0A=O,  V(a0e)=O  (1.55) 

■while  equations  (1.28)  and  (1.26)  yield,  respectively, 

Dh°-¥“1hODW-W“1hA0A¥  +  VAhA+¥_1nh°=O  (1.56) 

and 


1-^,  /  rp  1  2  U  \  T ,  t  1 _  A  rr—  B  ,  1  2U\  1,  0 _  —pr  3  t 

D(vT  e  )=W  eAB0  h  ,  dA(vT  e  )=^V  h  £AB9  (1.57) 
and  equation  (1.27)  translates  into 

l-gh°=0,  di-gh  A=0  (1.58) 

_ A 

Note  that,  in  fact,  the  hie  derivative  with  respect  to  9  of  any 
S  —tensor  vanishes. 

In  the  static  case  Einstein’s  equations  are  just 


' 


DW=-Q+MW 


(1.59) 


-l 


2w  ,oACn'B-w  'nn^-zw  3aAwaBw 
+  W_2V AdB’W+-jW_1R  7 AB  (1.60) 

vBn®-aAn=o  (i.ei) 

K-o3+nABo4B+2ivz=-2p. 


(1.62) 


CHAPTER  II 


UNIQUENESS  THEOREMS 

2.1  Historical  uniqueness  theorems 

Now  that  we  have  developed  some  formalisms  to  deal  with 
our  isolated  equilibrium  (stationary)  stellar  models  we  would 
like  to  get  an  idea  of  the  types  of  approaches  that  may  lead 
towards  the  uniqueness  results  which  we  mentioned  in  the  first 
chapter,  namely  that  in  the  static  case  the  solutions  must  be 
spherically  symmetric  and  in  the  slowly  rotating  case,  for  a 
given  mass  m,  surface  temperature  Tb,  and  equation  of  state 

p(p)  the  solutions  near  spherical  symmetry  form  a 
3-parameter  family  parameterized  by  the  angular  momentum. 

To  orient  ourselves,  we  will  look  at  the  approaches  used  in  the 
historical  uniqueness  theorems,  then  at  the  general  features 
which  any  uniqueness  theorem  must  have  and  finally,  outline 
the  type  of  approach  which  we  will  take. 

Proving  that  a  static  perfect  fluid  solution  is  spherically 
symmetric  in  nonrelati vistic  theory  is  considerably  simpler  than 
in  the  relativistic  theory.  (Static  here  means  that  the  velocity 

3 

of  the  fluid  vanishes.)  The  metric  of  tR  is  assumed  to  be 
Euclidean  so  that  the  physics  is  described  by  one  function,  the 
gravitational  potential  UN,  which  satisfies  a  Poisson  equation 
and  is  connected  to  the  pressure  via  the  equation  of  hydrostatic 
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equilibrium  dp+  p  dUN=0  .  The  proof  that  the  solution  must  be 


3 

spherical  if  p>0  in  some  compact  domain  of  IB  and  if  UN — >0 
as  x — >oo  was  first  given  by  Carleman  (1919)  for  the  case  of 
uniform  density  models  and  by  Lichtenstein  (1919)  for  arbitrary 
density  models.  Such  a  result  can  be  shown  to  follow  from  the 
theorem  that  a  rotating  Newtonian  stellar  model  must  have  a 
plane  of  mirror  symmetry  which  is  orthogonal  to  the  rotation 
axis  of  the  star.  This  latter  theorem  was  first  proven  by 
Lichtenstein  (1918,1933)  for  the  case  of  uniform  density  stellar 
models  in  rigid  rotation.  It  was  generalized  by  Wavre  (1932)  to 
the  case  of  stationary  axisymmetric  barotropic  ideal  Newtonian 
fluids  in  differential  rotation  and  further  generalized  by 
Lindblom  (1977)  to  the  case  of  barotropic  ideal  fluids  which  have 

stratified  flows  (  v  =  (vx,vy,0),  say).  Wavre’s  proof  involved  the 
Green’s  function  for  the  Laplace  operator  in  the  gravitational 
field  equations  so  it  cannot  be  generalized  to  relativistic  stellar 
models  where  the  field  equations  are  nonlinear.  Lindblom  (1978) 
had  hoped  that  his  proof  could  lead  to  a  generalization. 

However,  his  proof  relies  on  constructing  chords  parallel  to  the 
axis  of  rotation  (z-axis,  say)  with  end  points  on  the  same  level 
surface  of  the  potential  UN  .  Then  a  plane  is  taken  through  zm, 
the  midpoint  of  those  chords  with  the  maximum  z  component 
and  all  functions  are  decomposed  into  even  and  odd  parts  with 
respect  to  reflections  through  this  plane.  The  maximum 
principle  is  then  used  on  the  odd  part  of  the  gravitational 
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potential  UN  which  is  shown  to  satisfy  AUj^O  for  all  z^zm . 

Such  a  coordinate  dependent  proof  is  clearly  difficult  to 
generalize  and  no  one  has  yet  succeeded  in  acheiving  that  end. 

In  the  slowly  rotating  Newtonian  case  the  question  of 
uniqueness  does  not  seem  to  be  as  settled.  A  stationary  viscous 
heat  conducting  Newtonian  fluid  stellar  model  must  have  a 
plane  mirror  symmetry  and  be  axisymmetric.  However  most 
discussions  of  equilibrium  configurations  of  slowly  rotating  stars 
in  the  Newtonian  theory  start  with  an  assumption  about  the 
particular  form  of  the  star  boundary,  such  as  being  ellipsoidal 
or  perhaps  more  complicated,  as  wrell  as  with  assumptions  about 
the  equation  of  state  such  as  the  matter  being  homogeneous  (c.f. 
Chandrasekhar  (1969)).  Neither  ellipsoidal  configurations  nor 
constant  density  make  much  sense  in  relativity.  Lichtenstein 
(1933)  has  shown  that  the  Maclaurin  ellipsoids  are  the  only 
constant  density  configurations  in  a  neighborhood  of  the 
spherically  symmetric  configuration 

There  is  considerable  hope  for  obtaining  uniqueness 

results  in  general  relativity  as  the  general  relativistic  case  is 

less  degenerate  than  the  Newtonian  one.  In  fact  it  is  the 

0  0  0  0 

linearization  of  the  magnetic  part  of  the  gravitational  field  h 
which  does  not  vanish  on  the  spherical  background,  as  we  will 
see,  so  that  there  is  more  hope  of  proving  that  there  exists  a 
finite  dimensional  family  of  slowly  rotating  configurations  with 
a  given  equation  of  state,  given  mass  and  temperature.  We  will 
run  across  difficulties  in  the  existence  part  but  w7e  do  obtain  a 
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uniqueness  theorem  for  small  J.  (The  vector  field  h  is  called 

magnetic  because  it  is  divergenceless  and  its  curl  is  related  to 
a  mass  current  as  is  seen  from  equations  (1.28)  and  (1.24) 
respectively.) 

It  is  considered  physically  evident  that  a  static  general 
relativistic  stellar  model  must  be  spherically  symmetric. 

Indeed,  that  it  is  true  in  the  Newtonian  case  suggests  that  a 
contradiction  must  involve  strong  fields  accompanied  b}7  their 
significant  nonlinearities.  But  Israel’s  (1967)  proof  that  static 
black  holes  must  be  spherical  shows  that  the  non  linear  ities  of 
strong  fields  are  unlikely  to  prevent  the  spherical  symmetry  of 
nonsingular  stellar  models.  More  precisely,  Israel’s  theorem 
states  that  a  static,  asymptotically  flat,  regular  spacetime 
corresponding  to  a  vacuum  solution  of  Einstein’s  equations  must 
be  a  positive  mass  Schwarzsc-hild  solution  if  (i)  past  and  future 
event  horizons  (c.f.  Hawking  and  Ellis  (1973),  p.129)  exist  and 
intersect  in  a  connected  compact  spacelike  topological  2— sphere 
and  (ii)  U  has  no  critical  points  exterior  to  the  event  horizons 
Hawking  (1972)  showed  that  condition  (i)  was  necessarily 
satisfied  even  in  the  stationary  case.  Muller  zum  Hagen  (1973) 
and  later  Robinson  (1977),  by  a  much  simpler  and  more  elegant 
proof,  showed  that  condition  (ii)  could  be  removed. 

The  methodology  of  these  proofs  is  to  derive  expressions 
from  the  field  equations  with  the  form  of  a  divergence  equaling 
a  definite  signed  quantity.  By  integrating  the  divergences  are 
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converted  into  surface  integrals,  the  one  at  infinity  vanishing 
because  of  the  asymptotic  boundary  conditions  and  the  one  at 
the  central  singularity  thus  having  a  definite  sign.  The 
inequalities  obtained  then  yield  the  desired  uniqueness  result. 

In  particular,  Israel  constructs  some  expressions  from  a 
(2  +  l)-dimensional  formulation  of  the  vacuum  Einstein  equations 
with  the  divergence  part  being  the  derivative  of  some  scalars 
with  respect  to  U  .  The  inequalities  obtained  then  imply  that 

d Ap  and  t}AB-^nyAB  vanish.  Robinson  used  a  3-dimensional 
formulation  of  the  vacuum  Einstein  equations  to  obtain  an 

expression  containing  R^R1^  where 

Rijk=2V[kRj]i+i(7i[kaJ]R) 

is  the  conformal  tensor  which  vanishes  iff  7^  is  conformally 
flat  (c.f.  Eisenhart  (1926)  p.89).  Since  one  can  show  that 


RijkR‘ik=e“4U[8w4(nAB-^n7 


AB 


)(Q 


AB  1  n— AB\  ■  —  AEa  wa  wl 

)+y  9AV'aB^J 


(c.f.  Lindblom  (1978)  p  107), a  conformally  flat  solution  of 
Einstein’s  equations  must  have  W,  the  magnitude  of  the 
gravitational  field  strength,  being  a  function  of  U  only.  (The 

proof  actually  broke  up  into  two  cases,  one  yielding  that  R^ 
vanished  and  the  other  yielding  directly  that  W  was  a  function 
of  U  only.) 

Avez  (1964)  had  shown  using  Morse  t  h  e  o  r  3^  that  when  W 
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was  a  function  only  of  U  the  spacetime  must  be  spherically 
symmetric  (even  for  the  nonvacuum  case).  (He  had  thought  that 
he  need  not  assume  that  U  had  only  one  (nondegenerate) 
critical  point  but  there  is  an  algebraic  error  in  the  proof  of  his 
lemma  2,  p.297.  Kunzle  (1971)  generalized  his  result  by  showing 
this  assumption  could  be  eliminated.)  In  fact  when  W  is  a 
function  only  of  U  he  showed  that  U  can  have  only  one 
(nondegenerate)  critical  point  so  that  E  is  then  dif f eomorphic 

3 

to  tR\  Both  of  these  statements  rely  on  a  result  from  Morse 

theory  that  two  surfaces  S  and  SC2  are  diffeomorphic  if  there 

is  no  critical  value  of  U  in  [c^Cg]  (Milnor  (1963)  p.12).  Thus 
Robinson  was  done. 

The  stationary  black  hole  case  was  attacked  by  Carter 
(1971)  in  a  similar  manner  using  an  expression  containing  a 
divergence  and  signed  quantities  but  with  an  important 
difference.  This  expression  was  obtained  from  the  linearized 
equations,  not  Einstein’s  equations,  so  the  proof  is  not  global  as 
the  above  black  hole  theorems  were.  In  fact  Carter’s  result 
states  that,  the  possible  families  of  solutions  of  Einstein’s 
vacuum  equations  which  are  stationary  and  axisymmetric  and 

depend  dif ferentiably  (at  least  C  )  on  some  parameters  form 
discrete  families  wThich  can  depend  on  at  most  2  parameters. 
Hawking  (1972)  showed  that  the  geometry  exterior  to  the  event 
horizon  was  necessarily  axisymmetric  for  stationary  black  holes 
and  Robinson  (1975)  then  showed  that  only  one  such  family 
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exists,  namely  the  Kerr  family  with  J<m2 

As  we  will  see  later  we  also  approach  the  problem  through 
an  analysis  of  the  linearized  equations  and  our  result  in  the 
stationary  case  will  be  of  the  same  rigor  as  the  stationary  black 

hole  theorems.  In  other  words,  both  are  "local”  rather  than 
global  results. 

There  have  been  two  main  attempts  to  prove  that  static 
stellar  models  are  spherically  symmetric.  Lindblom  (1980)  has 
attempted  to  follow  the  approach  of  obtaining  an  appropriate 
expression  from  the  field  equations  which  contained  a 

divergence  term  and  R  ijl:^  as  well  as  other  signed  quantities. 

This  proved  possible  to  do  for  the  case  of  constant  density  but  is 
highly  dependent  on  this  assumption  as  all  the  expressions  must 
contain  the  density  (and  possibly  derivatives  of  it)  in  them.  A 
fairly  thorough  search  for  such  expressions  for  more  general 
equations  of  state  has  been  carried  out  but  with  no  success 
(Lindblom  (1978)).  The  other  approach  is  the  investigation  of 
the  linearized  equations  in  the  (2+l)-dimensional  formalism 
done  by  Kiinzle  (1971).  Here  he  considers  the  problem  of 
proving  that  the  static  spherically  symmetric  stellar  models  are 
isolated  from  other  possible  stellar  models  by  proving  that  there 
are  no  non  — trivial  static  perturbations  of  a  static  spherically 
symmetric  stellar  model  which  leaves  the  central  pressure  pc 
and  central  gravitaional  potential  Uc  unchanged.  These  are 
not  the  appropriate  physical  constraints  to  put  on  the 
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perturbation  but  it  'will  turn  out  in  our  approach  that  the 
above  result  can  be  related  to  a  perturbation  with  the  physical 
constraints  of  constant  mass  ,  constant  surface  temperature  and 
fixed  equation  of  state  p(p) 

Before  passing  to  an  outline  of  our  approach  we  should 
mention  one  other  paper.  Marks  (1977)  claims  to  prove  the 
general  result  that  static  stellar  models  must  be  spherical  but 
his  proof  is  fallacious.  He  assumes  without  justification  that 
the  curvature  of  the  equipotential  surfaces  is  constant.  We  are 
unaware  of  any  uniqueness  or  existence  results  for  the  case  of 
slowly  rotating,  stationary  stellar  models. 

2.2  Our  approach 

Let  us  first  examine  some  of  the  general  features  which  a 

99  99  99 

static  stars  are  spherical  theorem  or  a  stationary  stars  with 

99 

fixed  m,  J,  Tb  and  p( p)  are  unique  theorem  must  possess. 
That  the  matter  must  be  a  fluid  which  cannot  support  stresses 
is  evident  from  considering  a  static  cubic  crystal  or  a  solid 
dumbbell  slowly  rotating  about  the  axis  joining  the  two  bulbs  in 
an  asymptotically  flat  space.  In  fact  the  proof  must  use  the 
fluid  matter  assumption  in  some  global  way  for  some  small 
amount  of  solid  matter  arranged  in  a  lattice  or  the  shell  of  a 
dumbbell  can  destroy  the  uniqueness  because  of  the  long  range 
effects  of  the  gravitational  field.  Thus  the  proof  must  depend 
on  the  fluid  properties  of  the  matter  throughout  the  region 
occupied  by  it 
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The  boundary  condition  at  spacelike  infinity,  namely 
asymptotic  flatness,  is  also  crucial.  If  it  was  replaced  by  a 
non  — spherically  symmetric  boundary  condition  the  static 
theorem  would  clearly  fail  and  since  one  can  only  expect 
uniqueness  in  the  stationary  case  for  slowly  rotating 
configurations  which  are  close  to  spherically  symmetric  this 
type  of  theorem  would  also  fail.  Thus  the  proof  must  also  use 
the  asymptotic  boundary  conditions  in  an  essential  way. 

Furthermore  the  differentiability  across  the  star  boundary 

c  a  c  o 

must  play  a  role.  If  U,  y  and  h  are  too  differentiable  the 
boundary  conditions  will  be  overdetermined  making  the 

existence  of  a  solution  unlikely.  If  they  are  not  differentiable 

enough  one  may  lose  the  uniqueness.  It  is  physically  clear 
that  this  differentiability  must  depend  on  the  equation  of  state 
and  later  in  this  chapter  we  will  analyze  this  in  some  detail. 

What  kind  of  proof  can  fulfil  these  criteria?  As  we  have 
seen,  the  black  hole  uniqueness  theorems  use  the  boundary 
conditions  at  infinity  by  integrating  over  all  of  E  and  using 
divergence  terms  to  convert  some  of  the  integrals  into  surface 
integrals.  However  here  one  does  not  have  to  deal  with  star 
boundary  conditions  or  properties  of  the  fluid.  As  we  noted 
before  and  is  evident  from  the  expressions  of  this  sort  obtained 
by  Lindblom  (1978)  any  proof  along  these  lines  is  going  to  be 
heavily  dependent  on  the  equation  of  state  and  likely  to  be 
possible  for  very  few. 

Despite  the  complications  due  to  the  presence  of  a  source, 
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however,  there  is  a  heuristic  line  of  reasoning  which  gives  some 
hope.  If  the  equation  of  state  is  fixed  and  the  field  equations 
are  supplemented  by  the  conditions  for  rigid  motion  one  obtains 
a  system  of  equations  whose  linearization  may  be  expected  to  be 
elliptic,  once  the  coordinate  freedom  is  factored  out.  The  set  of 
solutions  of  these  linearized  equations  together  with  appropriate 
boundary  conditions  should  then  be  finite  dimensional. 

Our  approach,  then,  is  essentially  a  continuation  of  the 
deformation  method  used  by  Ktinzle  (1971),  namely  an  analysis  of 
the  linearized  equations  of  a  static  (stationary)  perfect  fluid 
with  a  fixed  equation  of  state  p( p),  fixed  surface  temperature 

Tb  and  fixed  total  gravitational  mass  m  (and  fixed  angular 

3 

momentum  J  )  on  an  arbitrary  background  solution  over  tR 
and  reasonably  close  to  the  unique  spherically  symmetric 
solution  with  the  same  p( p),  Tb  and  m.  This  is  done  in  the 
spirit  of  Fischer  and  Marsden’s  (1975a,  1975b)  work  on 
linearization  stability.  Linearization  stability  roughly  means 
that  for  every  solution  of  the  linearized  equations  on  some 
background  solution  there  is  a  corresponding  solution  of  the 
nonlinear  equations.  When  this  is  not  true  the  system  is  said  to 
be  linearization  instable  on  that  background. 

We  give  here  a  general  indication  of  the  argument  in  the 
static  case  since  the  stationary  case  is  very  similar  and  will 
wait  until  we  have  discussed  the  mathematical  machinery  used 
before  we  get  precise.  First  we  want  to  restrict  to  the  study  of 
spacetimes  with  the  appropriate  asymptotic  behaviour  and  the 
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differentiability  conditions  appropriate  to  the  equation  of  state 
p{p)>  which  we  can  consider  fixed.  (The  generality  comes  in 
being  able  to  use  the  same  types  of  spaces  for  different 
equations  of  state  with  a  parameter  determining  the 
differentiability.)  From  the  asymptotic  expansion 

h  m/]xj  +  0(|x|  )  it  is  apparent  that  we  can  also  consider  m  to 

be  an  asymptotic  boundary  condition  so  that  we  can  restrict  to 
spaces  which  will  contain  all  the  solutions  of  mass  m.  We  will 
make  the  appropriate  choice  of  a  Banach  space  X  of  3  metrics 
and  potentials  on  X  writh  the  desired  differentiability  and 
asymptotic  behaviour  from  among  the  weighted  Sobolev  spaces 
(Cantor  (1975a),  (1975b)).  Then,  using  Einstein’s  equations  to 
write  a  nonlinear  differential  operator  mapping  between  Banach 
manifolds, 

^:cr=(7,U)'— >(Rij-2aiUajU+2p,AU-M) 

from  X  into  a  space  of  symmetric  covariant  tensors  and 
functions  which  will  be  isomorphic  to  some  suitable  weighted 
Sobolev  spaces,  the  solution  set  of  Einstein’s  equations  with  the 

desired  asymptotic  and  differentiability  conditions  is  just  ^  *(0). 

The  coordinate  freedom  can  then  be  handled  by  the  fact  that 
the  set  of  diffeomorphisms  on  X  that  approach  the  identity  at 
infinity  operate  in  a  natural  way  by  pull-back  on  X  and  that 
the  corresponding  orbits  are  submanifolds  of  X  .  If  we  could 

show  that  the  set  £  *(0)  is  contained  in  the  orbit  through  the 
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spherically  symmetric  solution  we  would  have  proved  that  these 
solutions  are  obtained  from  the  spherically  symmetric  one  by 
diffeomorphisms,  or  in  other  words,  that  the  set  of  classes  of 
physically  equivalent  solutions  consists  only  of  the  spherically 
symmetric  one. 

We  do  not  prove  this  general  a  result,  however,  for  several 
reasons.  First,  the  tangent  map  £'  of  jzf  is  not  surjective  so, 

as  wTe  will  see,  we  are  unable  to  prove  that  X  !(0)  is  a 
submanifold.  We  are  thus  unable  to  tackle  the  problem  of 
existence  of  solutions  and  in  the  stationary  case  must  assume 
that  rigid,  rotating  perfect  solutions  with  nonzero  angular 
momentum  exist.  Secondly,  in  order  to  extend  our  analysis  of 
the  linearized  equations  on  the  spherical  background  to 
backgrounds  away  from  spherical  symmetry  we  must  use  a 
slicing  theorem  for  the  orbits  of  the  dif feomorphism  group  to 
obtain  an  elliptic  operator  from  the  tangent  map  £’  to  which 
we  can  apply  a  theorem  of  Nirenberg  and  Walker  (1973).  The 
slice  y  obtained  represents  the  set  of  equivalence  classes 
locally  and  Nirenberg  and  Walker's  theorem  deals  with  operators 
which  are  dif ferentiably  close  to  some  operator  so  both 
theorems  force  us  to  consider  only  solutions  which  are  close,  in 
the  Banach  topology  sense,  to  the  spherically  symmetric- 
solution  Furthermore,  since  Nirenberg  and  Walker’s  theorem  is 
a  proof  by  contradiction  there  is  little  hope  of  obtaining  an 
upper  bound  on  how  far  away  from  spherical  symmetry  our 
result  is  still  valid  and  thus  little  hope  in  finding  an  upper 
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bound  on  J,  in  the  stationary  case,  for  which  the  uniqueness 
result  we  obtain  holds.  There  is  little  doubt  that  such  an  upper 
bound  exists,  though,  since  we  expect  branch  points  as  in  the 
classical  Maclaurin-Jacobi-ellipsoid  series. 

With  all  these  qualifiers,  what  can  we  actually  show?  Since 
the  slice  represents  the  physically  distinct  spacetimes  (in  a 
one  — to  — one  correspondence)  we  still  have  a  physically 
interesting  result  if  we  restrict  consideration  to  it.  In  the 
static  case  we  show  that  there  are  no  non-constant  C1  curves 
of  solutions  in  going  through  the  spherically  symmetric 

solution  in  £  (0)  which  have  the  same  mass  and  Tb.  (We 
already  consider  the  equation  of  state  to  be  fixed,  the  mass  has 
been  fixed  by  the  choice  of  Banach  space  and  Tb  is  just  a 
constant  -which  we  will  always  consider  to  be  fixed.)  In  the 

stationary  case  we  show  that  there  are  no  non  — constant  C1 

rv 

curves  of  solutions  in  tf  (the  corresponding  slice)  which  have  a 
constant  angular  momentum  J  (and  where  again.  p( p),  m  and 
Tb  have  been  fixed). 

The  spherical  solutions  are  clearly  very  important  in  our 
framework  so  in  the  next  section  wTe  will  look  at  them  We  also 

need  to  investigate  the  equations  of  state  which  are  physically 
appropriate  and  determine  what  differentiability  conditions  they 
impose  before  beginning  our  study  of  the  weighted  Sobolev 


spaces. 


2.3  Spherically  symmetric  solutions 

As  we  have  already  noted,  when  the  strength  of  the 
gravitational  field  ¥  is  a  function  only  of  U  then  U  has  only 
one  (nondegenerate)  critical  point  at  the  center  x  ,  Uc=U(xc),  E 

3 

is  dif feomorphic  to  tR  and  the  spacetime  is  spherically 
symmetric.  If  M  is  spherically  symmetric  then  W  is  clearly  a 
function  only  of  U  so  U  has  only  one  critical  point  unless  M 
is  flat.  The  stability  of  Morse  functions  with  a  finite  number  of 
critical  points  (c-.f.  Golubitsky  and  Guillemin  (1973),  p.72,79) 
indicates  that  the  property  of  isolated  nondegenerate  critical 
points  is  an  open  one.  Since  as  we  have  seen  above  we  are 
forced  to  restrict  to  solutions  which  are  close  to  the  spherically 

symmetric  solution  (in  at  least  a  C1  sense  of  close  in  the 
function  space  in  which  U  lies)  we  can  always  consider  U  to 
have  only  one  (nondegenerate)  critical  point  and  can  take  E  to 

3 

be  dif  feomorphic  to  so  standard  Euclidean  coordinates  can 
be  used. 

The  following  results  about  spherically  symmetric 

spacetimes  are  standard  (c.f.  Kunzle  (1971))  The  line  element  of 
E  can  be  written  as 

ds2=¥~2(U)dU2+r2(d6T+sin2edy2)  (2.1) 

where  6  and  are  angular  coordinates,  r  (U):=2/R  and  the 
static  field  equations  (1.59-1.62)  imply 
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rd(W2)+4W2dr=2MrdlI  (2.2) 

Wdr=(l+r2W2+pr2)1/2dU  (2.3) 

These  can  be  integrated  to  give 

W  =  m(U)/r2  (2.4) 

where 

m2(U)=2  Ju  Mr4dU 

and  m(U)  agrees  with  the  total  gravitational  mass  rn  in  the 
exterior  region.  In  the  vacuum  region  (2.3)  implies 

U=— sinh  1(m/r)=— log  (m/r+(l+m2/r2)1,  2).  (2.6) 

Note  that  r  is  not  the  Sch warzschild  radial  coordinate  r'.  In 
this  coordinate 


U  =  ^log  (1  — 2m/ r' 


(2.7) 


in  vacuum.  We  could  also  write  the  line  element  a; 


,  2  .  2,  ,  ,  2  2/  ,  _ 2  .  2_  ,  2n 
ds  =A  (r)dr  +r  (d 6  +sm  ddcp  , 


(2.8) 


—1  2  2—1/2 
where  A=W  dU/dr  so  that  in  vacuum  A  =  (l  +  m  /r  )  *'  or  in 

terms  of  a  Riemannian  normal  coordinate  system  (x  )  based  at 
the  centre  which  are  related  to  the  polar  coordinates  (r ,6, if) 
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in  the  usual  way,  i.e. 


ds2=(6ij— B(r)xiXj)dx1dx^ 


(2.9) 


_ 2  2 

Here  B  =  r  (1-A  )  which  in  vacuum  is 


T3  2  -4,,  2  ,  2,-1 

B=m  r  (1  +  m  /r  )  . 


(2.10) 


Having  chosen  a  function  p( p)  (subject  to  conditions 
which  we  will  state  in  the  next  section)  the  equation  of 
hydrostatic  equilibrium  dp  +  (p  +  p)dTJ  =  0  determines  p(U)  (and 
hence  p(U),  p(U),  M(U)  )  uniquely  up  to  a  constant  which  can, 
for  example,  be  chosen  to  be  the  value  Ub  of  U  on  the  star 
boundary,  i.e. 


(2.11) 


Since  U,  dU/dr  and  p  must,  be  continuous,  as  is  seen  from 
Einstein’s  equations  (1.29,  1.30)  with  the  assumption  that  p(p)  is 


piecewise  C°,  it  follows  that  (2.2)  and  (2.3)  can  be  integrated  for 


r(U)  and  W(U)  in  terms  of  their  values  rb=-m/(sinh  Ub)  and 


-2 

Wb=mrb  so  that  the  whole  solution  is  determined  by  the  two 


constants  m  and  Ub.  (A  subscript  b  will  always  refer  to  the 
value  of  a  scalar  on  the  surface  of  the  star,  even  in  the 

stationar}^  case.)  The  center  is  then  defined  by  that  value,  Uc, 
of  U  for  which  r  vanishes.  But  since  it  is  a  critical  point  of  U 
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we  must  also  have  Wc=0  which  is  likely  to  impose  an 
additional  condition,  determining  U  ^  in  terms  of  ni . 

At  least,  this  is  what  happens  for  specific  very  simple 
models,  like  the  interior  Sc-h warzschild  solution  (where  p{p)  =  Po 
=const.  ).  We  conjecture  therefore  that  to  given  p(p)  and  m 
there  exists  a  unique  spherically  symmetric  static  stellar  model 
provided  p{ p)  satisfies  the  conditions  given  in  the  next  section 

and  0^m^mcrit.  (Clearly  m  must  be  bounded  from  above  to 
prevent  the  formation  of  a  black  hole,  but  as  we  noted  in  §1.1 

the  value  of  mcrit  depends  on  the  equation  of  state.)  A  proof 
along  the  above  lines  would  have  to  depend  on  a  very  close 
analysis  of  the  nonlinear  system  {2.2,  2.3)  and  may  not  be  so 
easy  to  prove  for  a  general  equation  of  state.  The  theorem 
which  we  prove  in  chapter  4  is  a  linearized  version  of  this 
statement.. 

If  this  conjecture  is  true  then  we  have  for  every  p{ p)  a 

curve  m1— wrhich  tends  to  the  flat  solution  (7, ,=  <5. ,,U=0)  for 
m — >0  .  This  curve  will  be  continuous  when  we  topologize  X  . 
Unfortunately,  it  fails  to  be  differentiable  at  m=0  .  This  fact 
prevents  us  from  applying  to  this  situation  a  powerful 
technique  of  Cantor  (1979)  which  requires  some  elliptic  operators 

9  0  9  0 

to  be  dif f erent iably  close  to  ones  with  constant  coefficients 
(This  theorem  is  theorem  3.7  of  the  next  chapter.  We  will 
discuss  this  more  later  when  we  have  developed  the  tools  to  do 

so.) 
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2.4  Equation  of  state 

We  now  want  to  consider  the  conditions  on  the  equation  of 

state  p:[0,pc]  p  1  >p( p)  which  are  sufficient  to  prove  our 

uniqueness  theorems  for  stellar  model  configurations  with  fixed 

m,  Tb,  p(p)  and  J  in  the  static  and  stationary  cases.  In 
addition  to  the  physically  obvious  conditions  0^  p  <coj  p  <°° 
made  in  the  thermodynamic  argument,  we  also  make  the 
standard  assumption  that  p^p(p).  This  is  implied  by  Hawking’s 
dominant  energy  condition,  namely,  for  every  timelike  vector 

field  Ta|5 'xaXp^Q  and  is  a  non— spacelike  vector.  This 

can  be  interpreted  as  saying  that  to  all  observers  the  energy 
density  appears  non  — negative  and  the  local  energy  flow  vector  is 
non  — spacelike  (Hawking  and  Ellis  (1973)  p.91).  Since  dp/dp  can 
be  interpreted  as  the  inverse  square  of  the  (adiabatic)  velocity 
of  sound  in  the  Newtonian  limit  we  assume  also  that  dp/dp^O. 

We  are  only  interested  in  non-singular  solutions  so 
equation  (2.11)  and  the  corresponding  equation  obtained  from  the 
stationary  hydrodynamic  equilibrium  equation  (l.ll), 

log  (T/Tb)= jbptpHpkdp,  (2.12) 

imply  that  the  integral  on  the  right  hand  side  must  be  finite 
for  all  pE[0,pJ  which  is  true  iff  limp/p  =  0,  re.  the  highly 

r  L  CJ  p - 

relativistic  limit  p= p  cannot  occur  at  extremely  low  densities, 
as  one  expects.  This  implies  that  there  exists  e,  d  with  0^£<1, 


h 
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0<d<°o  such  that  lim  p  p  e=d 

p  *o 

Let  us  first  look  at  situations  in  which  the  matter 
distribution  is  not  compact.  We  do  this  for  the  spherically 
symmetric  case  as  we  are  always  interested  in  solutions  close  to 
it.  Assume  there  is  no  compact  boundary  hypersurface  in  E 

given  by  p  =  0  so  that  p  has  an  asymptotic  expansion  p  =  0(r  k) 
at  oo.  Using  the  asymptotic  boundary  conditions  (1.33)  and  the 
equation  of  hydrostatic  equilibrium  (1.32)  we  find  that 

M  =  0(r  F/  [l  £/)  at  oo.  Together  with  Einstein’s  equation  AU=M 
this  implies  that  £^4/5  .  (If  e=l/(t+ 1)  with  4  then 

p=(mdr  ^E+l)  1)  !+0(r  1  2)  at  infinity.)  Since  the  physically 
reasonable  stellar  models  can  be  considered  to  have  compact 
matter  distributions  and  since  in  our  proofs  wre  use  the  exact 
solutions  for  vacuum  obtained  in  the  previous  section,  we  will 
have  to  make  the  assumption  that  equations  of  state  such  that 

£^4/5  must  be  of  a  form  which  results  in  a  compact  matter 
distribution. 

Now  let  us  look  at  the  star  boundary  conditions,  where 
p  =  0.  If  £=0  p(0)>0  so  p  is  discontinuous  at.  the  star 

boundary.  (We  will  always  consider  p  and  p  to  be  identically 
zero  in  the  vacuum  region.)  In  order  to  investigate  the  junction 
conditions  further  let  us  assume  that  (l- l)/t<e^t/{t  +  l)  for  some 

1^0  (noting  that  since  £^0  £  =  0  if  i  —  0  )  and  that  p( p)  is  C 

piecewise  on  (0,pc]  when  £  is  in  this  range.  Since  U  has 
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only  one  critical  point  and  even  in  the  stationary  case  we 
remain  close  to  the  spherical  solution  T  also  has  only  one 
critical  point  in  its  domain  which  is  the  interior  region.  From 
(2.12)  we  see  that  in  the  stationary  case  p  is  determined  as  a 
function  of  T  since  we  consider  Tb  given  and  fixed.  Thus  the 
continuity  of  derivatives  of  the  mass  density  across  the 
boundary  of  the  star  can  be  investigated  using  U  (  T  )  in  the 
static  (stationary)  case  to  calculate  the  limit  of  these 
derivatives  as  the  star  boundary  is  approached  from  the 
interior.  If  we  assumed  as  much  differentiability  for  p(p)  as 

being  piecewise  C1  on  [0,pc]  then  the  right  hand  side  of  the 
equation  dp/dlogT=p(p)+p  would  be  Lipschitz  in  p  on  closed 
intervals  of  [0,pc]  and  a  standard  existence  and  uniqueness 
theorem  for  first  order  differential  equations  (c.f.  Hartman 
(1973))  would  imply  that  p(0)>0  (  e  =  0  )  in  order  for  p  not  to 
vanish  in  the  interior  region.  Thus  when  p(0)=0  we  must 

always  have  that  limdp/dp  =  °o,  which  is  the  reason  for  the 

p 

differentiability  being  specified  only  on  (O.pJ  The  amount  of 
differentiability  demanded  in  this  range  is  purely  for 
mathematical  reasons;  in  order  to  have  the  star  boundary 
conditions  determine  the  proper  weighted  Sobolev  space  there 
cannot  be  less  differentiability  in  the  interior  region  than  there 
is  at  the  star  boundary.  This  restriction  would  be  very  difficult 
to  remove  in  our  formalism  but  it  is  not  a  very  serious  one  as 
all  equations  of  state  could  be  approximated  arbitrarily  closely 
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by  a  C  equation  of  state. 

_ g  p 

Let  us  also  assume  that  p  p  is  C  in  [0,5)  for  some 
small  6  so  that 

lim(d’p/dp1)p1  c~dE--(E— i+1)  for  O^i^f. 
p— >o 

We  look  first  at  the  junction  conditions  in  the  static  case.  Let 

CO 

z  be  a  C  coordinate  defined  in  a  neighborhood  of  the  star 
boundary  such  that  z=l  corresponds  to  the  p  =  0  hypersurface 
and  p(z)>0  for  z<l.  The  differentiability  of  m  (and  p  and  p) 
can  then  be  investigated  using  the  relation 

dfM/dz^=  E  £  f!(m1!-ms!Il(i+ms_i+1)) 

s— 1  mp — hm8+s=f  “  i=l 

•  dsM/dUs(((dU /dz)d/dz)midU/dz)"'(((dU/dz)d/  dz)msdU/ dz) 

This,  together  with  the  following  argument,  will  show  that  the 

limit  as  z — >1~~  of  d?M/dz*  is  given,  for  (t-l)/l<c^l/(l+\.)  by 

(dU/dz|z=c)C  times  the  limit  as  p — >0^  of  the  following 
expression. 

d*M/dUe=(-l)^e-2U  £  21[((p+p)d/dp)J(p  +  p)+2((p  +  p)d/dp)^1(p+p)] 

i+  }=t 

A  long  calculation  yields 

limp1  ^  1^(d1M/dU  )=(  —  1)  ye  bd  (i— 1)) 

p-^0 
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for  and  (l  —  l)/i<£^i/(i  +  l).  Recall  that  lim  dkM/dz  =0 

z — >0  + 

for  all  k  since  M=0  in  the  exterior  region.  Thus  M  (and  p  ) 
l—l  l 

are  C  (piecewise  C  if  e  =  t/(t  +  l))  across  the  star  boundary. 

We  are  justified  in  using  U  as  a  coordinate  like  this  since  we 
know  it  is  at  least  C1  so  we  can  calculate  the  first  derivative  of 

M  and  find  that  M  is  C°  (if  t> 0  ).  Einstein’s  equations  then 

2  3 

imply  that  U  can  then  be  taken  to  be  C  and  7  to  be  C  so 
we  can  take  the  second  derivative  of  M,  and  so  on.  Thus  for 

(l— l)/t<E^l/{l+l)  U  is  C*+1,  7  is  C(+2  and  if  E  =  t/(l  + 1)  then  U 

1+2  {+3 

is  piecewuse  C  and  7  is  piecewise  C  The  equation  of 
hydrostatic  equilibrium  shows  that,  p  is  as  differentiable  as  U. 

In  the  stationary  case  a  completely  analogous  derivation 
may  be  done  using  log  T  in  place  of  U,  resulting  in,  for  0^i^£, 

limp(i+1)c'ie2Ud1M/d(log  T)i=^di+1a--(ic-(i-l)) 

p  >0 

Using  Einstein’s  equations  (1.22-1.24)  the  same  type  of  argument 
shows  that  U  and  7  have  the  same  differentiability  as  in  the 

static  case  and  h  is  C  ,  piecewise  C  if  £=l/(t+ 1). 

Our  assumptions  on  the  equation  of  state  are  then, 

(i)  p(p)^0  on  some  interval  I=[0,pJ  , 

(ii)  dp/dp^0  on  1  , 

(iii)  the  solution  is  nonsingular,  i.e.  there  exists  e,  d  with 
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0^£<1,  0<d<oo  such  that  limp  Ep  =  d  , 

p-^O 

(iv)  For  with  e=0  if  1=0,  p(  p)  is  C? 

piecewise  C*  1  on  (0,pc]  and  p  cp  is  CC  in  [0,6)  for  some  small 
<5>0, 

(v)  if  a^4/5  the  equation  of  state  is  such  that  the  finite  mass 

spherical  solution  with  this  equation  of  state  has  a  compact 
matter  distribution. 

The  results  obtained  in  this  section  are  summarized  in  the 
following  theorem.  The  case  £=£  =  0  was  given  in  Klinzle  and 
Savage  (1980b). 

Theorem  2.1:  If  the  equation  of  state  satisfies  (i)  to 
(iv)  above  then  the  corresponding  stellar  model 
consists  of  a  compact  region  DCE  with  p^p^O  in  D 
and  an  exterior  vacuum  region.  For  all  compact 
matter  distributions  the  differentiability  of  M  (and 
p)  across  the  boundary  3D  of  D,  defined  by  p  =  0,  is 
determined  by  £  and,  for  (l  —  1) / / (l  +  l)  with  l> 0 

M  is  Ce_1  (and  piecewise  if  s  =  t/(t+l)  )  while  M  is 
piecewise  C°  if  s  =  (= 0,  with  the  £th  derivative  normal 

to  0D  going  like  const  p  as  p 

(We  are  now'  including  equations  of  state  of  the  form 
p  =  Apq  for  0<q<l,  the  finite  stellar  models  in  which  p  tends  to 
0  on  the  boundary,  wrhich  were  excluded  for  mathematical 
convenience  in  Kunzle  and  Savage  (1980b).) 
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Let  us  then  denote  by  Sp(p)  (  Sp(p)  ),  or,  for  short  Sp  (  Sp  ), 
the  set  of  static  (stationary,  axisy mmeti  ic)  solutions  of 
Einstein’s  equation  with  a  fixed  equation  of  state  p( p) 
satisfying  conditions  (i)  to  (v)  above,  with  a  fixed  surface 
temperature  Tb  and  such  that  U  has  only  one  critical  point. 

Let  (  2p>m  )  denote  the  corresponding  restriction  to  solutions 

with  a  given  constant  mass  m  .  Eventually  one  would  like  to 

prove  that  Sp  consists  only  of  the  spherically  symmetric 

solutions  and  that  Spm  consists  of  a  family  of  solutions  which 
close  to  the  spherical  solution  is  parameterized  by  the  angular 
momentum  J  (and  for  larger  J  perhaps  bifurcate  m  a  similar 
manner  as  the  Maclaurin-  Jacobi-ellipsoid  sequence  in 
Newtonian  theory).  As  we  have  already  noted  we  do  not  know  a 
priori  how  large  these  sets  are  or  what  kind  of  topology  and 
differentiable  structure  they  can  be  given,  but  they  can  be 
regarded  as  the  inverse  image  of  a  differentiable  map  on  a 
bigger  set  that  can  be  provided  with  a  fairly  natural  Banach 
manifold  structure. 

Unfortunately  there  is  considerable  arbitrariness  in  the 
choice  of  the  Banach  manifold  structure  for  a  set  of  tensor 

fields  on  a  non-compact  manifold.  We  try  to  make  the  most 
appropriate  choice  among  the  weighted  Sobolev  spaces  which  are 
presented  along  with  some  general  Banach  manifold  theory  in 

the  next  chapter.  In  chapters  4  and  5  we  will  see  that  physical 
and  mathematical  considerations  combine  to  limit  this  choice 
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remarkably.  However  one  cannot  exclude  the  possibility  that  a 
manifold  of  different  data  and  field  equations  in  a  different 
form  might  perhaps  lead  to  stronger  results  than  we  get 


CHAPTER  111 


MATHEMATICAL  PRELIMINARIES 

3.1  Differential  geometry  in  relativity 

The  use  of  differential  topology  and  differential  geometry 
techniques  in  relativity  theory  is  certainly  not  new.  For 
example,  the  incompleteness  theorems  of  Hawking  and  Penrose 
from  which  the  existence  of  black  holes  may  be  inferred  — 
under  reasonable  mathematical  assumptions  (c.f.  Hawking  and 
Ellis  (1973)) —  and  the  geometric  analysis  of  spatial  and  null 
infinity  by  conformal  mappings  use  techniques  from  the  study 
of  the  topology  and  geometry  of  finite  dimensional  manifolds. 

IFheeter  fl962)  was  the  first  to  point  out  the  relevance  of 
infinite  dimensional  manifold  theory  to  relativity  with  the 
introduction  of  superspace  S,  the  quotient  space  of  riemannian 
metrics  on  a  given  three  dimensional  manifold  obtained  by 
identifying  metrics  which  can  be  obtained  one  from  another  by 
a  coordinate  transformation  (  S  is  a  metric  space  but  does  not 
have  a  manifold  structure,  c.f.  Fischer  (1970).)  The  universe  can 
then  be  viewed  as  an  evolving  geometry  and  thus  as  a  curve  in 
S,  allowing  a  dynamical  theory  of  relativity.  This  application  is 

regarded  as  soft  in  that  infinite  dimensional  manifolds  are 
involved  mostly  as  a  language  convenience  and  as  a  guide  to  the 
theory’s  structure.  Brill  and  Deser  (1968)  gave  the  first 
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substantial  hard  theorem  using  infinite  dimensional  analysis, 
albeit  in  an  informal  way,  when  they  showed  that  any 
non  — trivial  perturbation  of  Minkowski  space  leads  to  a 
spacetime  with  strictly  positive  total  gravitational  mass.  Their 
technique  is  to  show  that  on  the  space  of  solutions  to  Einstein’s 
equations  the  mass  function  has  a  non-degenerate  critical  point 
at  Minkowski  space.  Later  investigations  of  the  positive  mass 
conjecture  have  used  both  infinite  dimensional  manifold  theory 
(Choquet— Bruhat  et  al.  (1979))  and  finite  dimensional  manifold 
theory  (Schoen  and  Yau  (1979)  and  Jang  (1979)). 

Our  approach  also  uses  infinite  dimensional  manifold 
theory  and  is  based  on  the  problem  of  linearization  stability 
which  has  been  investigated  b}r  Fischer  and  Marsden  (1975a, 
1975b).  The  reason  for  this  is  as  follows,  where  we  will  restrict 
to  the  static  case  as  the  reasoning  is  the  same  in  the  stationary 

case.  We  can  regard  Sprr_  as  the  inverse  image  cf  0  of  a 
differentiable  map  £:X — >Y,  obtained  from  Einstein's  equations, 
between  suitable  Banach  spaces  X  and  7  where  a  point  atX 
will  uniquely  specify  a  spacetime,  he  also  knew  that  there  is  a 

unique  spherically  symmetric  solution  crE^  (0)  (by  including  a 
specification  of  Ub  in  the  space  X  if  necessary)  and  we  want  to 

answer  the  question  what  does  the  solution  set  £  (0)  look  like0 

To  answer  this  question  globally  is  a  formidable  problem  and 
has  been  done  only  in  the  special  case  where  p  =  const.  by  a 
similar  approach  to  that  taken  in  the  black  hole  theorems 


. 
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(Lindblom  (1980)),  as  we  have  noted.  In  order  to  have  a  more 
general  approach  (and  the  generality  is  proven  by  its 
applicability  to  the  stationary  case  as  well)  we  attempt  to 
answer  the  question  locally,  near  <7,  by  considering  solutions  in 

Spm  which  are  close  to  cr  in  X.  However  £  is  formed  from 
Einstein’s  equations  and  is  thus  highly  non-linear  so  a  direct 
approach  is  unfeasible.  The  usual  approach  in  such  situations  is 
to  linearize  the  equations,  solve  these  linearized  equations  and 
assert  that  they  are  an  approximation  to  the  true  solution  of 
the  non-linear  equations.  More  explicitly,  for  o  near  <j  write 

<t(A)  for  a  parameter  A  and  expand  as  a(A)=(j+A(J1+A  <J2+‘"  ■  The 
approximation  to  first  order  is  cr+Acv  where  a1=(da/d A)|x=0. 
Demanding  that  o( A)Cdzf  ^0),  i.e.  £{_o{ A))  =  0,  we  find  that 

^'(cr)ai=0  where  ^'(cr):TaX — is  the  tangent  map. 

The  implicit  assumption  that  the  solution  to  the  linearized 
equations  is  an  approximation  to  the  full  equations  is,  however, 
not  always  valid.  In  fact  Fischer  and  Marsden  (1975a)  have 

shown  that  if  the  universe  is  toroidal,  T3X®,  where  TJ  denotes 
the  flat  3-torus,  this  assumption  is  not  valid  for  the  flat  space 

Einstein  equations  R„M(g)  =  0  .  For  a  more  visible  example  of  such 

2  2  2  —  1 

e  situation  consider  f:&  — *&,  (x.y)1  -y  .  Then  f  (0) 

=i(x,y)ly=±xi  but,  f'(0,0)h=0  for  all  h=(h1,h2)eT(0>0)«2=?R2  Only 

by  going  to  the  second  order  condition  f  (0,0)(h,h)=2(h.2—  hg)— 0 
<=>h2=±h1  do  we  get  a  true  approximation.  Situations  in  which 
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this  assumption  is  valid  will  be  called  linearization  stable. 
Definition  3.1:  Let  E  and  F  be  topological  vector 
spaces  and  L:E  >F  a  differentiable  mapping.  We  say 
L  is  linearization  stable  at  x0<EE  iff  for  every  h<EE 
such  that  L'(x0)h=0  there  exists  a  differentiable  curve 
x(t)£E  with  x(O)=x0,  L(x(t))  =  L(x0)  and 

x'(0)(:=dx/ dt|t=0)=h. 

The  vacuum  Einstein  equations  have  been  well  studied  in 
this  manner  by  considering  the  induced  metric  g  and  the  second 

fundamental  form  fi  on  a  spacelike  hypersurface  of  the 
spacetime  as  initial  data,  which  satisfy  some  nonlinear 
constraint  equations,  for  some  evolution  equations.  By  using 
various  elliptic  operator  techniques  Fischer  and  Marsden  (1975a) 
have  proven  the  following  theorem. 

Theorem  3.1:  Suppose  there  is  such  a.  spacelike 

hypersurface  N  of  the  spacetime  with  g  and  0 
satisfying  (l)  there  are  no  infinitesimal  isometries  7i 
on  both  g  and  fi  (with  7)  vanishing  at  infinity  if  N 
is  not  compact),  (ii)  if  Q=0  and  N  is  compact  then  g 
is  not  flat,  (iii)  if  D^O,  tr(fi)=  trace  of  Q  is  constant 
on  N  if  N  is  compact  and  tr(Q)=0  if  N  is 
noncompact,  (iv)  if  N  is  noncompact,  g  is  complete 

and  asymptotically  flat.  Then  near  N,  B/Ul,(g)=:0  1S 
linearization  stable. 

(See  Choquet-Bruhat  et  al.  (1977a)  for  a  more  recent, 


. 
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simpler  proof  based  on  some  results  of  weighted  Sobolev  spaces.) 

The  toroidal  example  fails  because  condition  (ii)  fails  for  N  =  T3  . 
The  corollary  that  in  Minkowski  space  the  vacuum  Einstein 
equations  are  linearization  stable  was  obtained  independently  by 
Choquet  —  Bruhat  and  Desei  (1973).  (In  this  case  the  linearized 
equations  are  the  weak  field  approximations  used  to  study 
gravitational  waves.) 

We  will  be  following  the  spirit  and  not  the  details  of  this 
work  on  linearization  stability.  In  fact,  as  we  will  see,  our 
tangent  map  is  not  surjective  so  we  cannot  show  that  our 
equations  are  linearization  stable.  However,  the  differential 
structure  which  we  choose  for  the  Banach  space  X  will  allow  a 
slicing  for  the  action  of  the  dif f eomorphism  group  on  X  neai  cr 
which  together  with  some  theorems  on  elliptic  operators  enables 
us  to  prove  a  theorem  about  the  physically  distinct  solutions  in 

^(0)  near  cr.  X  must  clearly  be  an  infinite  dimensional  space 
so  in  the  next  section  we  take  a  look  at  infinite  dimensional 
Banach  spaces.  The  following  section  will  then  present  the 
weighted  Sobolev  spaces  introduced  by  Cantor  which  we  use  to 
model  X  on  as  well  as  two  theorems  due  to  Cantor  (1979)  and 
Nirenberg  and  Walker  (1973)  for  elliptic  operators  which  are 

dif f erentiably  close  to  elliptic  operators  with  constant 

coefficients  which  we  use  to  extend  our  analysis  of  ^'{o)  to  a 
neighborhood  of  cr  in  X  .  It  is  these  latter  two  theorems  which 
together  with  the  slicing  theorem  allow  us  to  prove  our  results 
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about  curves  of  physically  distinct  solutions,  even 
have  not  shown  linearization  stability  (so  that  we 
existence). 


though  we 
cannot  prove 


3.2  Banach  spaces  and  infinite  dimensional  manifolds 

A  Banach  space  E  is  a  complete  normed  vector  space, 
while  a  Banach  manifold  N  is  a  (Hausdorff)  topological  space 
with  a  maximal  atlas  of  charts  )(<,<Ja,Ua)(  such  that  each  tpa. 

Ua  >VaCE  where  Va  is  open  in  E,  is  a  homeomorphism.  In 
particular  a  Banach  space  or  an  open  subset  of  a  Banach  space 
is  a  Banach  manifold.  SCN  is  a  submanifold  of  N  if  we  can 

write  E=F©G  (topological  sum)  and  for  every  xES  there  is  a 
chart  <^:UCN— >VCE  of  N  where  xQJ  such  that 
$*(UnS)=Vn(FXjw})  where  wEG  .  In  other  words,  the  chart  <p 


flattens  out  S  making  it  lie  in  the 
closed  subspace  F  of  E  is  such  that 
another  closed  subspace  of  E,  we  say 


subspace  F . 

E  =  E®G  where 
that  F  splits 


When  a 
G  is 

(c.f.  Marsden 


(1974)). 


Such  splittings  clearly  always  occur  in  the  case  where  E 
is  finite  dimensional  or  a  Hilbert  space  (where  the  inner  product 
can  be  used  to  find  the  orthogonal  complement  to  F  ).  However, 


intuition  from  finite  dimensional  spaces  sometimes  fails  in 

infinite  dimensional  spaces.  For  example,  in  iRn  a  closed, 
bounded  subset  with  non-empty  interior  is  always  compact 
whereas  in  an  infinite  dimensional  Banach  space  it  is  never 
compact  under  the  norm  topology.  More  importantly  for  our 
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work,  as  we  will  soon  see,  an  arbitrary  closed  subspace  F  of  an 
infinite  dimensional  Banach  space  does  not  always  split. 
Fortunately,  though, the  generalization  from  differential  calculus 

on  iR  to  differential  calculus  on  Banach  spaces  is  remarkably 
smooth  with  many  theorems,  including  the  inverse  and  implicit 
function  theorems,  similarly  phrased  (c.f.  Choquet-Br uhat  et  al. 
(1977b)).  These  latter  two  theorems  are  important  for  obtaining 
many  substantial  results  in  Banach  space  theory,  seme  of  which 
are  the  following. 

Let  N  and  R  be  Banach  manifolds  and  f:N — >R  a  C’1  map. 

W e  are  interested  in  when  S=f  ^Jq)  is  a  submanifold  of  N  and 

when  f  is  linearization  stable  on  S  (i.e.  we  would  like  £  1(0)CX 
to  be  a  submanifold  with  its  tangent  space  at  cr  determined  by 
£'(g)~ 0  for  a  in  a  neighborhood  of  a  ).  Suitable  conditions  are 
given  by  the  following  theorem  (c.f.  Fischer  and  Marsden  (1975a)) 

Theorem  3.2:  Let  x0€!N  and  f(x0)=y0  .  Suppose  that 
f'(x0)  is  surjective  and  that  kernel  f'(x0)  splits  Then 

j  1. 

f  (yc)  is  a  C  submanifold  near  x0  and  f(x)  =  y0  is 
linearization  stable  about  x0. 

Proof:  Work  in  a  chart  UCE  for  N  and  write  E  =  E1©E2 
where  E2=ker  f'(x0).  Consider  the  map  <f>,  defined  near  x0  to 
EjXF  (where  R  is  modelled  on  the  Banach  space  F  )  by  <f*(x1,x2) 

=(x1,f(x1,x2))  .  Since  (f|WxE2)'(xo)  is  an  isomorphism,  <F(x0)  is  an 
isomorphism  so  the  inverse  function  theorem  implies  that  ^  is 
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a  local  diffeomorphism.  Thus  <t>  1  gives  a  chart  for  f  \y0)  near 

x0  modelled  on  E{  and  Tx  f_1(y0)=ker  f'(xn)  .  Now  h€Tv  N  is  a 

first  order  deformation  iff  hGker  f'(x0)  and  since  f_1(y0)  is  a 

submanifold,  there  exists  a  curve  x(X)£f_1(y0)  which  is  actually 
tangent  to  h  .  Thus  f  is  linearization  stable  about  x0  R 

There  are  thus  two  properties  of  ^'(cr^T^X — »Tj£^Y  in 
which  we  are  interested,  surjectivity  and  the  splitting  of  its 
kernel.  Clearly  surjectivity  depends  on  the  image  space  Y  — it 
must  be  large  enough  but  not  too  large.  We  will  see  later  that 
in  our  case  it  is  too  large  but  it  is  by  no  means  obvious  how  to 
pick  a  submanifold  of  Y  which  will  be  the  right  size  The 
splitting  of  the  kernel  is  not  as  formidable.  In  fact  much  work 
has  been  done  on  orthogonal  decompositions  of  symmetric 
covariant  tensor  fields  S2  (of  which  the  metrics  are  an  open 
subset)  over  both  compact  and  noncompact  manifolds  (Berger 
and  Ebin  (1969),  York  (1974),  Cantor  (1979)).  In  trie  case  of  a 
noncompact  manifold  though,  the  decomposition  must  have  some 
dependence  on  the  space  S2,  such  as  its  asymptotic  properties 
In  the  next  section  we  will  see  see  some  spaces  which  are 
suitable  for  modelling  the  asymptotically  flat  metrics  on  Y,  but 
there  will  be  a  problem  with  the  rate  of  falloff  at  infinity 
required  for  such  a  splitting.  This  can  fortunately  be  resolved. 

But  why  are  we  interested  in  this  splitting  property  if  we 
can  not  show  linearization  stability  anyway?  The  answer  is  that 
this  sort  of  splitting  together  with  some  properties  of  weighted 
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Sobolev  spaces  (and  in  the  stationary  case  another  splitting 
theorem  for  vector  fields  as  well)  enables  one  to  show  that  near 
cr  the  orbits  due  to  the  action  of  the  diffeomorphism  group  are 

stacked  in  the  sense  that  in  a  neighborhood  V  of  cr  there  is 
a  submanifold  t/,  called  a  slice,  which  passes  through  cr  and 
such  that  all  the  orbits  in  this  neighborhood  V  pass  through  cr 
only  once.  Thus  represents  the  physically  distinct 

9  9  99 

spacetimes  which  are  close  to  cr  .  Thus  given  any  physical 

curve  of  solutions  in  £  ( 0)0 if  the  behaviour  of  the  tangent  to 
such  a  curve  near  <j  can  be  investigated  with  the  use  of  the 
theorems  about  elliptic  operators.  By  choosing  the  space  X 
well,  whether  there  is  such  a  nonconstant  curve  will  tell  us  how 
unique  the  spherically  symmetric  solution  is,  at  least  locally,  in 
a  similar  manner  to  that  of  the  stationary  black  hole  theorem 

Before  investigating  spaces  to  model  X  on,  we  state  a 
similar  result  to  that  of  the  submersion  result  of  theorem  3.2 

which  we  use  in  some  later  proofs. 

Theorem  3.3:  Let  f.N — >R  be  injective  and  closed.  If 

f'(x):TxN — >Tf(-x)R  is  injective  and  its  image  splits  for 
each  x6N  then  f(N)  is  a  submanifold  of  R. 

Proof:  (c.f.  Lang  (1972)  p.27) 

3.3  Weighted  Sobolev  spaces 

The  definitions  and  theorems  (except  theorem  3.14)  in  this 
section  come  from  Cantor’s  (1975a, 1975b, 1979)  work  on  the 


ft 


65 


weighted  Sobolev  spaces  which  he  introduced.  (The  relevant 

norm  was  suggested  by  some  results  of  Nirenberg  and  Walker 
(1973).)  As  will  become  readily  apparent,  these  spaces  allow  both 
the  asymptotic  properties  and  the  junction  conditions  at  the 
star  boundary  to  be  readily  incorporated  into  our  argument. 

As  noted  in  §2.2  this  is  clearly  a  necessary  requirement  for  any 
uniqueness  argument. 

Let  p^O  and  let  I!  Ilp  denote  the  Lp  -norm  on  the  set 
C^°(©n,iRm)  of  C°  maps  with  compact  support.  Then  MPj(5= 

OO 

l)  is  defined  to  be  the  completion  of  C'0  with  respect  to 
the  norm 


2=  £ 

p's'4  M*. 


I<t(x) 


(A!- 


c' T)“f 


(s££J,<5e&) 


(3.1) 


where  ct2(x):  =  1  +  Ix!2,  =  ,  \fu.\:—T,/j.i  and  D;  f:  — 

al/x!f /(dx^K--dx^a)  .  If  f€C°°(iRn,iKm)  then  define  also 


a(f):={h:lRn 


>lRm|(h-f)€M 


(3.2) 


s.p 

The  MP0  spaces  are  not  just  the  usual  Sobolev  spaces  W 

OO 

which  are  defined  to  be  the  completion  of  C0  with  respect  to 
the  norm 


P,s=,rJlD"fl!p- 


lyups 


(3-3) 
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Cantor  originally  introduced  the  weighted  Sobolev  spaces 
in  order  to  apply  the  techniques  of  global  nonlinear  analysis, 
such  as  that  used  by  Fischer  and  Marsden,  to  physical  problems 
over  unbounded  regions  in  space  where  the  boundary  conditions 
include  the  asymptotic  behaviour  of  the  solutions  as  |x! — >oo, 
such  as  classical  fluid  mechanics  with  steady  flow  at  infinity,  or 

isolated  relativistic-  ftubd-S:  This-  rs  dome-  by  specs  tying-  an  My  f- 
space  of  solutions  which  satisfies  the  desired  asymptotic- 
conditions  and  studying  the  behaviour  of  the  particular 
differential  operators  given  by  the  physical  problem  on  these 
spaces,  which  is  what  we  want  to  do.  For  bounded  regions  and 
compact  manifolds  the  Sobolev  spaces  tV  '  are  commonly  used 
(c.f.  Ebin  and  Marsden  (1970)).  However  the  properties  of  a 
differential  operator  can  be  very  different  on  the  different 

spaces.  For  example  the  Laplace  operator  from  W  '  to  W  ~  is 

not  surjective  but,  as  we  will  see,  it  is  surjective  from  to 

Mg_2(5+7  for  certain  p,  s  and  6. 

The  asymptotic  behaviour  can  be  seen  by  letting  f  admit 

an  asymptotic  expansion  f(x)=|x!  E  (|x!  f  (x))  with  x  o.f  =0, 
fVo  .  Then  f€Mgi(5  iff  n>d  +  n/p  .  Any  f€M^  w^hich  does  net 

admit  an  asymptotic  expansion  can  clearly  stay  large  as 
| x I — >o°  on  an  angular  portion  A0(x)  wThich  goes  to  zero  as 
|x I — or  on  annular  rings  whose  ‘width  approaches  zero  as 
| x | — >oo.  If  f  is  C°  however  neither  of  these  situations  can 


. 


. 


arise. 


The  inclusion  maps 

Ulv6~ >ul2.6  for  si^s2  and  for  6^8Z  (3.4) 

are  easily  seen  to  be  continuous.  The  Sobolev  imbedding 
theorem,  which  is  one  of  the  main  hard  theorems  of  Sobolev 

space  theory  and  states  that  'Ws'p — >Ck  continuously  for 
s>k  +  n/p  (c.f.  Adams  (1975)),  shows  that 

Mp6 — >Ch  is  continuous  if  5^0  and  k+n/p<s  (3.5) 

since  Ilf lip_S|(5^llf lip>s  for  any  fEM^.  Recall  that  Ck  is  the 

00 

completion  of  C0  with  respect  to  the  norm 

Ilf max  sup  €S,n|DMf(x)i. 

k — 1  k  p 

Note  that  a  C  and  piecewise  C  "  function  may  still  be  in  Mg 
for  s^k  (if  its  asymptotic  properties  are  right)  but  not  in 

<+1.6- 

It  is  easy  to  see  that  partial  differentiation  induces  a 
continuous  map 


ak:M 


p 

s,di 


»M 


p 

s— 1,6+1 


which  is  also  linear  and  thus  C 


(3-8) 


If  p>l,  s>n/p,  O^k^s,  6,5'^- 0  then  any  pointwise 


! 


' 
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multiplication  tR  XtR  ^?R  induces  continuous  maps 


M^(»n,0®M^kia+k(finfam,)-^Map_ka+k(«n,iRm',)l  (3.6) 


s-k,<5  +  k' 


Mj|4(«nl»m)eMy»n,am,)-^M^+a,(iRn,«ir-’ 


s,  (5  +o 


(3.7) 


This  is  so m etimes  referred  to  as  the  Schauder  ring  property. 
Let  us  prove  (3.6).  It  suffices  to  show  that  (f,g) — >DM(fg)  is 

continuous  from  M£4©M£_w+t— >M£4+k+w  for  all  lyul^s-k  .  Using 
Leibnitz’  rule  one  gets 


d+k+l/u.!^^^ f  _q:  ^ ^  ^  ||/  p r  \ /  _<5+k+l//-t'L n~v 


D  (fg)l  p ^  £  ll(<j  D  f)( 


ll/fh/i 


O 


D  g)l 


for  some  positive  constant  C,  .  It  is  known  that  if  q+t>n/ 


-/P 


then  pointwise  multiplication  is  continuous  from  Wq,pelfl,p — >L 


(c.f.  Adams  (1975)).  It  is  easily  seen  that  D^fEM^, 6+]l/]  and  that 

<5  +  k+|ju— i/L  yu— i'  _  s-k-|u-i/|,p  .  I,,,,  , 

cr  D  gtl  so  since  s-|i/|+s-k-|/x-i/|^s>n/p  there  is 


a  constant  C2  such  that 


|  ^  d+k+lyulj-^  /u 


(fg)iip^Ca 


p,s,<5 


p,s-k,<5+k 


Before  proceeding  to  more  theorems  about  these  spaces, 

perhaps  we  should  pause  to  note  the  relation  between  the 
spaces  when  the  underlying  space  is  a  Riemannian  manifold 
N=(LRr  \g)  where  g  is  complete  instead  of  a  Euclidean  space 
En=(tRn,e)  where  e  is  the  Euclidean  metric  and  to  note  that 
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these  theorems  can  be  applied  to  appropriate  tensor  fields. 

Letting  I:  i!  be  the  norm  generated  by  g  on  a  tensor  space  T^N 
and  d\  be  the  volume  element  generated  by  g  we  can  define 

^s^TqN),  in  a  natural  way,  to  be  the  completion  of  c”(TqN)  with 
respect  to  the  norm 


for  p^l,  sEK\  6E&.  (Similar  definitions  can  be  made  for  Ws,y 

and  C  tensor  fields.)  Replacing  cr(x)  with  (l  +  d(x,0)  )  '  ~  where 
d  is  the  distance  function  generated  by  g  gives  an  equivalent 
norm  under  suitable  assumptions  on  g  .  In  fact,  from  Cantor 

(1979)  we  have  the  following. 

Lemma  3.4:  When  N  is  such  that 


lim  supla^(x)D/i(g-e)(x)|=0  for  l/d-m  where 

Ixl — »oo 

^  R*  k  ID 

m^rnax(k,s+2),  the  W  ,  C  and  MS(5  norms  generated 
by  g  and  those  generated  by  e  are  equivalent. 

9  f  9  9 

This  allows  one  to  treat  e  as  a  background  metric  for 
N.  Cantor  (1979)  uses  this  to  transfer  results  obtained  about 
the  dif f eomorphism  group  on  flat  space  (Cantor  (1975b))  to 
results  about  the  dif f eomorphism  group  on  N  .  We  will  simply 
state  the  relevant  definitions  and  theorems. 

Let 


l)^=^<EMsi(5(I)ly  1  exists  and  y 


(3.10) 


1 
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where  I  is  the  identity  map  on  &n  .  This  is  the  set  of 
diffeomorphisms  asymptotic  to  the  identity  at  infinity.  Such 
diffeomorphisms  will  not  destroy  the  asymptotic  behaviour  of 
tensor  fields  provided  they  fall  off  fast  enough  as  can  be  seen 
from  the  following  theorem. 

Theorem  3.5:  Let  p>l,  s>l-t-n/p,  <5^0.  Then 
is  a  topological  group  under  compositon  (right 
composition  is  smooth)  and  a  smooth  Banach  manifold 

(in  fact  it  is  an  open  subset  of  Mgtf(I)).  For  <5^6 

composition  from  5  ©2)^ — is  continuous  for  all 

k^s  and  if  f:tRn — is  such  that.  cr'M  Cl  is  a 

bounded  map  for  all  \/u]^s  then  composition  from 

(5  6~ — ^^k<5  continuous  for  all  k^s 

Proof:  This  corresponds  to  theorems  1.2,  1.6,  1.7  and  1.9  of 
Cantor  (1975b)  and  theorem  2.10  of  Cantor  (1979). 

For  convenience,  let  us  denote  both  Ms>{?  vector  fields  and 
1-forms  by  6  and  symmetric  covariant  2-tensors  by 

.  Since  specifying  a  spacetime  involves  specifying  a  positive 

3 

definite,  asymptotically  flat  metric  y  on  we  are  clearly 

going  to  be  interested  in  the  set 

ftl6=\g\g-e£Sl&  and  g  is  positive  definite!, 
the  set  of  Riemannian  metrics  which  are  asymptotic  to  the 
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identity.  has  a  natural  Banach  space  structure  for  s>n/p, 

being  an  open  cone  in  Sg5(e). 

Since  the  set  of  metrics  g  2  5?^  6  which  satisfy  the 
conditions  of  lemma  3.4  for  N  =  ( ttl  ,g)  is  dense  and  on  this  dense 
subset  2);>(S(N,N)  does  not  depend  on  g,  from  now  on  we  will 

consider  2^  to  denote  2)^5(En,En)  (l  e.  the  set  of 
diffeomorphisms  near  the  identity  with  respect  to  e  ). 

If  we  are  to  gain  an}7  insight  into  the  physicall_v  distinct 
solutions  which  correspond  to  the  solutions  of  the  linearized 

equation  t^'(a):TcrX — we  must  have  some  way  of  factoring 
out  the  infinite  dimensional  space  of  solutions  corresponding  to 
the  coordinate  freedom  on  X  .  Recalling  that  a  positive  definite 
metric  y  on  X  is  part  of  the  specification  of  a  spacetime  aEX 

dD  P 

it  is  clear  that  the  action  of  the  diffeomorphisms  on  j\s6  and 
their  resultant  orbits  are  going  to  determine  the  action  and 
corresponding  orbits  for  X.  With  this  in  mind  we  will  now 
present  the  relevant  theorems  dealing  with  the  action  of  the 

dif feomorphism  group  on  beginning  with  the  following 

theorem  of  Cantor  (1979) 

Theorem  3.6:  Let  p>l,  s>l+n/p  and  5^0.  Then  2^ 
has  a  continuous  action  on  2* s_1)f5+1  given  by 


A 


P 

5  —  1,6- 


p 

s— 1,6+1 


:(^.g  )’ 


. 


i.  ■  i.  i 


i  -  *  <■  •  *  ‘#1 
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00 

Moreover  A  is  a  C  function  of  g  and  if  gCR^_1+M+1 

then  p — >p  (g)  is  a  C  "  function  on  Xrs6. 

Proof:  Setting  g=e+h  with  hGSg_j  ^  we  have  p  (g )  =  p  (e) 

(h).  In  coordinates,  p\e)ii={dpl: / dx'Xa^1  / dx])bk:  and  /(h),j 

=(9^  /3x  )(3^  /dx^)(h°<^)kl  .  Using  the  Schauder  ring  property 

together  with  cp — >dp  being  smooth  from  into  M^_1i5t1  so 

that  dp  /dx  =dg+Bs  with  Bg  in  Mg_1<5+1,  it  is  seen  that  p — >p  (e) 

is  smooth.  For  h£Mg_1+M+1  theorem  3.5  implies  (h,$p) — >h*p  as  a 

function  from  Mg_lj(5+1X<^1  ^ — *Mg_1(5+1  is  C°'  in  h  and  C'L  in  p, 
completing  the  proof  e 

How  can  the  coordinate  freedom  be  factored  out  of 

g_l  (5  +  l?  Ideally  one  would  like  to  deal  with  the  quotient  space 

^s-u+i/^si •  the  ort>it  space  of  the  group  of  dif feornorphisms 
acting  on  the  Riemannian  metrics,  linearizing  the  equations  on 
the  quotient  space.  But  such  quotient  spaces  do  not  in  general 
have  a  manifold  structure.  In  fact  superspace,  the  quotient 

space  obtained  from  the  space  of  C  Riemannian  metrics  Ji 

and  C  diff eomorphisms  Jj  on  a  compact  oriented  manifold,  is 
not  a  manifold  (Fischer  (1970)).  This  results  because  for 
symmetric  metrics — i.e.  metrics  invariant  under  some 
continuous  coordinate  transformation — possessing  a  non— trivial 

isometry  group  I  =  \pC3)\p  (g)= g(  the  appropriate  dif feomorphism 

b 

group  to  form  the  quotient  space  with  is  $/Ig  since  an 


I 
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isometric  metric  (one  of  the  same  functional  form  as  g) 

generated  by  the  action  of  a  diffeomorphism  is  not  considered  to 
00  00 

be  new  .  Thus  in  the  quotient  space  where  the  crbits  are 
identified  to  points  the  neighborhood  of  a  symmetric  geometry 
will  not  be  homeomorphic  to  neighborhoods  of  generic- 
geometries  (geometries  with  no  symmetries).  However 

superspace  can  be  viewed  as  a  stratification  cf  manifolds  of 

geometries  with  given  symmetries,  the  more  symmetric 

geometries  forming  a  manifold  boundary  to  a  manifold  of 

geometries  of  less  symmetry. 

Because  we  are  dealing  with  dif f eomorphisms  asymptotic 

0  0 

to  the  identity  the  problem  of  symmetries  is  pushed  out  to 

0  Q 

infinity  and  in  fact  the  isometry  group 

l  (g)-g  for  g^«^s-ic5+i^  trivial  so  the  above  problem 

does  not  arise.  However  another  difficulty  resulting  from  the 
differentiability  conditions  and  therefore  appearing  also,  as  well 

S  D 

as  the  symmetry  problem,  in  the  case  of  W  metrics  and 

dif feomorphisms  (Ebin  (1970))  is  present.  This  is  related  to  the 

slicing  theorems.  Roughly  speaking,  when  there  are  no 

nontrivial  isometries  of  g  a  slice  at  g  is  a  subspace  which  is 
00  0  0 

orthogonal  to  the  orbit  through  g  and  which,  together  with  a 
small  neighborhood  of  the  orbit,  fills  out  an  open  neighborhood 

CO 

of  g.  When  everything  is  C  the  orbits  through  any  g  are 
C  submanifolds  and  a  slice  always  exists  (Eoin  (1968)).  Fischer 
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(1970)  uses  this  stacking  of  orbits  to  show  that  geometries 

with  the  same  symmetry  form  a  manifold.  In  the  Ws,p  or  Mp, 
spaces  of  metrics  the  orbits  are  submanifolds  but  are  in  general 

only  C  submanifolds  and  only  at  the  more  differentiable 
metrics  does  a  slice  exist.  There  thus  seems  little  hope  in 
establishing  a  useful  structure  for  the  quotient  space 

How  then  can  the  coordinate  freedom  be  dealt  with?  It 
turns  out  that  by  carefully  choosing  the  functions  with  which 
we  define  the  topology  of  X  we  can  make  the  unique  spherical 
solution  cr^X  have  enough  differentiability  so  that  a  slice  tf 
does  exist  at  cr  .  Thus  although  no  global  statements  can  be 
made  we  can  still  deal  locally  with  a  submanifold  of  physically 
distinct  spacetimes. 

Let  us  begin  the  presentation  of  the  appropriate  slicing 
theorems  by  stating  one  of  the  two  main  theorems  about  elliptic 

99  9? 

operators  which  are  dif ferentiably  close  to  elliptic  operators 
with  constant  coefficients  that  we  use  in  our  proofs.  We  will 
refer  to  this  theorem  as  Cantor’s  isomorphism  theorem.  Recall 

first  that  a  differential  operator  A=  X  a,.(x)DM  on  iRn  is  an 
elliptic  kth  order  differentia]  operator  if  there  are  no  real 

solutions  y<EtRn  of  X  a  (x)yM=0  for  all  xGitn 

1/7=1: 

Theorem  3.7:  Let  n>k  and  Aoc=  X  a-  Dp  a 

l/x|=k  h 

homogeneous  elliptic  operator  with  constant 


coefficients  on  tRn  and  A(x)  =  A00+  £  b  (x)D^  an  elliptic 

l/ubk  ^ 

operator.  Then,  if  p>n/(n-k),  0<i5<n-k-n/p,  s^k, 

s>n/p  and  b^GMg_k+|^j k_|/i|  with  continuous  outside 

a  compact  set  on  tRn  for  l/zl^n/p,  A  maps 

continuously  into  M^_k  6+k  with  closed  range  and  finite 
dimensional  kernel.  Moreover,  suppose  that  either 


(i)  £  lib  !!  _ kk_|„|<£  for  sufficiently  small  s,  or 

(ii)  There  is  a  continuous  curve  c  from  [0,1]  into 

the  space  of  bounded  linear  operators  between 

and  Mg_k(J+k  such  that  c(0)  =  Aoo,  c(l)=A  and  for  each 
t€£[0,l],  c(t)  is  an  injection  and  satisfies  the 
hypothesis  of  the  theorem. 


Then  A  is  an  isomorphism  of  onto  M[_k  6+k 

Proof:  Let  fCM^  .  By  Leibnitz’  rule,  for  |/?!^ s-k 


lla^+k+!^D^(bLtD/zf)llp^Ci,  £  _ ll(crk  '^Vb  J(c7 


k-l/til+lvl^r,  \/  (5+1  til+l|3— yL  fi  — v+ u 


D 


\v\*\p\ 


f)iL  and  it  i; 

P 


readily  seen  that  and 


P 

s+|«/|— |/S|—  |yu.|,(5 — |i/|+|jSl+l/u! 


Since  (s-k-|n|+lM!)+(s+l3y!“'/?l_lMD^s>n/ 


multiplication  from  M^_k+|/xlik_!/x|XM[’_|/i,  6+\g-^^s-k_d+k  is  continuou 
and  there  exists  a  constant  C2  such  that 


fi;p,s-k,6+k^C2'*b/j 


p,s-k.k-IA  !i  p,s,A 


(3.11) 


Thus,  as  shown  in  Cantor  (1979),  if  fR  denotes  functions  whose 
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support  is  in  JxIxPRi  with  llfEllpAj=l  then  r1  A  -  A^f  Rl  lps_k,5+k 

^1  Rlli?J|b,»'‘fR!ip,-1I.5+k=0  and  since  b^ec”  1  for  iplbn/p  all 

the  are  continuous  outside  a  compact  set  so 

lim  supla  (x)l=0.  This  satisfies  the  conditions  of  theorem  1.4 

Ixl— »°c  M 

of  Cantor  (1979),  where  b^GMg_k  k_|^|  with  s>k+n/p  was  assumed, 
the  latter  condition  ensuring  that  multiplication  was  continuous 
and  that  the  b„  were  continuous. 

h1 

The  outline  of  his  proof  is  as  follows.  That  A  has  finite 
dimensional  kernel  follows  directly  from  theorem  4.1  of 
Nirenberg  and  Walker  (1973).  The  closed  range  result  is 

obtained  by  showing  that  for  all  f€  \ker  A,  lifi!  a 

^cIlAfll  _ka+k  for  some  constant  C,  using  a  similar  argument  to 
the  one  we  use  in  the  proof  of  theorem  3.14.  The  isomorphism 
results  (i)  and  (ii)  follow  from  lemma  1.1  of  Cantor  (1979)  which 

shows  that  ATO:Mg5 — >Mg_kl5+k  is  an  isomorphism  together  with  (i) 
the  relation  (3.11)  or  (ii)  a  lemma  essentially  due  to  Sc-hauder 
which  states  that  for  two  continuous  linear  maps  A0,  A1 
between  Banach  spaces  E  and  F  such  that  A0  is  an 
isomorphism  and  such  that  there  is  a  continuous  curve  c  from 
[o,l]  into  the  space  of  bounded  operators  from  E  to  F  with 
C(0)=A0,  c(l)  =  A]  and  c(t)  an  injection  with  closed  range  for 

each  te[0,l],  then  A2  is  an  isomorphism.  E 

This  theorem  is  used  in  three  ways.  First,  we  use  the 
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closed  range  portion  of  the  theorem  directly,  together  with 
theorem  3.14  (essentially  due  to  Nirenberg  and  Walker)  in  §  4.4 
to  extend  the  results  about  the  kernel  of  an  elliptic 
inhomogeneous  operator  which  we  obtain  from  £'{c)  to  the 
kernel  of  the  corresponding  operator  obtained  from  ^'{o)  for  a 
in  a  neighborhood  of  <j.  (  £'{a)  is  not  elliptic.)  Secondly,  the 
isomorphism  property  of  some  specific  elliptic  operators  is  used 
to  obtain  the  vanishing  of  some  tensor  fields  which  satisfy 
certain  homogeneous  elliptic  equations.  The  following  corollarie 
are  useful  in  this  regard. 

Corollary  3.8.  Let  A  =div  od  be  the  Laplace-Belt rami 

operator  (where  div  £=— V1^  for  a  1-form  £  ).  If 
p>n/(n— 2),  s^2,  0^<5<n-2-ri/p  and  g  £5^  then 

A  :MP* — >MP  p  ,.0  is  a  continuous  isomorphism.  (The 

g  S,G  S  tZ, ,0  i  A 

subscript  g  will  be  dropped  when  the  metric  is  clear.) 

Proof:  (Cantor  (1979))  That  Ae:MP)(5  >Mp_2(5+2  is  an 
isomorphism  was  shown  by  Cantor  (1975a)  using  estimates  of 

Nirenberg  and  Walker  (1973).  For  ftM^,  Agf=(divg°d)i 

=-\g\-1/ZEd/dxX\g?y2g1)df/dx1)  where  |g|=det(gij)  .  It  is  easily 
i.j 

verified  from  this  together  with  that  Ag  satisfies  the 

n  2  ,  i2  t 

hypothesis  of  theorem  3.4.  Writing  AOQ—~,^d  /  dx  we  fina  that 
each  At=A00+t(Ag-A00)  is  an  elliptic  operator  with  no  lowest 
order  term  and  therefore  satisfies  the  maximum  principle  (c.f 


1 
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Friedman  (1969)  p.88).  In  particular,  each  At  has  no  non-trivial 
bounded  solution  and  so  it  is  an  injection  on  M  ^  ^  .  Thus 
condition  (ii)  of  theorem  3.7  holds  and  Ag  is  an  isomorphism. 

m 

Corollary  3.9:  Let  Kg:X^— >S^_1 5+1  be  given  by  £— -^g  . 

If  p>n/(n— 2),  s>2  +  n/p,  O^dCn— 2— n/p  and 

then  divg°Kg:X^ — -2,6+2  an  isomorphism. 

((divgc)i=-VJcij  for  a  symmetric  2-tensor  c.) 

Proof:  (Cantor  (1979),  Choquet  — Bruhat  et  al.  (1977a))  The 
proof  is  similar  in  spirit  to  the  above.  Since  (div  °K  £). 

=ai(,9ifj+5j(i)  =  XjeXs-2(6+2  has  a  solution  given  by 

£j=Ae  (xJ-Yai9jAe1X'i)eXs-2,6+2’  dlve°Ke  1S  an  isomorphism.  It  is 

easily  verified  that  div  +  K  satisfies  the  hypothesis  of  theorem 

6  © 

3.4.  Using  a  metric  gR=f Rg+(l-fR)e  where  fR  is  a  C  function 
which  equals  one  on  a  ball  of  radius  R  and  vanishes  outside  a 

ball  of  radius  R  +  l  Cantor  shows  that  A  =div  °K — *div  oK  in 
the  operator  norm  and,  using  some  technical  lemmas,  that  A 

©  P. 

is  an  injection.  (See  Choquet  — Bruhat  et  al.  (1979)  for  some  other 

methods  of  proof.)  E 

The  third  usage  of  theorem  3.7  is  to  obtain  splitting 
theorems,  or  more  specifically,  to  obtain  a  splitting  of  the 
tangent  space  of  X  at  a  .  This  splitting  turns  out  to  give  one 
subspace  as  the  tangent  to  the  orbit  through  <j  so  the  other 


. 
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orthogonal  subspace  can  be  taken  to  be  the  tangent  of  a  slice 

submanifold.  The  two  splitting  theorems  we  require  follow 

directly  from  corollaries  3.8  and  3.9  and  the  following  lemma. 

Lemma  3.10.  Let  E,  F  and  G  be  Banach  spaces  and 

f •  E  * F ,  h . F  * G  bounded  linear  maps.  Then  if  h°f  is 

an  isomorphism,  we  have  F=fE©ker  (h). 

Corollary  3.11:  If  p,  <5,  s  and  g  are  as  in  corollar}T 
3.8  then 


■^s-i|t5+i_d(Mgi(5)©/g 

where  /g=^x!?_1>(5+1!divg£  =  0  j. 

Corollary  3.12:  If  p,  d,  and  g  are  as  above  with 
s>2+n/p  then 


®s-i,a+i=Kg(Xs,a)®  Jg 
where  Jg=!c£Sj_w+1ldivgc=0j. 

The  first  corollary  was  first  given  by  Cantor  (197S).  The 
second  corollary  is  often  called  the  Berger-Ebin  decomposition 

as  it  is  one  of  the  decompositions  for  symmetric  covariant  W 
2-tensors  on  a  compact  manifold  that  they  obtained  using  some 

elliptic  operator  theory  (Berger  and  Ebin  (1969)).  (Note  that  divg 
is  the  formal  adjoint  of  Kg  ).  This  decomposition  was  used  by 

S,  p 

Ebin  (1970)  in  proving  the  slice  theorem  for  the  Riemannian 

metrics  on  a  compact  manifold  in  an  analogous  manner  in 


i 
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which  corollary  3.12  is  used  below. 

^  e  are  now  ready  to  investigate  the  orbits  of  2^  and  the 
slices  (when  they  exist)  of  ^l_16+v 

Theorem  3.13:  Let  p>n/(n-2),  s>2+n/p  and 

0^(5<n— 2— n/p  .  Let  2)^  act  on  2£g_1>(5+1  .  Then 

(i)  For  any  g£2^_u+1  the  only  isometry  of  g  in  3)vs6 
is  the  identity  I. 

(ii)  If  g^Rs_1+ki(5+1,  the  orbit  through  g, 

(0% =\<p  (g)l is  a  ck  submanifold. 

(iii)  If  g£R^4+1.  there  is  a  neighborhood  V  of  I  m 
&s,6  and  a  slice  of  the  action,  i.e  a  submanifold  tf  of 

+i  containing  g  such  that  (<p, g)1 — >A(^,g)  is  a 
homeomorphism  of  ¥>'•!/  onto  a  neighborhood  W  of  g 

in  fts-i,6+i  and  ®gn  ^  =k\- 

Proof:  (Cantor  (1979))  (i)  Since  ker  K  Cker  (div  °K  )=j0^,  the 

C>  too 

set  of  infinitesimal  isometries  is  trivial  so  the  isometries  of  g 
are  isolated.  Cantor  then  shows  that  the  isometry  group  of  g 
has  no  nontrivial  compact  subgroup.  The  idea  behind  this 
argument  is  that  if  there  were  such  a  subgroup  it  would  have 

finite  order  so  there  would  be  an  isometry  <p€L3)^  (  9 )  such 

that  9?k=I  )  for  some  k  .  The  asymptotic  properties  are  then 
used  to  show  that  there  is  an  entire  neighborhood  of  fixed 

points  of  9?  .  But  any  C1  isometry  which  fixes  a  neighborhood  is 


* 
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the  identity  so  the  subgroup  is  trivial.  Thus  all  of  the  orbits  of 

any  isometry  ^£<2^  must  be  unbounded  since  otherwise  p 
would  generate  a  compact  subgroup  (Kobayashi  and  Nomizu 
(1963)  p.47,48).  Along  an  orbit,  however,  the  displacement 
function  of  an  isometry  remains  constant  and  thus  for  cp  the 
displacement  must  vanish  and  <p=\. 

This  result  for  the  infinitesimal  isometries  generated  bv  a 
Killing  vector  field  was  obtained  independently  by  Choquet  and 
Choquet-Bruhat  (1978)  for  the  asymptotically  flat,  n  =  3  case, 
requiring  only  that  the  Killing  vector  field  goes  to  zero  at 
infinity.  We  will  have  occasion  to  use  this  result  as  well  as 
corollary  3.9  in  showing  that  certain  Killing  vector  fields  vanish 
as  this  result  does  not  require  as  much  differentiability. 

(ii)  Since  there  are  no  non-trivial  isometries  in  3)^,6  the 
map  Ag:  2)*g  — given  bY  Ag(p)=/(g)  13  injective  and  for 

g^s-i+M+i  il  is  °k  •  For  and  near  «*o.  <p°n  is  near 

I  and  <p\g)={(p°pQl)*0<fo(g)  so  it  is  sufficient  to  show  that  for  k^l 

A  '(I)  is  injective  and  its  image  splits  (theorem  3.3).  Let 
IGTjS^X^  and  let  <pt  be  the  flow  of  £  .  It  is  fairly 
straightforward  to  show  that  then  (Cantor  (1979) 

appendix).  Then  Ag'(I)(=(d/dt)(^t(g));=Kg^  =  i>^g.  From  the 

Berger  — Ebin  decomposition  (corollary  3.12)  and  ker  Kg=jO^  it 
follows  that  Ag  is  a  C  -imbedding  so  <Dg  is  a  C  submanifold 


. 


•• 
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(iii)  It  is  clear  that  there  is  an  open  neighborhood  Z  of  0 


in  Jg  such  that  g+c£^_1)(5+1  for  c£Z.  Then 


ygHg+clcezj  (3.12) 

is  an  imbedded  submanifold  of  Defining 

A:3)l'6Xtfg  +  — >cf  (g+c)  we  see  its  tangent  at  (I,g) 

is  given  b}" 

A'(I,g):  Xgj(5  ©Jg  >Sg_1(5+1:(((c)1 — >Kg£-l-c  which  in  view  of  theorem  3.6 
and  corollary  3.12  is  continuous,  onto  and  injective.  The  inverse 
function  theorem  and  its  corollaries  (c.f.  Lang  (1972))  give  the 
required  homeomorphism.  E 

In  order  to  investigate  curves  of  physically  distinct 
solutions  of  Einstein’s  equations  in  the  slice  we  need  to 
investigate  the  solutions  to  the  linearized  equations  off  of  the 
spherical  background.  This  is  clearly  impossible  to  do  directly 
by  solving  the  equations  on  all  possible  backgrounds. 

Fortunately,  however,  Nirenberg  and  Walker  (1973)  have  shown 

$9  cc 

that  for  elliptic  operators  which  are  dif f erent lably  close  to 
elliptic  operators  with  constant  coefficients  the  dimension  of  the 
kernel  is  locally  non-increasing.  In  other  words  a  sufficiently 
small  perturbation  of  such  an  operator  gives  an  operator  whose 
kernel  is  never  larger  than  the  kernel  of  the  unperturbed 
operator.  Thus  the  dimension  of  the  solution  space  ker  iZ  (a) 
will  not  increase  for  <j  close  enough  to  <j  This  is  used 
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together  with  Cantor’s  isomorphism  theorem  and  the  fact  that 
ker  £'{o)  is  O-dimensional  (l-dimensional)  in  the  static 
(stationary)  case,  which  we  prove  in  chapter  4  (5),  to  obtain  the 
desired  uniqueness  theorems. 

The  relevant  theorem,  stated  below,  is  the  second  of  the 
two  main  theorems  about  such  elliptic  operators  which  we 
require.  It  is  essentially  an  adaptation  of  theorem  4.2  of 
Nirenberg  and  Walker  (1973)  to  the  weighted  Sobolev  spaces  so  we 
will  refer  to  it  as  Nirenberg  and  Walker’s  theorem.  The  proof, 
given  in  detail  for  the  first  time  below,  follows  closely  the 
method  of  proof  given  by  Nirenberg  and  Walker  which  was  also 
followed  by  Cantor  in  obtaining  the  closed  range  result  of 
theorem  3.7. 

Theorem  3.14:  If  n,  k,  p,  s,  <5,  A and  A  are  as 
in  theorem  3.7  then  there  exists  an  e>0  such  that  if 

A  =  A„+  2  b  (x)D" 
l>u!^k  h 

is  another  elliptic  operator  for  which 

Hb/x_b  J‘P,s-k,k-l/A<£  (lM:^k) 

then  the  dimension  of  ker  A  is  less  than  or  equal  to 
the  dimension  of  ker  A. 

Proof:  Denote  the  dimension  of  ker  A  by  q  and  suppose 
the  theorem  is  false.  Then  for  all  integers  i>0  there  exists  an 

elliptic  operator  A^M^  ^Mg_M+k  given  by  Aco+^E^b/Li(x)D 


: 


1 


' 

- 


such  that  b^l!piS_kjk_|^|<i  1  for  |/r|^k  and  such  that  dim 

ker  Aj>q  .  A  positive  number  £0  can  be  found  such  that  any 

subspace  of  with  dimension  greater  than  q  contains  an 

element  of  norm  1  whose  distance  from  ker  A  is  at  least  £0. 

Choose  such  a  fi  in  each  ker  Aj,  llfjl  s ,5=1  .  Equation  (3.11)  of 
theorem  3.7  implies  !l Af i!!p s_  k,i+k=ll(A-Ai)tillp.3-W+k<C1i“llili!ip,a.6  for 
some  Cj  so  HAfJ!  _k>(5+k — >0  .  Therefore,  if  a  subsequence  of  {fj 
exists  which  is  Cauchy  in  M^,  since  A  is  closed  the 
subsequence  converges  to  an  element  in  ker  A  .  But  all  the  jfj 
lie  a  distance  a0  away  from  ker  A  .  Thus  to  prove  the  theorem 
we  need  only  find  a  Cauchy  subsequence  of  jfj. 

_  XI  CO 

Let  <^R:tR  — be  a  C  function  with  compact  support 
such  that  9?r(x)=1  if  Ixl^R,  9?R(x)=0  if  lxl^2R  and  |DM9?R(x)l^l  for 
all  fi.  Write  fi=^Rfi+(l“^R)fi  and  note  that  {fj  is  Cauchy  if 
\<pni  j  and  KI-^rW  are  Cauchy  for  some  R.  Let  BR=|x£  iRn 
I  Ixl^R}. 

For  all  R,  the  sequence  {<pRf  j  is  bounded  in  Ws,p(B2R)  . 

This  together  with  the  Rellich  compactness  theorem  stated  below 

implies  it  has  a  subsequence  (taken  to  be  all  of  |^RfJ  )  which 

converges  to  an  element  in  LP(B2R). 

Lemma  3.15:  (Rellich  compactness  lemma)  Let  0 ^ r < s  . 

Then  any  bounded  sequence  in  WS’P(B2r)  has  a 
convergent  subsequence  in  W  P(B2R). 


M 
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Proof:  (c.f.  Friedman  (1969)  p.31) 

Now  llA(^Rfi)llpiS_k^ll^RAfillp>s_k+llA((^lRfi)-^RAfilipiS_k  and 

^A^pf^-^pAfJ  has  support  in  B2R  and  is  bounded  in  Ws~k+1,p(B2R) 
since  the  highest  derivative  term  of  fj  cancels.  The  Rellich 
compactness  lemma  together  with  the  fact  that  ^RAfJ  is 
clearly  Cauchy  implies  ^A^f^  is  Cauchy  in  LP(B2R). 

As  long  as  lim  sup!ak  '^'b  (x)l  for  l/rl-k  is  small  enough  or 

|x| — »co  e 

b  (x)  is  uniformly  continuous  on  iRn  for  |gi|=k  there  is  a 

constant  C2  such  that  Ilf l!p  g^C2(llAfi!p+ Ilf llp)  (c.f.  Nirenberg  and 
Walker  (1973),  equation  (4.1)).  Since  the  former  condition  is 
always  satisfied  (c.f.  proof  of  theorem  3.7)  (and  the  latter  is  if 

s>n/p+l  since  then  is  C1  for  lgil=k  )  we  have 


A(^Rf  j)llp  +  H^Rfi  ^Rfjt:p) 

so  |^RfJ  is  Cauchy  in  Ws,p(B2R)  and  it  follows  that  )<pRfJ  is 
Cauch}7  in  Mpd  for  each  R. 

Since  A^M^ — >Mp_M+k  is  an  isomorphism  (lemma  1.1  of 

Cantor  (1979))  llfllp.Si^C3llA0C)fllpiB_ki<+k  for  all  fGM*a  .  Thus,  for  all 

R^l, 

ll(l-^R)f1llp,s.5£C3(!lA(l-^R)fillp.s-  +  A)(l  <^R)f i^'p,s-k,5+k  ■ ' 

As  seen  in  the  proof  of  theorem  3.7,  linr  li(Aoo- A)(l-^R)f  1!!ps_k  (5+k 

K  * 


*0  so  choosing  R  large  enough  we  have 


>l(l-^R)fillp.=.^C4llA(l-^R)fillps_k4+k 
£C4(ii(l“^R)Af1l.ps_k  <?  +  k+ll  A(l-f'R)fi-(l-^R)Afi!:piS_ki(j  +  k). 

The  sequence  |A(l— 9?R)f1-(l-f R) Af J  has  support,  in  B2R  and  is 
bounded  in  Ws  k"rl,p(B2R),  so  is  Cauchy  in  Ws  k,p(BgR)  .  By  passing 
to  a  subsequence  we  can  assume  it  is  Cauchy  in  MP_k  (5+k.  Sine 
|(l— 95R)AfJ  is  clearly  Cauchy  in  Mp_k(5+k 


|(l-^R)fJ  is  Cauchy  in 
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CHAPTER  IV 


ON  THE  UNIQUENESS  PROBLEM  FOR  STATIC  PERFECT  FLUIDS 

4.1  Introduction 

The  results  in  this  chapter  for  the  case  of  a  discontinuous 
matter  density  at  the  star  boundary  first  appeared  in  Kunzle 
and  Savage  (1980b).  The  extension  to  more  general  equations  of 
state  is  new.  As  indicated  in  §2.2  we  conjecture  that  the  set 

S  ^  ^  of  static  bounded  perfect  fluid  solutions  to  Einstein’s 
equations  with  a  fixed  equation  of  state  p(p)  satisfying 
conditions  (i)  to  (v)  of  §2.4  (and  fixed  surface  temperature  Tb  ) 
with  U  having  only  one  (nondegenerate)  critical  point  consists 
of  the  spherical  solutions.  Recall  that  we  are  going  to  consider 

the  set  S_  as  the  inverse  image  of  a  differentiable  map  YuX — >Y 

r 

between  some  suitable  weighted  Sobolev  spaces.  Since  all  static 
solutions  can  be  described  by  a  positive  definite  metric  y  and  a 
gravitational  potential  function  U  on  E,  which  by  the 
argument  about  spherical  solutions  can  be  taken  to  be 

dif feomorphic  to  tR  ,  we  want  to  consider  Sp  as  a  subset  of  an 
appropriate  Banach  space  formed  from  the  set  )(7,U)|. 

In  view  of  the  asymptotic  conditions  (1.35)  we  would  expect 

U  to  be  in  and  y  to  be  in  for  0^d<l-3/p  and  some 

s  and  p.  But  the  linearization  of  £  on  a  background 
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i.e.  £'{&),  does  not  approach  a  system  with  constant  coefficients 
as  o  tends  to  the  trivial  flat  space  solution.  In  other  words, 

for  a  fixed  p(p)  a  curve  m:[0,mcrit] — >u=(7,U)GSp  which  tends  to 

the  flat  solution  (7^=5^, U=0)  as  m — >0  is  continuous  but  is  not 
differentiable  at  m=0  .  Thus  Cantor’s  isomorphism  theorem 
cannot  be  applied,  as  one  would  expect  since  £'(o)  is  physically 
expected  to  have  a  1— dimensional  kernel  on  this  larger  space 

corresponding  to  a  curve  of  solutions  parameterized  by  the 
mass.  For  this  and  other  technical  reasons  which  will  become 
apparent  when  we  deal  with  the  appropriate  slicing  theorem  it 
is  more  convenient  to  fix  the  total  mass  and  consider  only  the 

subset  S  Imass(c0=m}  .  Then  for  any  two  cq,  cr2C;Spjn 

the  difference  U1— U2  has  a  fall  off  rate  faster  than  l/|x|. 

We  expect  Ub=U(p=0)  to  be  determined  by  m  and  p(p) 
but  as  noted  in  §2.3  we  cannot  yet  be  sure.  Therefore  we  will 

describe  the  set  Spm  by  the  triples  cr=(y, U,Ub)  .  There  is  clearly 
a  unique  spherically  symmetric  solution  in  this  set,  which  we 
will  denote  by  a  ■  By  choosing  the  appropriate  Banach  space 

structure  for  S  p  m  we  will  be  able  to  investigate  the  behaviour 
of  curves  of  physically  distinct  solutions  to  Einstein’s  equations 
with  constant  mass  m,  fixed  equation  of  state  p(p)  and  fixed 


surface  temperature  Tb. 


m  I  II  >•  -  '■  ■  I  ■ .  i  yfi 
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4.2  Choosing  a  Banach  manifold  structure  on  the  set  j(7,U,Ub)J 
We  are  going  to  want  to  apply  a  slicing  theorem  based  on 
theorem  3.13  to  the  space  which  we  pick.  But  for  0^6<1— 3/p 

theorem  3.13  is  applicable  only  for  g^Rg1+c5  (for  some  s  and  p  ) 

whereas  the  asymptotic  conditions  (1.33)  imply  This  is 

due  to  the  fact  that  on  decomposing  tensors  which  fall  off 

as  1/r  is  usually  impossible;  1/r  fall  off  is  generally  required. 
Fortunately  there  is  enough  coordinate  freedom  to  avoid  this 
problem. 

In  cartesian  coordinates  the  asymptotic  conditions  for  the 
general  stationary  metric  give 

7ij=<5ij+7ijlxl  14-0(|x|  2),  hi=h||x|  ~+0(|x|  3),  l^-mlx/Votlxl-2) 

with  xkdk7|j=0,  xkdkhq=0  .  But  it  is  well  known  that  for 

stationary  spacetimes  the  quantities  y\-  and  h)  can  be  made  to 
vanish  by  a  suitable  coordinate  transformation  that  is 
asymptotically  Euclidean  (i.e.  consisting  of  a  rotation  and  a 

translation).  For  example  y\-  and  h)  can  be  shown  to  vanish 

i  i  j  k 

by  demanding  that  the  coordinates  x  be  harmonic  so  7 
(Misner,  Thorne  and  Wheeler  (1972)  p.456)  or,  in  the  static  case, 
by  demanding  that  the  coordinates  are  asymptotically  harmonic, 

7ljrk=0(lxf3)  (Beig  (1979)).  (This  latter  result  can  probably  be 


* 


extended  to  the  stationary  case.)  Choquet-Bruhat  et  al.  (1979) 
have  dealt  with  this  problem  by  restricting  to  the  set  jg£ 

17  rij=0(  which  near  flat  spacetime  can  be  shown  to  be  a 

submanifold.  Whether  this  holds  true  far  from  flat  space  is  not 
obvious  and  we  will  not  take  this  route  but  as  we  will  see  at  th 

end  of  chapter  5  this  may  be  one  way  of  attempting  to  obtain  a 

surjective  map.  It  is  not  necessary  to  make  a  restriction  to 
harmonic  coordinates  since  the  desired  falloff  can  be  obtained 

explicitly  using  the  (2  +  1)  —  dimensional  formalism  where  U  is 
treated  as  an  intrinsically  defined  coordinate.  This  argument 
depends  essentially  on  the  surfaces  of  constant  U  being 
topologically  2  — spheres  and  involves  making  a  coordinate 

transformation  on  these  2— spheres.  From  now  on  we  will 

3 

therefore  consider  the  coordinate  system  of  S  =  3K  chosen  in 
such  a  way  that 

7ij=(5ij+0(|xr2),  h^Oflxr3),  U=-m/|x|+0(|x|  c)  (4.1) 

Further  coordinate  transformations  that  are  asymptotic  to  the 
identity  will  then  not  destroy  this  behaviour. 

Another  problem  arises  due  to  the  differentiability  at  the 
star  boundary.  Recall  that  for  an  equation  of  state  such  that 

lim pp_£=d>0  for  /(l+l)  that  M£Ck  for  k <1  but 

p— >o 

l  —  l 

limld*M/dUe|=oo  unless  s=l/(t+l)  in  which  case  M  is  C 

p— >0 

piecewise  Cl  (theorem  2.1).  Using  the  derivatives  calculated  in 
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§2.4  it  is  easy  to  verify  that 

^M^p,el(5+n^ci+c2 JocP(U+1)C  °P(p+p)  2dp  (4.2) 

for  some  positive  constants  and  C2  and  any  n>0  since  the 

only  unbounded  integral  is  the  i  normal  derivative  to  the 
star  boundary.  The  integral  on  the  left  side  of  (4.2)  is  bounded 
if  p<(l-£)/(£-(£+l)e)  .  In  order  to  apply  the  theorems  of  the 
previous  chapter  we  must  take  p>3  .  Thus  for  a  such  that 

(3£-l)/(3£+2)<a^/(£-fl)  M€Mj*i(5+n  (and  M^Mj+1^+n  )  for  3<p 
<(l-e)/(£-(£+l)a)  and  for  all  p>3  if  e=t/{t+ 1)  ).  If  (f-1  )/t<e 

^(31— 1)/ (31  +  2)  and  p>3  then  MGM^_ld+n  (and  M^M^+n  ).  We  can 
now  choose  the  equation  of  state,  and  thus  e,  to  be  fixed  once 

and  for  all  so  that  MEMg_3(5+n,  or  some  fixed  s^3, 

for  any  p>3  satisfying  the  above  conditions  for  an  upper  bound 
which  depends  on  the  value  of  e,  and  for  any  large  n. 

This  is  partly  why  we  have  chosen  p(p)  to  be  C  piecewise 

£+1 

C  in  the  interior  region.  Thus  even  though  p  is  uniquely 
determined  by  dp+(p+p)dU=0  when  Ub  is  given,  if  p(p)  is 

piecewise  C1  on  (0,pj  our  formalism  does  not  allow  the 
incorporation  of  different  differentiabilities  in  different  regions 
of  E.  The  other  reason  is  that  in  order  to  obtain  a  slice  at  the 

spherically  symmetric  solution  we  have  to  define  the  topology  in 
terms  of  a  smoothed  out  potential  which  is  a  little  more 
differentiable  at  the  spherical  solution.  Allowing  less 


1' 

■ 
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differentiability  of  p(p)  would  mean  that  a  different  smoothing 
potential  would  have  to  be  constructed  for  different  equations 
of  state.  This  is  therefore  not  a  fundamental  barrier  but  would 

be  messy  to  deal  with.  However,  as  we  noted  before,  it  is  not  a 
great  physical  restriction  to  assume  as  much  differentiability  as 
we  want  for  the  equation  of  state. 

Now,  using  the  asymptotic  conditions  (4.1)  and  taking  0^(5< 
1-3/p,  Einstein’s  equations  (1.29,  1.30)  and  the  weighted  Sobolev 

space  properties  (3. 6-3. 8)  indicate  that  U6Mg_1(5  but  that  we 

could  choose  <5+1  •  ^  we  t°°k  the  set  ^s-i  <5+i®-^s-i  <5  to  m°3el 

the  space  j(7>U)(  on,  the  action  of  the  dif feomorphism  group 

on  this  wmuld  then  only  be  C°,  as  we  will  see,  even  though 

for  7£Rg(5+1,  the  action  of  on  j^_1(5+1  is  C1  .  A  C°  action  is 

not  enough  to  obtain  a  slice  theorem  so  we  must  describe  the 

topolog}^  of  our  set  K7,U)(  in  terms  of  a  modified  potential  U 

which  we  make  to  be  in  for  the  spherical  solution. 

The  potential  U  is  not  in  because  the  normal 

derivative  of  order  s-2  of  M  at  the  star  boundary  is  not  in 

M06+n-  as  can  be  seen  from  equation  (1.30).  Define 

U=U-f  p2_E  (4.3) 

for  some  function  f  on  E  .  Let  [Vvu-  -VvuU]s=  lim_Vvu--  VvuU 

u  >ub 

-  lim  Vvu--VvuU.  As  seen  in  §2.4  this  is  proportional  to  the 
u->u; 


) 


1 
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difference  in  the  s  normal  derivative  of  U  with  respect  to  a 

CO 

C  coordinate  at  the  star  boundary.  Doing  similar  calculations 
to  those  of  §2.4,  using  dp/ dU=— 2M,  one  finds  that 


[Vvu-VvuU]s=lim(2W  2+fde  2U(1  °(2-£))ds  2M/dUs  2  (4.4) 

p — >o 


so  that  U2M^  provided  lim_=(2(2-e)d W2(Ub)d  2Ufc(1  e))  1  and  f>0 

u  >ub 

for  U^Ub  .  Since  it  is  not  easy  to  choose  such  an  f  that  is  a 
simple  functional  of  U  and  7  and  regular  in  the  whole  interior 

domain  (recall  W(Uc)=0  )  we  let  f=f0=(2(2-e)dW2(Ub)e~2U^1  E') _1= 

-q 

const,  so  that  for  the  spherically  symmetric  solution  U£Mg5, 

7£RL+s  •  This  is  also  convenient  for  the  stationary  case  so  in 
the  following  we  will  always  consider  f0  to  be  a  fixed  constant 
determined  by  the  unique  spherically  symmetric  m. 

For  the  remainder  of  this  chapter  we  take  0^<5<1—  3/p 
determined  by  the  asymptotic  conditions,  s^3>2+n/p  fixed  by 
the  equation  of  state  and  p>3  with  an  upper  bound  determined 
by  the  equation  of  state  as  seen  previously.  Define 


©p 

J  s-u+r 


)c=(7lj,U,Ub)l7ije^P_1>(5+1,U 


1J 


's— 1,5+1* 


U€MSP. 


1,5+1 


(4.5) 


where  U  and  U  are  defined  by  (4.3)  in  terms  of  U  and  U, 
respectively.  Since  U  and  U  determine  each  other  for  given 
Ub  we  will  indiscriminately  write  <r=(7,U,Ub)  and  c=(7,U1Ub)  . 
We  can  calculate  with  either  U  or  U  but  the  topology  is 


1 
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defined  in  terms  of  the  former.  Xl_16+i  is  clearly  a  Banach 
manifold  with  any  asymptotically  flat  solutions  of  Einstein’s 

equations  in  b— i.<5+i  having  the  same  mass  since  for  any  oq, 

a2^^s-u+i^^p  we  have  found  U1— U2  (  =UL — U2  in  vacuo)  to  fall 
off  at  a  rate  greater  than  l/r.  In  this  sense  the  fixed  mass 
just  becomes  a  boundary  condition  which  we  incorporate  into 

the  space  5)g_1(5+1,  but  note  that  X^_ls+l  is  much  larger  than  the 
space  of  asymptotically  flat  solutions. 

We  are  now  ready  to  investigate  the  action  of  the 

diff eomorphism  group  Xvs6  on  and  to  see  under  what 

conditions  there  is  a  slice  of  the  action.  The  following  theorem 

is  a  direct  analogue  for  the  space  ^s-ia+i  theorems  3.6  and 
3.13. 

Theorem  4.1:  (i)  <2)^  acts  continuously  on  ■2)g_1>(5+1  by 


A:2)^X^_u+1-^5>P_u+1:(^>(7,U,Ub))^(9P^7)>Uo^,Ub). 


s— 1,5+1' 


Moreover,  Auch — >A(^,u)  is  C°°  and,  if  <J^XVS_ 


s— l+M+i’ 


Aa:<p'—>A(cp,cr)  is  C  '  for  k  =  0  or  1. 

(ii)  If  £r£^>g_1+kt(S+1  (for  k=0  or  l)  then  the  orbit 
(0 a:=\A(cp,c7)\cpCXli6\CXl_li6+l  is  a  Ck  -submanifold. 

(iii)  If  aC5)g(5+1  there  is  a  neighborhood  V  of  I  in 
3)v z  and  a  slice  of  the  action,  i.e.  a  submanifold  X  of 

S  ,0 


s-1,5+1 


containing  a  such  that  (<p,c r')  >k(tp,cr')  is  a 


■ 


' 


n  .  . 
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homeomorphism  of  VX1/  onto  a  neighborhood  W  of  u 
in  ^s_u+1  and  ®aC\tf=\(T\. 

Proof:  (i)  Let  <j=(7,U,Ub)G5>g_1>(5+1.  Since 
theorem  3.6  gives  (i)  as  far  as  the  metric  is  concerned.  The 
function  U  is  explicitly  known  in  the  vacuum  region  (from 

(2.4))  and  it  is  clear  that  <r(x)<'liU'+<5  ^DMU  is  bounded  for  O^lgil^s  . 
Thus  theorem  3.5  shows  that  3)vs  tf+1(U) — >M^_1  d+1(U):(^,U)1 — » 

U°9?  is  continuous  (and  linear  whence  C  in  U  ).  But  if 
then  U£Mgjfi+1(U)CMgi5  .  Therefore  (9kU^Mg_1 6+1  and  (^,dkU)1 — » 

(aku)o^  is  continuous.  Hence  Aa  is  C1  in  cp  .  Since  the  action 
on  Ub  is  trivial  this  establishes  (i). 

(ii)  Since  there  are  no  nontrivial  isometries  in  the 

map  Aa  is  injective.  If  k=l  the  tangent  map  TA^  at  cp=l  is 

given  by  — >(£  ^y  ,£^1]  ,0)  where  £  is  an  Mgd  vector  field  on  E 
and  (^^7.^^U.0)^Sg_1(5+1©Mg_k(5+1©iR'A'TCT5)5_li(5+1  ■  But  Cantor’s  version 
of  the  Berger-Ebin  decomposition,  corollary  3.12,  for 
gives  the  splitting  Sg_li5+1=K7(X3i(5)©J7  where  J7=)cHSg_1(5+1ldiv7c=0(. 

Thus,  for  given  (c,u)GSg_lid+1©Mg_1(i+1  there  are  unique 

^HXgi(5  and  lCSg_15+1,  with  div7l=0  such  that  c=K7(£)+l  .  If  we 
let  v=u-^U  then  (c,u)=(/f7,^fU)+(l,v)  which  gives  the  direct 
sum  decomposition  of  Sg_li(5+1©Mg_1)(5+1  (since  (c,u)=0,  div7l  =  0  and 
\=£  7=>div7°K7(0=0  whence  £=0  whence  1=0  and  v=0  ).  Thus 


. 


96 


the  image  of  TjA^  splits  and  is  a  C1  —imbedding. 

(iii)  Choose  an  open  neighborhood  B  of  0  in  J  such  that 
7+1  is  positive  definite  for  1EB.  Now 

H(7+l.U+u,Ub+k)lieB,  uCM^+1,  k£iR i 

is  clearly  an  imbedded  submanifold  of  ^s-i  6+i  •  Defining 

><^)s-i,(5+i:(^’(7  +  l.U  +  u,Ub+k))1 — >(cp  (7  +  1),^  (U  +  u),Ub+k)  we 
have  for  its  tangent  at  (bn),  TMA:X^©J7©Mf_1>(5+1©l-^ 

Ss-i1<5+ieMs-iJ6+i0^:(^L-v+)1  >(K7^  +  L,,+^U+V,/c)  which  is  clearly 

continuous  and  onto  in  view  of  corollary  3.12.  It  is  injective,  for 

K7£+1=0  implies  (=0,  1=0  whence  <£^11=0  .  The  inverse  function 

theorem  and  its  corollaries  then  imply  (iii).  ■ 

In  particular,  the  tangent  space  to  ^s-16+1  cr+Pg^ 
splits,  i.e. 

T^'-W+i=T,©„®T„y.  (4.6) 

(This  theorem  has  a  trivial  extension  to  a for  all  k+2 

but  such  differentiable  spacetimes  cannot  be  in  Spn,  for  our 

fixed  equation  of  state  so  we  are  not  interested  in  them.) 

We  would  like  to  know  the  differentiable  structure  of  the 

set  Sp  £+1  of  static  solutions  of  Einstein’s  equations  with  a 

fixed  mass  m,  fixed  surface  temperature  Tb  and  a  fixed 
equation  of  state.  From  Einstein’s  equations  (l. 29,1.30)  it  is 


97 


clear  that  this  set  can  be  characterized  by  the  inverse  image  of 
0  under  the  map 

^'•^a-1,6+1  >^s-3,(5+30^s-3,5+3=:c^ 

(7ij.U,Ub)— »(Rir2aiUajU+2^7ij,AU-M),  (4.7) 

recalling  that  p  and  M  are  uniquely  determined  functions  of 
U  if  p(p)  and  Ub  are  given. 

Physical  intuition  leads  to  the  conjecture  that  £  1(0)=&o., 

i.e.  all  solutions  in  Sp  mn5)g_1 5+1  are  obtained  from  a  given 
spherically  symmetric  one  by  dif f eomorphisms.  According  to 
the  linearization  stability  technique  we  should  show  that  the 

tangent  map  ^'(u)=T(7^:T(75)^_1jC5+1— >T^((T  )=2  is  surjective  and  its 
kernel  splits  at  the  spherical  solution  cr  .  It  would  then  follow 
_ ^ 

th  at  X (0)  is  a  submanifold  with  its  tangent  space  at  cr 

tangent  to  the  orbit  .  This  result  would  not  be  as  strong  as 
the  above  conjecture  but  would  be  a  local  version  of  it.  As  we 
stated  before,  the  existence  of  a  Berger-Ebin  type  splitting  gives 
us  the  splitting  of  the  kernel,  namely,  equation  (4.6).  However, 
this  map  is  not  surjective,  as  we  will  see.  Its  kernel,  the 
tangent  space  to  the  orbit  through  a,  is  in  some  sense  too  large 
to  allow  it  to  map  onto  <2,  or,  in  some  sense,  2  is  too  large. 

But  it  is  by  no  means  obvious  what  sort  of  a  submanifold  of  2 

to  pick  to  try  and  make  £’{&)  surjective. 

However  all  is  not  lost.  The  set  of  equivalence  classes  of 
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^s-i.tS+i  a  neighborhood  of  <j  with  respect  to  the  action  of 

2)SP,  is  in  one  — to  — one  correspondence  with  the  slice  £  through 
cr  .  By  solving  the  linearized  equations  on  the  spherical 
background  and  finding  that  ker  0  and  then  using 

Cantor’s  isomorphism  theorem  3.7  and  Nirenberg  and  Walker’s 
theorem  3.14  to  conclude  that  the  kernel  is  zero  for  o  close  to 

cr  we  shall  show  that  if  there  was  a  C1  -curve  cr(A):(-£,£) — > 

Spmny  such  that  a(0)  —  <j  then  its  tangent  vector  would  have  to 
vanish  everywhere  so  that  in  fact,  cr(A)=cr  for  all  A  .  Less 

precisely,  there  are  no  non  — constant  differentiable  (Cl)  curves  of 
physically  distinct  solutions  of  Einstein’s  equations  with  fixed 
m,  Tb  and  p(p)  which  pass  through  the  spherical  solution. 

4.T3  "Linearization  of  £  to  find  ker  £'{<j) 

Let  (cir67ii.<5U,k)eT^^_w+1=S^w+1®MsP_u+1®S  .  For  <re5>P5+1 
and,  in  particular,  for  u-u  theorem  3.13  shows  that  one  has  the 
unique  decompositions 

‘WhTU  vi$lj=o,  <su=$+.qu  (4.8) 

where  4,  (J+1,  4,  ,qUEMp_1(i+1  and  k— of!",.)  (which  is, 

in  general,  not  equal  to  (dU)b  ).  Since  the  relation  between  U 
and  U  is  known  we  can  express  <5U  in  terms  of  <5U  by 


(5U=(l-4f0Mp)c5U+2(2-a)kf0p1  E(M-p)+(df0)p2  ' 


(4.9) 


) 


« 
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But  since  the  linearized  equations  are  considerably 

simpler  when  expressed  in  terms  of  <5U=:u=:£fU+^eM^_2 c 
and  k  we  will  perform  all  calculations  using  u,  c  and  k  .  We 
need  only  remember  that  the  differentiability  properties  at 

U=Ub  must  be  expressed  in  terms  of  U. 

The  variations  of  p  and  M  which  have  been  used  already 
in  obtaining  (4.9)  are  determined  in  terms  of  <5U  and  k  from 
the  integral  equation  (2.9)  obtained  from  the  equation  of 
hydrostatic  equilibrium, 

<5p=-2M6U+2k(M-p),  <5M=(dM/dU)<5U-k(2M+dM/dU).  (4.10) 

We  now  want  to  linearize  £  in  the  manner  of  the 
linearization  procedure  in  §3.1  where  we  saw  that  the  variations 
6(7  can  be  thought  of  as  the  tangent  to  a  curve  of  solutions.  In 

terms  of  £7^=0^  the  varied  Christoffel  symbols  are  then 

tfrt=!tr(VkCjr+V,Ckr“VrCjk) 

so  that  the  variation  of  the  Ricci  tensor  becomes 

«Rlj=t7rS(2VrV(loj)s-V1Vicrs-VrVsclj). 

Similarly,  although  we  do  not  use  the  variation  of  the  volume 
element  until  we  compute  the  variation  of  the  angular 
momentum  in  chapter  5,  note  that 

6det(7iJ)=det(7ij)c^. 
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A  lengthy  calculation  in  the  same  manner  then  gives  that 


y'(^)=(^'1(<T),^'2(CT)):T„5>sP_w+]=SsP_w+1®M5P_1,5+1eS 


P 

s— 3,<5+3 


©M 


P 

s_  3,6 +3 


is  given  by 


^'1(CT)(c.u,k)lj=|(VrVrclj+ViVjc") 
-V(ivrcj)r-3Rr(iCyRrscrs7ij 
+iRcij+cP(RiriR7ii)-2pcij 
+4Mu7ij+4U(i0jjU-4k(M-p)7ijl  (4.11a) 

^'2(a)(c,u,k)=Au-(dM/dU)u  +  ^Uiaic^ 

-U1Jcij+k((dM/dU)+2M)-UiVjc1J  (4.11b) 

where  Ui:=diU,  U^^VjdjU,  U1;=7irUr,  etc. 

Our  object  is  to  show  that  ^'(cr)(c,u,k)=0  for  (c,u,k) 
tangent  to  the  slice  if  implies  that  (c,u,k)  =  0  .  We  should  do 
this  first  for  a=cr  and  then  extend  to  a  neighborhood  of  <j  by 
the  use  of  Nirenberg  and  Walker’s  theorem  3.14,  But  this 
theorem  is  of  use  only  for  elliptic  operators  and  it  is  easily  seen 
that  the  operator  is  not  elliptic.  However  we  can  reduce 


d£'(<7)(c,u,k)=0 


(4.12) 


to  an  elliptic,  inhomogeneous  system  if  we  restrict 
consideration  of  this  system  to  the  slice  submanifold  if  .  To 


. 
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this  end,  define 

Li(a)(c,u):=^'1(c7)(c,u,k)+^K7odivyc+4k(M— p)7  (4.13a) 

and 

L2(a)(c,u):=^'2(o-)(c,u,k)+7(VU,div7c)-k((dM/dU)+2M).  (4.13b) 

It  is  straightforward  to  verify  that  L(cj)=(L1(o-),L2((j)): 

^s— 1,5+1® ^s— 1,5+1  >^s-3,(5+3®^s-3,(5+3  a  second  order  elliptic  operator 

and  that  the  equation  (4.11)  is  now  equivalent  to 

Li(a)(c,u)-yK7odiv7c=kq1(c7) 
L2(o-)(c,u)-7(VU,div7c)=kq2(o-) 

where 

q(a):=(q1(n),q2(a)):=(4(M-p)7,-(dM/dU)-2M).  (4.15) 

We  will  show,  for  nEl/,  that  the  inhomogeneous  equation 

L(cr)(c,u)  =  kq(a)  (4.16) 

for  (c,u)ESg_1  (j+1®Mg_1  S+1  implies  c,  u  and  k  vanish  provided  a 
is  close  to  <j  .  This  means,  in  view  of  (4.13)  and  (4.14),  that 
£'{o)  will  vanish  on  all  tangent  vectors  to  the  slice  1/,  which 

will  immediately  give  us  that  there  are  no  C1  curves  of 
solutions  in  £  \o)ny. 

First  we  will  show  this  on  the  spherical  background  a, 


(4.14a) 

(4.14a) 


h 
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then  we  will  apply  Nirenberg  and  Walker’s  theorem  3.14  to 
extend  it  to  neighboring  solutions. 

(We  will  continue  to  write  u  when  it  would  be  more 
correct  to  use  u,  keeping  in  mind  that  they  are  determined  by 
each  other  through  (4.9)  and  that  u  is  less  differentiable  across 
the  star  boundary.  Also,  in  the  case  of  equations  of  state  which 
give  s=3  we  will  regard  functions  such  as  dM/dU  as  piecewise 

C°  functions,  not  as  distributions.) 

The  splitting  of  TCT5Pg_u+1=TCT(PCT©TCTt/  allows  (4.16)  to  be 
investigated  on  the  two  subspaces  separately. 

Proposition  4.2:  The  equation  (4.16)  has  no  nonzero 
solution  (c,u,k)  on  T JO  . 

Proof:  Since  £{(p  a)=cp  <^(<t)=0  for  any  <p£L3)vs6  it  is  clear 
from  (4.13)  that  any  (c,u,k)(^T0.(^0.  is  a  solution  of  (4.16)  iff 

K  °div  °K  £=0  since  c=K^,f  for  a  unique  .  But  since  there 

7  7  7^  7^  * 

3 

are  no  Killing  vector  fields  on  iR  that  vanish  at  infinity 
(Choquet  and  Choquet-Bruhat  (1978))  it  follows  that  div7°Kyf=0 
which  implies  f=0  by  corollary  3.9.  Thus  (c,u)£Tcr(Dcr  is  zero, 
whence  L(<7)(c,u)=0,  so  k=0.  ■ 

Unfortunately,  on  T^T/  the  situation  is  considerably  more 
complicated.  We  are  unable  to  solve  equation  (4.16)  directly, 
even  with  the  extra  condition  div7c  =  0.  One  approach  used  to 

study  second  order  elliptic  differential  equations  on 
asymptotically  Euclidean  manifolds  is  to  make  a  conformal 
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transformation  to  a  compact  manifold  so  that  theorems  known 
in  the  compact  case  can  be  used  (Choquet-Bruhat  (1977)). 
However,  there  seems  to  be  no  such  compactif ication  of  E  which 

allows  the  coefficients  of  the  compactified  equations 
corresponding  to  (4.16)  to  be  non-singular.  Even  when  equation 
(4.16)  is  written  in  the  (2+1)— dimensional  formalism  and  a 
spherical  harmonic  expansion  is  done,  the  spherical  (f  =  0) 
portion  of  the  equations  results  in  several  coupled  second  order 
equations  which  can  not  be  solved  explicitly  even  in  vacuo  and 
which  seem  to  be  unamenable  to  a  maximum  principle  argument 
except  in  the  case  k=0,  where  it  can  then  be  shown  that  all  the 
spherical  variations  vanish. 

On  the  other  hand,  Kunzle  (1971)  has  linearized  the  static 
field  equations  (1.59-1.71)  in  the  (2+l)-dimensional  formalism 
with  the  requirement  that  the  central  value  of  the  potential  Uc 
and  the  pressure  pc  remain  constant.  This  means  in  our 
present  situation  that  k  =  0  since  it  follows  from  (2.9)  that 

k=<5Uc+(pc+pc)_1<5pc.  (4.17) 

These  linearized  equations  turned  out  to  decouple  into  a  three 
dimensional  scalar  second  order  equation  which  could  be  solved 
by  maximum  principle  arguments  and  into  equations  on  the 

compact  equipotential  surfaces  which  could  be  solved  using  a 
1  q m m a  which  is  valid  for  Rie m annian  2  —  m anifolds  with 

non-negative  curvature.  (An  error  in  the  proof  of  this  lemma 
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given  in  Kunzle  (1971)  is  corrected  in  appendix  1.)  The  maximum 
principle  arguments  which  Kunzle  gave  used  the  assumption 
that  p^3p  but  the  extension  to  the  case  p  is  straightforward 
and  is  presented  in  appendix  2  along  with  the  linearization 
arguments  of  Kunzle  which  we  use  a  little  later. 

We  will  therefore  reduce  the  present  problem  to  the  one 
treated  by  Kunzle  by  showing  that  (Sm=0  implies  that 
(5Uc=0  =  (5pc  or,  equivalently  in  view  of  (4.17),  that  k  =  0  and  Tc=0 

since  ^Uc='Tc+o£^U|u=u  =Tc  .  To  this  end  one  could  derive  (for 

(c,u,k)eTay  so  that  div^c=0  )  from  (4.16)  three  coupled  second 
order  scalar  equations  involving  u, 


■J>:=g  and  f:=$-W  Wu'ul 


(4.18) 


Since  only  spherical  variations  need  not  vanish  at  the  center 
and  the  general  variation  (on  the  spherical  background)  is 
obtained  by  superposition  of  spherical  harmonics  we  need  only 
look  at  the  1-0  component,  i.e.  we  can  confine  ourselves  to 
purely  spherical  variations.  However,  as  we  noted  above,  this 
system  is  still  too  complicated  and  we  must  proceed  differently 
to  make  some  progress. 

As  in  §2.3,  with  a  suitable  choice  of  coordinates  the 
Einstein  equations  for  a  spherically  symmetric  static  spacetime 
reduce  to  just  two  first  order  ordinary  differential  equations  for 
two  quantities  which  we  can,  for  example,  choose  to  be  (2.2)  and 
(2.3).  In  these  equations  r  denotes  the  curvature  radius  of  the 


■ 
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equipotential  surfaces.  We  now  consider  a  one-parameter  family 
of  spherical  perfect  fluid  spacetimes  described  in  terms  of  a 
fixed  radial  coordinate  r\  We  denote  by  r=F(X,r')  the  curvature 

radius,  where  X  is  the  parameter  of  the  curve  in  whose 

tangent  is  the  (spherical)  variation  (c,u,k)  at  one  fixed  solution, 

given  by  X=0,  and  we  choose  F(0,r')=r'  and  write  dF/9X|x=0=:f(r') 
The  coordinate  freedom  is  again  eliminated  as  previously:  We 
want  the  variations  to  take  the  form 

<5U=*+£fU,  <57ij=cij=«>ij+-£f7ij  (4-19) 

where  again  0,  but  now  <f=f(dr'/dr)d/dr'. 

Writing  (2.2)  and  (2.3)  in  terms  of  the  polar  coordinate  r' 
we  get 

dW2/dr'=(dF/dr')(2MWA“1-4W2F“1),  (4.20) 

dU/dr'=(dF/dr')WA-1,  (4.21) 

where  A  is  defined  through  (2.61).  Linearizing  these  and  using 
(4.12),  it  turns  out  that  f  drops  out  and  we  find  (after  some 

manipulations  using  V14ij=0  ),  for  X=0  where  r'  =  r, 

f2rZd't/dU-MrZi+{(l+pr8)^-l!r2(M-p)  (4.22) 

and 


W2r(dr/dU)d0/dU+2(l  +  2pr2)0-8Mr24  =  -8kr2(M-p)  (4.23) 


■ 
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_ 2 

where  0.=W  fCjU  UJ  .  These  equations  are  much  simpler  than 
those  that  result  from  (4.16)  for  spherical  variations  since  f 
drops  out  before  a  Lie  derivative  of  AU  has  to  be  subtracted. 
This  results  because  here  we  are  fixing  the  coordinate  system  in 
such  a  way  that  r  remains  the  curvature  radius.  In  (4.16)  the 
coordinate  r  is  allowed  to  vary  from  the  curvature  radius  and 
then  this  coordinate  change  is  factored  out. 

Lemma  4.3:  If  for  variations  (<5'7'ij,<5U,<5Ub)  on  the 
spherical  background  the  variation  of  the  mass 

vanishes,  6m=0,  then  also  <$UC,  6pc  and  <5Ub=k 
vanish. 

Proof:  Equations  (4.22)  and  (4.23)  can  be  explicitly 
integrated  in  the  vacuum  region  using  the  vacuum  solution 
(2. 4, 2. 5)  and  give 

d'^m  ^m  tanh  U  and  0=— 2m  1(5m  tanh2U  (4.24) 

where  one  integration  constant  was  put  equal  to  zero  in  view  of 
the  asymptotic  conditions  and  the  other  identified  with  the 

mass  change  (5m.  Requiring  T  to  be  in  Mg_2(5+1,  however, 

_o 

forces  dm  to  vanish,  since  tanh  U=m/r+0(r  )  for  r — >oo  . 

Since  07iTJ1U3=<f>i.U1UJ,  clearly  has  the  same  differentiability 

properties  as  4^  so  (Recall  that  we  are  on  the 

spherical  background  so  J •)  ^  easily  seen  that 


. 
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2  T) 

pr  eMs_2j(5+2  and  from  (2.3)  that  W(dr/dU)£M^_2i(5(l)  so  that  all 

the  terms  in  (4.23)  are  in  M^_2($+1  and  thus  in  Cs~3  except  the 

term  — 8Mr  (4/— k)  .  But  M^M^_2fi+2  because  of  its  behaviour  at 
the  star  boundary  while  all  derivatives  of  't  up  to  order  s— 3 
vanish  at  the  star  boundary  so  we  must  have  k  =  0  .  By  our 
assumption  (iv)  of  §2.4  on  the  differentiability  of  p  in  the 
interior  region,  M  is  more  differentiable  in  the  interior  so  we 

2  n 

see  that  Mr  'kEMg_2(5+2  and  equation  (4.22)  then  implies 

d1k/dUEMg_2 5+1  so  5+1CCs  2  Since  s^3,  'T  is  at  least  c\ 

In  the  interior  region,  then,  (4.22)  and  (4.23)  form  a 
regular  linear  homogeneous  system  with  continuous  coefficients 
in  (U  ,Ub]  with  zero  initial  conditions  at  Ub  so  that  </>  =  'T  =  0 
identically  (c.f.  Hartman  (1973))  and  therefore,  in  particular,  at 
the  center.  Equation  (4.17)  and  the  fact  that  ctU(xc)=0  now 
imply  that  6 Uc  and  <5pc  vanish.  ■ 

In  order  to  use  Kunzle’s  (1971)  result  the  relation  of  the 
linearizations  in  the  3  — dimensional  and  the  (2+l)-dimensional 

formalisms  must  be  determined.  Since  UEMg_1(5  has  only  one 
critical  point  that  is  nondegenerate  the  same  is  true  for  U 

sufficiently  close  to  UEM^^CC1  .  Thus  any  (thy)  sufficiently 

close  to  (U.y)  in  can  be  described  in  terms  of  (W.y)  in 

the  (2+l)-dimensional  formalism  of  §1.4  using  coordinates 

(x1):=(U,xA).  (Recall  that  y  is  the  induced  metric  on  the 
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U=const.  2— spheres.) 

Now,  W£M^+2  and  7 AB(U)£Ms_u+1(7AB),  7Ab  being  the 
metric  of  the  standard  unit  sphere.  Now  any  C1  curve  through 
a  with  its  tangent  in  can  also  be  represented  by  a 

curve  in  ^Li^+L;=rMs-2,«5+2©MJ-n(5+i(7AB)  through  (W,y)  whose 
tangent  at  (W,y)  is  in  S?_16+1  .  Let  (<57,c5U,0)eTCTF2)(5+1  and 

(w,c)€T(f-^ti(g_lifi+1  .  Noting  that 

u=x1,  71j=('91x1)(ajx1)w  2+(axK/ox1)(axL/axj)7KL 

and  linearizing  we  find  that  (5U  =  (5'x1  and  that 

(57ij=2w_2u(iaj)<5x1-2w"3uiujaw+2  7KLa(ixKaj)axL 

+aixKajxLd7KL 

where  6W  and  $yKL  contain  variations  due  to  the  change  in 

the  coordinate  system  (U,TA)  along  the  curve  in  ^bg_1(5+1  as  well 

as  the  tangent  to  the  curve  in  *^s-i  <5+1  (  ^7kl  ts  n°t  a  tensor), 
namely 

<5W=w+(auW)d 

^7kl“ckl+^  tl'L7KL+2  T  b^k7  L)A«5  x  . 

K  L 

(Recall  from  §1.4  that  Q=j  0KL  where  QKL  is  the  second 
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fundamental  form  of  the  U=const.  topological  spheres.  Also, 

_ A 

rBC  are  the  Christoffel  symbols  of  the  metric  7  .) 

Thus,  defining  <f:=((5U ,6xk)  in  the  (x1)  -coordinate  system 
so  that  (5U=('1=^'1Ui  and  we  have 

57ij=2V<ifj,+8ixKajxLcKL-SUiUj¥“3w. 

Since  the  central  potential  and  pressure  are  fixed  to  first 
order  and  a  variation  of  higher  order  will  not  affect  the  first 

order  quantities  w  and  c'KL  the  result  of  Kunzle  (1971),  given  in 

appendix  2,  yields  w=0  and  "c  (k^l)  for  some  vector  field 

yA  on  the  U  =  const.  hypersurfaces. 

One  can  choose  a  coordinate  system  (U,x  )  such  that 

71A=0  identically.  Then  (^x7)ia=^(iXa)=0  which  implies  auxA=0  . 
But  yA  must  vanish  at  infinity  since  we  require  so 

y A= 0  .  Thus  <5U  and  67^  are  merely  Lie  derivatives,  hence  so  is 
<5U,  and  we  conclude  that  if  (<57,<5U,k)  satisfy  (4.16)  and  are  in 

T^  they  must  all  vanish.  Summarizing  this  we  have 
Theorem  4.4:  The  operator  equation 

L(cr)(c,u)=kq(a) 

where  L  and  q  are  defined  by  (4.13)  and  (4.15), 
respectively,  and 
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(c,u,k)€!T0.5)s_lj5+1— Sg_1>(5+1ffiM 


s— 1,<5+1 


©& 


implies  (c,u,k)=0. 


4.4  Curves  of  solutions  in  y 

We  now  wish  to  extend  this  result  to  slices  supposedly 
representing  nontrivial  static  deformations  of  the  spherically 
symmetric  solution  in  order  to  show  that  there  are  really  no 
(sufficiently  regular)  such  deformations.  It  seems  natural  to 
apply  Nirenberg  and  Walker’s  theorem  3.14  to  get  this  result,  but 
there  are  still  a  few  problems.  In  order  to  apply  the 
Berger— Ebin  decomposition  to  the  variations  of  the  3-metric  and 

thus  for  the  proof  of  theorem  3.13  we  need  a  fall  off  rate  of  l/r 
or  higher.  However,  the  isomorphism  theorem  3.7  and 
Nirenberg  and  Walker’s  theorem  both  require  a  l/r  fall  off,  the 
former  being  required  in  order  to  have  a  closed  range.  Without 
the  closed  range  result  the  presence  of  the  inhomogeneous  term 
would  make  it  hard  to  extend  theorem  4.4  to  a  neighborhood  of 
cr.  We  must  therefore  now  investigate  the  operator  L(cr)  also 

on  the  larger  space  Sg_1;(5©Mg_15. 


Writing  L((j)  =  L00+^  E  b^D  it  is  readily  verified  from  (4.11) 
and  (4.13)  that  each  contains  derivatives  of  U  and  7  up  to 

order  2-|/rl  for  all  I jjl\^2  .  The  properties  (3. 6-3. 8)  of  the 

weighted  Sobolev  spaces  then  easily  show  that  for  oyyg_1(5+1, 
b^Ms-3+l//l<5+3-l/ul  for  all  If  in  addition  the  b^  are  smooth 
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enough,  (they  actually  only  need  be  continuous  for  |/r|=0  )  in  the 
asymptotic  region  the  conditions  on  b  for  theorems  3.7  and 

3.14  will  be  satisfied.  This  is  certainly  the  case  for  any 
since  we  have  only  compact  matter  distributions  and  the 
asymptotically  flat  solutions  to  the  field  equations  are  analytic 
in  the  vacuum  region  (Muller  zum  Hagen  et  al.  (1970)). 
Application  of  theorem  3.7  gives  the  following. 

Proposition  4.5:  For  uHjPg_1|(5+inSp  the  maP 

L((t).P.— Ss_1(5®Ms_1>(5  >Ss_3($+2©Ms_3(5+2=:Q  (4.25) 

is  a  continuous  linear  operator  with  finite 
dimensional  kernel  and  closed  range. 

If  we  try  to  solve 

L(cr)x=kq(<j)  for  kHtR,  xEP  (4.26) 

using  the  method  used  in  the  proof  of  lemma  4.2  we  find  that 
'F  and  </>  for  the  spherical  variation  are  uniquely  determined 
by  (5m  in  the  vacuum  region.  They  determine  therefore  k  on 
the  boundary  and  hence  the  value  6UC  .  Thus  there  is  a 
one-dimensional  solution  subspace  of  P©iK,  spanned  by  (x0,k0), 
say,  where  (x0,k0)  is  the  solution  corresponding  to  <5m=m  in 
the  vacuum  region. 

Proposition  4.6:  For  cr  in  a  neighborhood  of  a  in 
•^s-i<5+i  the  °Perat°r  L(ct):P — >Q  is  injective. 
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Proof:  Let  first  o-u  .  Putting  k=0  in  (4.23)  means  (5Ub=0, 
i.e.  (/oc+pc)<5Uc+6pc=0  and,  as  we  just  saw  above,  6m=0. 

Therefore  there  is  no  spherical  solution  for  which  the  central 
pressure  and  potential  stay  fixed.  Hence,  just  as  in  theorem  4.4 
there  are  no  nonspherical  variations  either.  The  operators  L(cr) 
and  L(<t)  on  P,  however,  satisfy  the  hypothesis  of  theorem  3.14 

since  llo--allp>s_lf5+1  small  implies  Hb/u(^)-b/U(a)lip>s_3+1/uLk_,Ml  is  small, 
so  that  the  kernel  of  L(<t)  can  have  no  higher  dimension  than 
ker  L((7)=|0|.  e 

But,  finally,  we  must  again  consider  equation  (4.16)  for  x= 

(c,u)£Pe=Sg_15+1©Mg_16+1  and  k  a  constant  to  be  determined.  By 

(3.4)  Pe  is  a  subspace  of  P  but  not  a  closed  one.  On  the  other 
hand,  for  <j,  we  know  that 

x£Pe,  L(cr)x=kq(a)=>x=0,  k  =  0. 

This  means  that  q(o-)^L(cr)Pe  while  according  to  the  above 
there  is  a  unique  xQOP  such  that 

L(cr)x0=q(cr).  (4.27) 

Now  suppose  that  a  is  in  a  neighborhood  of  a  such  that 

llL(a)  — L(ct)H<£1  (operator  norm) 

and  Hq(u)-q(cr)llpiS_3j(5+3<£2  (4.28) 


for  some  small  £q,  £g>0  .  Since  kerpL(cr)  0  there  is  at  most 
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one  solution  xEP  of 


L(o-)x  =  q(a).  (4.29) 

We  will  show  that  x  cannot  lie  in  P 

Since  L(or):P  >Q  is  injective  and  has  closed  range  it  has  a 
bounded  inverse  whence  there  exists  C>0  such  that 

llL(CT)llp,s_3.«+22cllxllps_w  for  all  x€P.  (4.30) 

From  (4.27)  and  (4.29)  we  have  q(cr)-q(cj)=L(o’)(x-x0) 

4-(L(cr)— L(cr))x0  which  leads,  together  with  (4.28)  and  (4.30),  to 

h-X0llp,s_w^(C-£1)_1|lL(ff)(x-X0)llps_3,!+2 

£:(C-£ir1(llq(<r)-q((7)llpiS_34+3+llL(<7)-L(cT)ll  llx0llpia_li4) 

<(C— sq)  ^IIxqII  _1i5£1+£2)  i.e.  the  solution  x  is  arbitrarily  close  to 

x0  in  P  for  small  enough  £q  and  On  the  other  hand  it  is 

not  hard  to  see  that 


d(x0,Pe):=inpf!lx0— yllps_w>0. 

(Observe  that  x0  is  explicitly  known  in  the  vacuum  region  and 
has  an  asymptotic  expansion  starting  with  a  term  proportional 

to  l/r,  while  for  yEPe  the  mean  falloff  rate  is  higher). 
Together  with  the  remarks  made  after  equation  (4.16)  this  leads 
to  the  following  uniqueness  result  for  static  spacetimes. 

Theorem  4.7:  Let  c:[0,l] — be  a  C1  curve  such 


■ 
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that  c(0)  =  <7  is  a  spherical  solution.  Then  c  is 
constant  if  the  slice  y  is  contained  in  a  small 
enough  neighborhood  of  a. 

Proof:  Since  is  a  C1  Banach  manifold  and  c(t)£Sp 


implies  that  ^(c(t))=0  whence  y'c^(dc/dt)=0  for  all  t.  But 
since  c  is  also  tangent  to  the  slice  d/7  this  is  equivalent  to 
L(c(t))(dc/dt)=kq(c(t))=0  which  implies  dc/dt  =  0  for  all  t. 


CHAPTER  V 


ON  THE  UNIQUENESS  PROBLEM  FOR  SLOWLY  ROTATING  STATIONARY 

FLUIDS 

U.l  Introduction 

In  this  chapter  the  advantage  of  the  Banach  space 
formalism  developed  for  the  static  case  will  become  apparent  as 
it  can  be  applied  also  to  the  stationary  case  resulting  in  a 
similar  type  of  uniqueness  result.  Namely,  there  are  no 

nonconstant  differentiable  (C1)  curves  of  rigidly  rotating,  perfect 
fluid  solutions  of  Einstein’s  equations  with  fixed  mass  m,  fixed 

surface  temperature  Tb,  fixed  equation  of  state  p(p)  and  fixed 
angular  momentum  J  which  are  close  to  the  spherically 
symmetric  solution,  i.e.  which  are  slowly  rotating  so  that  J  is 
small.  This  result  was  presented  for  equations  of  state  wrhich 
result  in  a  discontinuous  matter  density  at  the  boundary  of  the 
star  by  Klinzle  and  Savage  (1980c)  using  results  obtained  from 
the  linearized  (2+l)-dimensional  formalism  of  Einstein’s 
equations  for  stationary  spacetimes  which  was  given  in  Kunzle 
and  Savage  (1980a).  This  latter  paper  contains  an  analogous 
investigation  for  the  stationary  case  as  was  carried  out  by 
Kilnzle  (1971)  for  the  static  case.  In  other  words,  it  is  a 
linearization  keeping  the  central  potential  Uc  and  pressure  pc 
fixed  instead  of  the  more  physical  constraint  of  fixing  the  mass. 
We  will  be  able  to  use  the  relation  between  these  two 
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linearizations  which  was  obtained  in  the  last  chapter  to  obtain 
our  result. 

The  extension  to  more  general  equations  of  state  which 

satisfy  conditions  (i)  to  (v)  of  §2.4  follows  easily  as  in  the  static 
case  but  one  additional  restriction  must  be  made  for 

mathematical  reasons.  In  analysing  the  linearized  equations  we 
need  to  make  a  power  series  expansion  in  the  interior  region  so 
we  must  also  assume  that  the  equation  of  state  is  analytic.  We 

will  assume  then  that  S denotes  the  set  of  asymptotically 
Euclidean  stationary  axisymmetric  solutions  of  Einstein’s 
equations  with  a  fixed  surface  temperature  Tb  and  a  fixed 
analytic  equation  of  state  p{ p)  satisfying  conditions  (i)  to  (v)  of 
§2.4.  To  demonstrate  the  similarity  to  the  static  case  and  to 
avoid  a  proliferation  of  symbols  we  will  use  a  tilde  above  a 
symbol  used  to  denote  a  map  or  space  in  the  static  case  to 
denote  the  corresponding  map  or  symbol  in  the  stationary  case. 

Recall  that  Lindblom  (1976)  has  shown  that  equilibrium 
fluid  configurations  satisfying  our  conditions  for  a  stellar  model 
are  necessarily  axisymmetric.  However,  his  method,  which  is 
similar  to  Hawking’s  (1972)  proof  that  stationary  black  holes  are 
axisymmetric,  uses  analytic  extension  methods  which  are 
somewhat  difficult  to  compare  rigorously  with  our  asymptotic 
conditions.  Nevertheless  this  result  shows  that  we  can  restrict 

to  the  axisymmetric  case  without  a  serious  loss  of  generality  so, 
as  in  §1.3,  we  take  a  fixed  vector  field  p  on  E  which  behaves 

asymptotically,  in  cartesian  coordinates,  like  rj=x  d2~x  d1+0(|x|  )  . 


.  3 
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Thus  we  are  eliminating  some  coordinate  freedom.  The  angular 
m o m entum  J  is  then  given  by  (1.38)  or  (1.39)  where  cc  is  a 
1-form  such  that  H  =  da  (c.f.  §1.3). 

The  slicing  theorem  we  require  in  the  stationary  case 
again  depends  on  the  existence  of  a  Berger  — Ebin  type 
decomposition  but  it  also  depends  on  the  decomposition  of 
vector  fields  given  by  corollary  3.11.  However,  since  we  have  the 

asymptotic  properties  (4.1)  which  give  a  fall  off  rate  of  1/r2  for 
both  7  and  a  this  creates  no  problems.  In  fact,  using  the 
vacuum  field  equations  together  with  the  asymptotic  properties 
of  7]  and  (4.1)  it  is  seen  that 

hi-2j|xr3(4-3ixr2x3xi)+o(ixr4)  (5.1) 

while  the  vector  potential  a  is  given  (up  to  a  gradient)  by 

ai=2J|x|  3aik3xk+0(lxl  3).  (5.2) 

In  order  to  obtain  an  elliptic  operator  to  which  we  can 
appl}^  Cantor’s  isomorphism  theorem  and  Nirenberg  and  Walker’s 
theorem,  Einstein’s  equations  (1.22—1.24)  and  the  thermodynamic 
equilibrium  equations  (1.25—1.28,  1.11)  must  be  modified  so  that 
the  linearized  system  will  consist  only  of  as  many  second  order 
equations  as  there  are  unknown  functions. 

We  use  the  vector  potential  a  to  do  this.  Equation  (5.2) 

implies  «GXg_1(5+1  for  0^W<l-3/p  for  some  p  and  s  which  we 


■ 
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will  see  later  must  satisfy  s^2  and  p>3  so  that  corollary  3.11 
then  implies  a  will  be  unique  if  we  require 

div7«=-V1ai=0  (5.3) 

Using  a  instead  of  h,  where  h1=£lr'“i9ras,  makes  equation  (1.24) 
into  a  second  order  equation  and  eliminates  equation  (1.28)  as  it 
is  then  identically  satisfied. 

In  the  matter  region  D  we  have  from  (1.27)  that  <£ecx=0 
for  the  unique  a  defined  by  H  =  da  and  (5.3)  so  we  can  integrate 
the  rigid  rotation  condition  equation  (1.26)  and  find  that 

0Ja— uT  !e  2U=:-a=const.  (5.4) 

From  the  conditions  that  6  is  a  symmetry  of  the  spacetime, 

(1.25),  it  follows  that  6  vanishes  at  the  center,  whence  uc=eUc 
and  therefore 


rp— 1 

a  =  Tc  e 


-u„ 


(5.5) 


From  (5.4)  we  then  find 

u=Te2U(a  +  @J  a)  (5.6) 

and,  solving  for  T, 


T  2=e2U(a  +  0JcO 


2 


The  vector  field  6  is  not  defined  in  the  vacuum  region 


- 
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S\D,  but  it  is  an  infinitesimal  isometry  in  D  .  Since  it  is  also 

tangent  to  the  topological  2— spheres  U=const.  by  (1.25)  and 
vanishes  at  the  center  it  m ust  have  closed  orbits  and  thus 

represents  an  infinitesimal  rotation.  It  follows  in  the 
axisymmetric  case  that  unless  the  space  time  is  spherically 
symmetric  8  must  be  proportional  to  the  given  rotational 
generator  r),  i.e. 


6=br],  b=const.  in  D,  (5.8) 

since  there  are  no  closed  2-dimensional  subgroups  of  the 
isometry  group  of  a  Riemannian  2-sphere  (c.f.  Kobayashi  (1972) 
P-47).  Note  that  b/a  is  essentially  the  angular  velocity.  Recall 
that  7)  is  globally  defined  on  E  and  satisfies  (1.36)  and  (1.37). 

Since  we  assume  axial  symmetry  it  might  seem  reasonable 
to  eliminate  an  angular  variable  and  work  in  the  (r,6>)  -half 
plane.  However,  for  global  arguments  this  quotient  space  is  not 
particularly  convenient  and  we  could  not  immediately  apply  the 
Banach  space  techniques  which  we  used  in  the  static  case. 
Similarly,  if  we  were  to  consider  r\  as  a  fixed  vector  field  on  E 
when  parameterizing  the  set  of  equilibrium  configurations  the 
structure  of  the  group  of  dif f eomorphisms  leaving  nr\  invariant 
would  be  more  complicated.  So  it  seems  best  to  treat  rj  also 
as  a  variable.  In  view  of  the  fact  that  r]  is  a  Killing  vector 
field  it  satisfies  the  second  order  equation 


divrU^7=-V'(Virjj+VjT)i)dx:l=0. 


(5.9) 
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For  the  temperature  T  we  have  from  the  equation  of 
hydrostatic  equilibrium  (l.ll)  that  in  D 

T=Tbexp  J^(p(p)+P)^1dp  (5.10) 

where  Tb  is  the  temperature  on  the  surface  3D  of  the  star, 
defined  by  p=0  .  Since  for  each  different  value  of  Tb  of  the 
surface  temperature  we  get  an  otherwise  identical  model  with 
the  temperature  differing  at  each  point  of  D  by  a  constant 

factor,  just  as  in  the  static  case,  we  keep  Tb  fixed  once  and  for 
all. 

If  the  tensor  fields  7,  U,  a,  77  and  the  constants  a  and 
b  are  given  then  our  model  will  be  uniquely  determined  since 
8  is  given  by  (5.8)  in  the  domain  D,  T  by  (5.7)  and  finally  p 
and  p  as  functions  on  D  by  (5.10)  (or  equivalently  by  (l.ll)) 
since  we  also  consider  p  as  a  function  of  p  to  be  given.  For  a 

3 

fixed  manifold  £  which  we  take  to  be  dif f eomorphic  to  iR  and 
which  thus  can  be  considered  to  be  provided  with  a  fixed 

coordinate  system,  we  can  therefore  describe  the  set  S  ^  j  by  a 

subset  of  a  set  ZP  of  sextuples  a=(7,U,a,?ya,b)  which  satisfy  the 
appropriate  differentiability  conditions  (determined  by  p(p)  ) 
and  the  asymptotic  conditions  (4.1,  1.37). 

Then,  as  in  the  static  case,  we  can  consider  Sp(p)  to  be  the 

rxj  r>j  r^j 

inverse  image  of  zero  of  a  map  — >2  between  Banach  spaces. 


i 
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r^j 

We  will  take  the  map  y  to  be  given  by 

2:(7,U,a,77,a,b)— >(Rir2UiUj+^e4Uhihj(p+p)e"4UT20i0J 
+(2pe  +{p+p)T262e  4J)y1y  AU--J(p+3p)e-2U+^-e4Uh2 

-(p+p)e~4UhV,  e_4UVr(V[rak]e4U)+(p+p)e_6%T@k,  VrV(rpk)) 

(5.11) 

where  6,  T,  v,  and  p  are  defined  by  (5.8),  (5.7),  (5.6)  and 
(5.10),  respectively. 

We  have  chosen  this  particular  map,  so  as  to  have  as 
many  second  order  equations  as  unknown  functions.  It  is  clear 
that  the  first  three  components  of  the  equations 


2(ct)= o 

correspond  to  Einstein’s  equations  (1.22  —  1.24)  respectively,  and 
the  last  to  (5.9).  We  have  not  included  the  stronger  equations 
(1.36)  nor  the  gauge  condition  (5.3)  on  a  since  it  will  turn  out 
that  the  chosen  equations  already  determine  the  equilibrium 
configuration  uniquely  up  to  the  two  integration  constants  m 
and  J  and  arbitrary  dif f eomorphisms  (asymptotic  to  the 
identity)  and  arbitrary  gradient  fields  added  to  a  .  This  occurs 
because  the  gauge  freedom  of  a  is  factored  out  by  the  slice 
theorem  we  give  in  the  next  section  in  the  same  way  that  the 
coordinate  freedom  is  factored  out,  i.e.  by  restricting  to  the  slice 

rvj 

y. 
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5.2  A  Banach  manifold  of  nearly  spherical  stationary  spacetimes 

As  we  saw  in  §2.4  the  normal  derivatives  of  M  at  the  star 

boundary  can  be  investigated  using  lim(d1M/d(log  T)1)  and  this 

p— >0 

has  a  similar  behaviour  to  the  behaviour  of  lim(d1M/dU1)  in  the 

p— >o 

static  case.  Thus,  using  the  results  of  theorem  2.1  we  can  find 

—S' 

that  for  a  given  equation  of  state  such  that  limp  p=d>0  with 

p— 

(l— /(l +1),  M  has  the  same  differentiability 

properties  as  were  determined  for  the  static  case  in  §4.2.  Thus 

we  again  take  M£Mg_3(5+n  for  a  fixed  s^3  determined  by  the 
equation  of  state  and  some  p>3  (and  less  than  some  upper 
bound  determined  by  the  equation  of  state)  and  any  large  n. 

The  equations  (5.11)  then  show  that  we  can  take  U£Mg_1(5, 

aeXs-i,6+i>  1,<5+1  and  where  we  fix  r)=d p  in 

terms  of  a  polar  coordinate  system  related  to  the  asymptotically 

cartesian  system  in  the  usual  way  so  that  rj  has  the  form  (1.37) 
for  |x| — >°°  and  0^<5<l-3/p  is  fixed  by  the  asymptotic 
conditions. 

Since  the  equation  of  state  p(p)  and  the  surface 
temperature  are  regarded  as  fixed,  for  every  value  of  m  (in 
some  interval  [0,mcrit)  )  there  is  a  unique  spherically  symmetric 
model  (y,U,Ub)  since  Ub  is  also  given.  Since  0=0  in  the  static 
case  we  have  b=0  so  the  first  expression  in  (5.11)  shows  that 


. 


■ 
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7  —  ^s.tf+i  •  From  (5.7)  we  see  a=Tb1e  Ub  so  that  giving  a  is 
equivalent  to  giving  Ub  in  the  static  case.  Thus  for  a  given  m 
we  can  describe  the  unique  spherically  symmetric  solution  by 
cr=(7,U,cx=0,77,a,b=0)  .  Again,  all  of  the  components  of  a  have 
more  differentiability  than  those  of  the  general  <r  ’s  except  for 

U,  so  we  must  again  make  use  of  U  =  U-f0p2  f  to  define  the 
topology,  where  the  constant  f0  is  defined  as  in  the  static  case 
in  terms  of  the  unique  spherical  solution.  (See  the  discussion 

below  equation  (4.3).)  Thus  while  UEMg_1(5. 

In  analogy  to  the  static  case  we  define 

^)s,<5=^a  =  (T.U,a,h.a)b )l7e^s-1,<5+1>  UFM^jU), 

a£Xg_ita+i,  a,be®(  (5.12) 

where  we  will  again  use  U  to  do  the  calculations  but  keep  in 
mind  that  it  is  U  that  defines  the  topology.  Note  that 

rsj  rNJ  p 

O^s-id+i  and  again>  any  two  solutions  o  v  a 2^^s-i,6+i  wil1  have 
the  same  mass,  but  they  do  not  necessarily  have  the  same 
angular  momentum. 

The  action  of  the  dif feomorphism  and  gauge  group  on 

3»ia+i  can  now  be  investigated.  We  will  see  that  in  order  to 
obtain  a  slicing  theorem  along  the  same  lines  as  in  the  static 
case  it  will  be  necessary  to  let  <5  —  0  but  this  is  the  largest  space 
allowing  the  slowest  fall  off  rate  so  we  do  not  lose  anything. 


Let 


^b=®bXMbH0.X)l?£2>b  x€M*  (S3,S)i 


(5.13) 


and  define 


((^.x).(7.0.a,7),a,b)) — >(ip  a  +  dx.lh.a.b) 


(5.14) 


where  p  denotes  the  pull  back  by  the  dif f eomorphism  p  of 

the  appropriate  covariant  tensor  field.  Equip  'S'l 6  with  the 
product  differentiable  structure.  It  also  carries  a  group 

structure  as  the  semidirect  product  of  2)^  and  the  Abelian 

group  with  respect  to  the  natural  action  of  <2)^  on  £  by 

pull  backs. 

Theorem  5.1.  If  p>3  and  s>2+3/p  (and  <5=0  )  then 

P  rsJ  p 

(i)  A  defines  a  continuous  action  of  on  UPs-i.v 

A G  is  C°°  and  if  <7e^_1+kl  then 

A a:(<P,x) — >A ({<P,X)>°)  is  Ck  for  k=0  or  k 

(ii)  If  aG5)g_1+ki  (for  k=0  or  1)  then  the  orbit 
S(7:=|A((^1x)^)I(^,x)^^s,^C^>s-i1i  is  a  °k  submanifold. 

(iii)  If  then  there  exists  a  neighborhood  V  of 

(1,0)  in  and  a  slice  of  the  action,  i.e.  a 

submanifold  S  of  ^_u  containing  cr  such  that 

r\,/ 

((cp,x),v') — >A({p,x),v')  is  a  homeomorphism  of  VXf/ 
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onto  a  neighborhood  W  of  a  in  and  €(a3?=\(7\. 

Proof:  (i)  The  proof  is  analogous  to  theorem  4.1  for  7  and 
U  .  Assume  first  that  5€[0, 1-3/p)  as  in  that  theorem.  Since 

X  ><^X  is  smooth,  (<p,a)  ><p  /  dx1(a°cp)k  is  continuous  by 

theorem  3.5  and  linear  in  a,  so  ((<^,y),cx) — ^  a+dy  is  C  in  a 

and  continuous  in  (<p,x)  ■  If  «^X^(5+1  then  aka£Mg_1(5  and 
* 

(cp,dka) — >cp  dka  is  continuous.  For  the  corresponding  proof  for 

( <p,r] ) — ^7?*  1?7  it  is  necessary  to  set  5  =  0  to  satisfy  the  condition 

o M  m  77  is  bounded,  required  in  theorem  3.5,  because  of  the 
asymptotic  behaviour  of  rj. 

(ii)  Ac  is  injective  since  there  are  no  non— trivial 

isometries  in  3)^?,  and  there  are  no  constants  in  .  If  k=l 

the  tangent  map  TA^  at  cp= I,  y  =  0  is  given  by 

where  fEXg0,  fCM^ 0  .  Cantor’s  version  of  the  Berger-Ebin 
decomposition  (corollary  3.12)  together  with  the  splitting  given 

by  corollary  3.11,  Xg_u=d(Mg0)©/7  where  £7=ju€Xg_uldiv7cu=0} 

imply  that  the  image  of  T(I0)Aa  splits  so  that  Aa  is  a  C1 
imbedding. 

(iii)  Choose  open  neighborhoods  V  of  0  in  Jy  and  W  of  0 
in  $  such  that  7+c  is  positive  defenite  for  cCV  .  Then,  if 
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y=|(7+c)U  +  u)a+/S,7?+C.a+A)b+B)|cev, 

/5eW,  ^£Xg_li0,  A,B€*j 

is  clearly  an  imbedded  submanifold  of  ^p_.  .  .  The 

S  1  j  l 

h o m e o m o r p h i s m  is  then  obtained  by  use  of  the  deco m position 
theorems  as  was  done  in  theorem  4.1.  B 

rsj  T~v  fVJ 

In  particular,  the  tangent  space  to  at.  splits 

so  for  <j  we  have 


T^s-u-T^.+T^.  (5.15) 

We  conjecture  that  the  set  of  solutions  of  Einstein's 

f**SJ  Cso 

equations  Sp  in  some  neighborhood  V  of  a  in  ^s-iv  namely 

^  (O)HV,  is  equal  to  U  & ^fW  where  K  is  a  one— dimensional 

submanifold  of  passing  through  a  and  parameterized  by 

the  angular  momentum.  As  in  the  static  case,  however, 

r^j  c^j 

£\o)=l„£  is  not  surjective  at  o=a  so  we  are  unable  to  show 

r-o  _ 

that  is  a  -submanifold.  Thus  we  cannot  show  that  there 

exist  solutions  with  nonzero  angular  momentum  nor  obtain  as 
strong  a  uniqueness  condition  as  we  would  like.  But  since,  in  a 

rx_f 

neighborhood  of  a,  the  set  of  equivalence  classes  of  with 

respect  to  the  action  of  is  in  one-to-one  correspondence 

with  the  slice  if  through  a  we  show  that  there  can  be  no 
non-constant  C1  curve  in  <?HSp  passing  through  a  particular 


1 


I  f 


. 
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solution  (existence  assumed)  with  some  angular  momentum 
unless  the  curve  passes  through  solutions  with  different  angular 
momenta.  This  will  be  done  by  a  similar  procedure  to  that  used 
in  the  static  case,  solving  the  linearized  equations  on  the 
spherical  background  and  extending  the  result  to  a 
neighborhood  of  cr  by  the  use  of  Cantor's  isomorphism  theorem 
and  Nirenberg  and  Walker’s  theorem.  Since  Nirenberg  and 
Walker’s  theorem  is  proven  by  contradiction  one  can  not 
determine  the  size  of  a  neighborhood  for  which  this  extension  is 

valid  and  therefore  there  is  little  hope  of  obtaining  an  estimate 
for  a  lower  bound  for  the  value  of  J  which  leads  to  a 

bifurcation  of  the  solution  space  such  as  one  has  in  the  classical 
Newtonian  series  when  the  Maclaurin  sequence  bifurcates  to  the 
Jacobi  ellipsoids. 


5.3  Linearization  of  £'  to  find  ker  £'(o) 

Let  (c=<57,u=(5U)/?=(5c(,('=(5?7;A=(5a)B=(5b)GT0.5)3_1  v  Let  u=oU 
and  recall  that  u  and  u  are  related  through  equation  (4.9) 
with  u  being  less  differentiable  than  u  .  Again  we  do  all 
calculations  with  u  keeping  this  in  mind.  The  variation  of  the 
pressure  can  be  seen  to  be  given  by 

dp=(p+p)[-T(l+2e~2UT202)u+^T3e~2U0r0s$rs— T2u/?r0r 
-T2(Dar-Te“SUer)(B77r+bfr)-T4A] 


For  o’  =  c7^^>s,i 


we  have  from  the  proof  of  theorem  5.1  the 


■ 
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decomposition 


c=^7+<i>,  div7(<3>)=0  (5.16) 

u=o^U+*  (5.17) 

|5=d£^a+df-t-o,  div7(co)  =  0  (5.18) 


+v 


(5.19) 


where  ££X^0,  $€S*_U,  ^£MSP_U>  wexj_u,  fEMj0l  ^€X^_li0  are 

all  unique.  Another  lengthy  but  standard  calculation  gives  that 


£'{<7)\la9 


p 

s— 1,5+1 


3,3®  ^  s— 3,3®  ^  s— 3,3®  ^s— 3,2 


(5.20) 


is  given  by 


£\{<J  )(c,u,^,^',A,B)ij= 


n  ,S  _  — 2U ;r  ,s 

8(i(5jr2pe  di0j 


iVrVrcii+TV1VjcPv(i(Vrcj)r)+[-3E^J)+Rrs7ii+iRq6Js 

+  (RirtR7.rth.hje457rS+e4Uhrh(i5p2(p+p)e“4UT+r 

-(p+p)T202e_4U6pp^e_6U(p+p)T20r95(4e2U7ij+TS(p'+3)(e27ireie  ,))]c„ 
+4U(13J)u+[2hihje4U+2(p+3p)e"2U7ij+(p+p)e_4UT!!(8eZ7ir4e10 
+  T(p'+3)(027ir0,0J))]u  +  h(1£j)”e4UVr/3s+-yT(p  +  p)e~4U[2 

+T2(p'+3)(e27ii-e1ej)]e7r+(p+p)^4UT2[27k(iei)+2(l'T"1 

+(p'+3)(^Tc<k-T2e"2U0k)(e27ireiei)](Bpk+bfk) 

+  vT2e^4U(p+p)[2e2UT“17ii+(p'+3)T(02-010j)]A, 


1J 
,  2U 

Je  7 


.  2U  o/!  \ 
e  «k-20k)7ij 


(5.21) 
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^'s(<r)(c,u, /?,<•, A, B)= 

Au+[2h2e4U+(p+3p)e-SU+(p  +  p)e-2U(4T2e2e-2U+i(p'+3)T2)]u 
-Ur(VSosr-iVrc:)+tihrhse4U-VrVsU-ih27''s 

~(p+p)T  e  (l+7r(p'+3)r)erea]crs+Hrse4UVr/?s+-2-(p+p)(p'+3)TTtie  "l0rflr 
+(p+p)e_!iUT[-2Te'2Uer+i(p'+3)(^c<r-e“2UTer)T](B77r+bfr) 

+7-(p  +  p)(p'+3)Tt>re  2l'A,  (5.22) 


if'g((j)(c,u,^,f,A,B)i= 


^VrVr(3  -^V1(Vr^r)+4UrV[r)Si]-[^R[+(p+p)e"6UT2(t)2(p'+3)-e£U)eier]/?r 


4U„r,TTs  4U-, 


+(p+p)uTe  6U((5^s+^e  2U(p'+3)T  rdru 

-(p-fp)uTe  6L(4+T(p'+3)0iu+(p+p)e  6UT2[uT  1yirJc^'Xioir6i 


+V(p'+3)(Te"2Uer-^ar)0i](Bpr+bCr)-(p+p)e"6V((p'+3)^-e2U)0iAf 

(5.23) 


c^'4(a)(c,u,/?,t,A,B)i= 


VrV(rC i)+ 2"^rVr(Vsc^)— •^-Vrpi(Vsc^)+V  rps  (Vrcsi— yVicrs)+^:(Vi?7r+Vrpi)VrCg 

(5.24) 


,r  s  k-. 


+p?W+7R3>iHKvrvVRYkP>rS. 


where  p 


_ 2 1 J  2  2 

dp/dp  and  T:=l+2e  T  6  . 

r*j 


As  in  the  static  case  £'{&)  is  not  elliptic  but  by 


restricting 


■ 
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£'(cr)(c,VL,p,t,k,B)=0 


(5.25) 


to  the  slice  we  can  obtain  an  elliptic  operator  to  which  we 
will  be  able  to  apply  Cantor’s  isomorphism  and  Nirenberg  and 
Walker’s  theorems.  To  this  end  define  (for  x:=(c,u ,/?,£")) 

L1(cr)(x)ij:  =  2,1(cj)(c)u)/?,t,A,B)ij+Y(K7odiv7C)ij 

-^T2e  4U(p  +  p)[2e2UT  17ij+(p'+3)T(0271j-0i0J)]A 
-(p+p)e“4UT2[27k(iej)+2(^T_1e2UCxk-20k)7ij 
+(p'+3)(nTak-e“2UT2@k)(027irei0j)]?7kB,  (5.26) 

L2(o)(x):  =  2'2(ct)(c,u,  /?,<•,  A,  B)  +  Ur(div7c)r 

—  y(p  +  p)(p'+3)TuTe  2UA 

—  (p  +  p)e  2UT[-2Te  0r+y(p'  +  3)(nar— e  T@r)T]77  B,  (5.27) 

L3(<j)(x)  3(cf){c,u,(3  ,k,B)  L+^d  yfi) 

+(p+p)e-6UT2[(p'+3)r>2-e2U]01A 
-(p  +  P)e-6UT2[^T-17ir+e2Uarei+ 
i;(p'+3)(Te“2U0r-^ar)0i]prB, 

L4(^)(x)i:=2'4((j)(cIu,iS,C,A,B)i-i7?rVr(div7c)i. 

rso  .  .  ro  rso  Cs-jl  \  /  \  •  l 

It  is  easy  to  verify  that  L(<r)=(  L1,L2lL3,L4)(o )  is  a  second 
order  elliptic  operator.  Equation  (5.25)  is  now  equivalent  to 


(5.28) 

(5.29) 


. 
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L1(^)(x)-^K7odiv7c=q1(a)(A,B)T 

(5.30) 

L2(a)(x)-VUJdiv7c=q2(o-)(A,B)T 

(5.31) 

L3(^)(x)-Yd(divr/?)=q3((7)(A.B)T 

(5.32) 

L4(cr)(x)+Y?7Jd(div7c)=q4(a)(A,B)T 

(5.33) 

where  the  qi(a)  are  the  obvious  row  matrices  of  tensor  fields 

obtained  from  equations  (5.26-5.30)  (  q4=(0,0)  )  and  (A,B)T  is  the 
transpose  of  (A,B)  .  We  also  write  q(u)  for 
(q1(CT),q2(a),q3(cj),q4((7)). 

We  will  show  that  for  aCtf  near  a  and  x=(cJu,|f?I£)ESg_1 1 
©Mg_ii©Xg_i i©Xg_^Q  the  inhomogeneous  equation 

L(a)(x)  =  q(o-)(A,B)T  (5.34) 

implies  that  c=u=jS=C=A  =  B=0  if  also  the  variation  <5J  of  the 
angular  momentum  corresponding  to  (c,u,/S,<f ,A,B)  vanishes.  This 
will  then  imply  that  the  solution  of  (5.25)  vanishes  if 
(c,u,/?,£ ,A,B)  is  tangent  to  the  slice  and  the  corresponding  <5J  =  0. 

We  will  do  this  in  the  same  manner  as  in  the  static  case, 
solving  (5.34)  first,  on  the  spherical  backgrouund  cj  and  then 
extending  the  result  using  the  theorems  3.7  and  3.14  about 
elliptic  operators. 

On  the  spherical  background  equation  (5.34)  simplifies 
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greatly,  decoupling  into 


Lli2(o-)(c,u)  =  q12(CT)(A,0)T 


(5.35) 


L3(^)(/3,^)  =  q3(cT)(0,B)T 


(5.36) 


L4(ct)(c,<^)=0 


(5.37) 


where  (5.35)  is  the  set  of  equations  (4.16)  obtained  in  the  static 
case.  Theorem  4.4  then  gives  us  that  c,  u  and  A  vanish. 

Putting  c=0  in  (5.37)  yields  the  equation  div  °K  <T  which  bv 

7  7 0 

corollary  3.9  implies  that  £  =  0.  Letting  equation 

(5.36)  then  becomes 


£Ukai(e4U/ck)=2(p+p)T2aB?7i 


(5.38) 


with  /c1  satisfying 


(5.39) 


These  equations  are  no  easier  to  solve  in  this  form  than 
(5.36)  is.  However,  if  we  write  these  equations  in  the 


(2+1)— dimensional  formalism  we  obtain,  with  k?:=k\J.v 


'j7a=Va  (A=l,2)  and  D  the  derivative  defined  by  (1.46),  which  by 
the  choice  of  coordinates  which  makes  71A=0  can  be  taken  to  be 


0 
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D/ca+4/ca-'W  2<9a/c°=2W  1(p+p)e  4UT’“aB  a  AB?7 B  (5.40) 

^[a/cb]=®  (5.41) 

D/c0-W"1/c°DW-W“17cAaAW  +  VATPA+W"1n/<c0=0  (5.42) 

where  all  the  tensor  fields  on  the  U  =  const.  surfaces  are 
obtained  from  the  unique  spherical  solution  a  .  As  is  seen  in 
appendix  3  these  equations  are  the  same  as  some  of  the 
equations  obtained  from  linearizing  the  field  equations 
(1.48—1.58)  in  the  (2+1)— dimensional  formalism  keeping  the 
central  potential  Uc  and  pressure  pc  fixed  (which  by  lemma  4.3 
we  know  is  the  case  here  as  well).  The  analysis  of  (5.40  —  5.42) 
which  we  give  below  was  first  given  by  Kunzle  and  Savage 
(1980a). 

_  A  —1/  \ 

Since  the  1-form  ka dx  is  closed  on  the  SC=U  (c)  spheres 
by  (5.41)  and  since  the  first  Betti  number  of  S  is  zero  (c.f. 
Goldberg  (1962)  p.89)  7cAdxA  is  exact  so  that  there  exists  a 
function  K(U,xA)  on  Sc  determined  up  to  an  additive  function 

of  U  only  such  that  70A=9AK  .  Letting  k:=/c°,  equations  (5.40) 
and  (5.44)  now  become 

aA(DK+4K-W_2k)=2TceUc_6U?r1(p  +  p)TAB8B  (5.43) 


AK=-Dk-W  \2n-MW  *)k 


(5.44) 
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since  T=TceU°U. 
r2(U)(dtf2+sin26>d/) 


To  solve  for  K  and  k  we  let  7ABdxAdxB= 
.  Now  expand  into  spherical  harmonics, 


K=  2  K{m(r)Y(m(e,«,),  k=  E  kc  (r)Yg 

t,m 

where  jj  Y£mY£,£,sin  dddd^=(5££,<5mm. .  Let  R'=D(log  r)  .  Equation 
(5.44)  then  becomes  0=r2Dk0+rSW_1(2O-MW_1)k0  for  1= 0  and 


Kem=r2(«C+l))“1[Eklm+(4R'-M¥  2)k(m]  (£*0).  (5.45) 


Since  v=d  (5.43)  yields  E  (DK(m+4KJm-W  2k(m)39Y(m= 

l,  m 

2BTceU'"6lJW_1(p+p)r2sin  6  and  (E  (DK(m+4K(m-W"2k(m)d?iYIm=0  so 

that 

DK^rn+^K^-W  2k^m=0  for  m^O  or  f^0,l  (5.46) 

and 

DK1+4K1-¥_2k1=-2B(47r/3)1/2TceUo_6UW“1(p+p)r2  (5.47) 

Note  that  K0  is  an  arbitrary  function  of  U  and  need  not  be 

A 

known  to  determine  /cAdx  uniquely. 

With  the  use  of  equations  (1.48)  and  (1.53)  on  the  spherical 

—2 

background  we  have,  from  the  equation  for  k0,  k0=dWr  for 


d  =  constant  .  By  doing  an  expansion  in  normal  coordinates  (y1) 
at  the  center  we  find  that  W  =  ^-M0|y|  +  O(|y|3)  and  r=|y|+0(|y|3)  so 


■ 
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that  in  order  for  k0  to  vanish  at  the  center  (since  h°=h191U=0 
at  the  center),  c=0,  whence  ko=0.  Equation  (5.46),  with  the 
use  of  equation  (5.45)  to  eliminate  Kgm  and  equations  (1.48)  and 
(1.53)  on  the  spherical  background  become 

w  DDk,m+MDkem-F,(U)k,m=0  for  m*0  or  1*0,1  (5.48) 

where  kfm=rW-1e2Ukgm  and 

Fe(U)=W2(5— 4R')+£(£+l)r~2+2M-pe-2U  (5.49) 

Since  the  three-dimensional  Laplacian  A  in  the  spherically 

symmetric  case  has  the  form  Af =WSDDf +MDf+  A f,  equation  (5.48) 
is  equivalent  to 


Akgm=Fg(U)kgm  for  m^O  or  1^0,1  (5.50) 

We  now  use  an  argument  similar  to  that  used  by  Klinzle 

g  — 2U 

(1971)  in  the  static  case  (c.f.  appendix  2).  Let  u=fr,  v=pr  e 
and  x=Mr  /(3u)  so  that  u,v,x^0  in  the  physical  domain,  and 
ur=const.  and  v  =  x  =  0  in  vacuo.  We  can  show  that  Fg(LT)^0  is 
equivalent  to  £2(£+l)2-t9u4+36u2x2+v2-t2(5£(£+l)  — 8)u2+12£(f+l)ux 
-2£(£+1)v+60u3x-26u2v-12uxv^0  .  Since  p,p^0,  we  have  that 

2  4 

0^v^2xu  so  that  the  above  inequality  is  implied  by  l  (£+l)+9u 

+12u2x2+2(5£(£+1)-8)u2+8u3x+v2^0  .  Since  5f(f+l)-8^2  for  f^l 
this  shows  that  Fg(U)>0  for  1^1  .  The  asymptotic  flatness 
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conditions  "with  Einstein’s  vacuum  equations  give  us  that 
U— mr  +  f,  ftMg_1(1  and  W=mr  +  g,  gEMg_22  in  the  asymptotic 
region.  Together  with  the  differentiability  properties  /c1CXg_11  it 

is  then  easily  seen  that  .  Thus  hm^supkfm=0  .  Also 

k(Uc)=(/c  Ui)(Uc)=0,  so  k^m=0  at  the  center.  These  boundary 
conditions  together  with  (5.50)  and  Fe(U)>0  for  £>1  yield  k£m=0 

for  m^O  or  £^0,1  since  k^m  is  C1  and  we  can  apply  a 
maximum  argument. 

Thus  the  only  nonvanishing  component  is  k1=k10  which.,  by 
(5.47)  and  (5.45),  satisfies 

W2DDk1+MDk1— F1(U)k1— — dcre  4U(p+p)  (5.51) 

where  dc=4-BV/47r3  ly/2TceU° .  In  vacuo  one  can  easily  show,  using 

-1  2 

the  exterior  Schwarzschild  solution,  that  W— m  sinh  U, 

r=— m  sinh  U  and  R'=-coth  U  .  The  only  solution  of  (5.51) 
satisfying  the  asymptotic  conditions  is  then  found  to  be 


k=4CeUsinh2U.  (5.52) 

(C  is  determined  by  B  through  the  interior  solutions  discussed 
below.)  This  agrees  with  what  one  obtains  by  linearizing  the  Kerr 
solution  on  a  spherical  background,  in  this  formalism.  (Hartle 
and  Thorne  (1968)  have  shown  that  the  empty  space  metric 
outside  any  slowly  rotating  perfect  fluid  agrees  with  the  Kerr 


. 
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metric  up  to  second  order  in  the  angular  velocity.) 

To  determine  the  solutions  in  general  we  introduce  a 

coordinate  z  defined  by  e  =e  c+r  z  where  r>0  is  chosen  such 
that  z(surf ace  of  star)  =  l  .  Note  that  z  behaves  like  a  radial 
polar  coordinate  near  the  center  since  U  has  a  positive  definite 
critical  point.  Equation  (5.51)  can  then  be  written  as 

ud2k/dz2+Trdk/dz+tk  +  det=0  (5.53) 

where  n=e2UW2,  rr  =  2T2ze2UM  — z  ^e21^— r4 z4)W2,  l——  4t4z2(2M  —  p 

— 4W2R'+5W2+2r  2)  and  t  =  8T4z2e  2Ur(M  — p)  .  Note  that  since  s^3, 

n^Mg_22CC°,  7T,pt€-Mg_31  and  K^Mg_20CC°.  As  we  noted  before,  we 

assume  that  the  equation  of  state  is  analytic  in  order  to  do  the 
following  expansions  near  the  center  in  terms  of  z  .  After  a 
lengthy  calculation  (and  dividing  by  a  common  factor)  we  obtain 

i/=z2+u4z4+0(z6)  7r=2z4-7T3z3+0(z5) 

t=-2+0(z2)  t=t3z3+t5z5+0(z f)  (5.54) 

where  the  coefficients  depend  on  Uc,  pc,  pc,  dp/dp(z=0),  etc. 
Considering  the  homogeneous  equation  to  (5.53),  we  see  that  z=0 
is  the  only  singularity  in  [0,1]  and  it  is  a  regular  singular 
point,  so  if  we  make  suitable  power  series  expansions  we  find 

that  the  general  solution  of  (5.53)  satisfying  the  regularity 
conditions  at  the  center  is 


, 
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k=/j.kh+Xki  (5.55) 

where  kh  is  the  solution  of  the  homogeneous  equation  with 
kh(0)=0,  dkh/dz(0)=l  and  kA  the  solution  of  (5.53)  with  ki(0)=0, 
dki/dz(0)=0  .  To  determine  X  and  fj,  we  use  the  fact  that 
dk/dtl  is  continuous  at  the  star  boundary  together  with  the 

boundary  conditions  k(l)=Ce  Ub(l-e2Ub)2=:C1,  d.k/dz(l) 

— — 2T2Ce  2'b(l-e2Ub)(l-f3e2lb)=:C2  (where  eUb—  eUc+r2  .  Thus  X  and  fi 
are  determined  by  /r.kh(l)  +  Xki(l)  =  C1  and  yudkh/dz(l)  +  Xdki/'dz(l) 
=C2  .  These  equations  will  uniquely  determine  /r  and  X 

provided  W  (1)^0  where  <M^(z)=khdki/dz-kidkh/dz.  Now,  W( z) 
obeys,  by  (5.53), 


vdW /dz+nW+tk^O  (5.56) 

with  W  (0)=0  . 

A  power  series  expansion  gives  a  unique  solution  for  (5.56), 

W=— z3/5+0(z5)  where  the  factor  t3  has  been  absorbed  into  dc 
which  is  just  a  fixed  constant.  Thus  for  small  z>0,  ^<0  and 

dW/dz<0  .  Now,  n,t>0  in  (0,1]  (see  (5.54))  and  kh(0)>0  in  (0,1] 
since  kh(0)  =  0,  dkh/dz(0)>0  and  if  z1G(0,l]  is  the  smallest  z 

such  that  dkh/dz(z1)=0,  then  dckh/dzc(z1)=-(i/i/)kh(z1)^0  since 

<,=-4t4z2F1(U)^0.  Thus  kh  has  no  maximum  in  (0,1]  .  Now 
suppose  z2  is  the  smallest  z  61  (0,1]  such  that  it  (z2)  — 0  .  Then 


139 


dtf/dz(z2)=-t(z2)k1(z2)y  1(z2)<0  so  that  there  exists  an  £>0  such 
that  14'  (z)>0  for  z£(z2-e,z2)  which  contradicts  W< 0,  dli^/dz-CO 
for  s m all  z  .  Therefore  W (1)^0  and  there  exist  unique  X  and 
/r  solving  the  equations  giving  the  boundary  conditions  of  k  . 
Thus  there  is  a  unique  k  determined  by  B  with  an  asymptotic 
form  given  by  (5.52)  and  vanishing  identically  if  B  =  0  .  From 
this  it  is  readily  verified  that  there  is  a  unique  solution  hc  to 
(5.38)  and  (5.39)  which  must  have  the  asymptotic  form 


/c^C/cm  3+0(r  4)  (5.57) 

and  which  vanishes  iff  B  =  0  (iff  C=0  ). 

We  can  now  use  this  to  investigate  the  solution  /?  of 

(5.36).  (Recall  that  <f=0  .)  Since  Ta^-i,i  sPllts-  i-e- 

(equation  (5.15)),  we  can  investigate  jS  on  the 

two  subspaces  separately.  Since  ^a-tdf  on  T JO  a  for  some 

^Xj0,  f£M^0  and  2(A(^x)(ct))=^*(2(ct)),  (5.32)  implies 

d(div  odf)  =  0  whence  f  =  0  by  corollary  3.8.  On  T^Y  we  have 

div^j5=0  so  that  (3  is  uniquely  determined  by  ic  and  hence  by 
B  .  This  follows  from  corollary  3.9  since  tc= 0  implies  B=0,  d/5  =  0 

so  that  equation  (5.36)  becomes  0=VrVrjffi-Rjr/?r=-(div?°K7j5?)i  . 

From  (5.57)  it  is  seen  that  /?  has  an  asymptotic  expansion  of 

0(|xP2)  at  infinity  and  vanishes  iff  B=0. 

Taking  a  variation  of  (1.39)  it  is  easily  seen  that  the 


140 


coi  responding  variation  in  angular  momentum  is 

5J=(l/l6rr)|xlm^/]x|=oonsta[l/?j]T)irlkdEli  (5.58) 

and  so  will  not  vanish  unless  /3=B=0  .  In  fact  k  has  an 
asymptotic  expansion 

k,  =2dJ|x|  3((5q  3 1 x I  ~x1x3)+0(|xl  4).  (5.59) 

Summarizing  the  above,  we  have  the  following  theorem. 
Theorem  5.2:  The  operator  equation 


L(cr)(c,u,/?,(f)=q(cT)(A,B)T 

r-u 

with  (c,u,j8,^,A,B)£To.jPg_1 1  implies  that  c,  u,  and 

A  vanish  and  that  /3,  B  are  uniquely  determined  by 
<5J,  vanishing  iff  dJ  =  0. 

5.4  Curves  of  solutions  in 

T 

Again,  the  presence  of  the  inhomogeneous  terms  q(cr)(A,B) 
does  not  allow  a  direct  application  of  theorem  3.14  so  we  proceed 
analogously  to  the  static  case  with  the  slight  complication 
introduced  by  having  two  constants  A  and  B.  From  theorem  3.7 

it  is  easily  verified  that  for  cr<E^_unSp  (so,  in  particular,  for  a) 
the  map 


L(<t  ):P:=Sg_lj0©Mg_lj0©Xg_li0©Xg_li0' 


"^s-3,2@Ms-3,20^s-3,2®^s-3,2~:Q 


is  a  continuous  linear  operator  with  finite  dimensional  kernel 
and  closed  range.  Just  as  in  the  static  case  this  follows  using 
properties  (3.6— 3.8)  of  the  weighted  Sobolev  spaces  together  with 

the  asymptotic  properties  of  crCS^. 

Equation  (5.35)  determines  c  and  u  uniquely  in  terms  of 

dm  for  xEP  while  (5.36)  determines  uniquely  in  terms 

of  <5J  and  c  since  (5.37)  determines  uniquely  in  terms  of  c, 

as  seen  above.  Thus  the  solution  space  of  L(a)x-q(cr)(A,B)T=0  in 
PetR©tR=|(x,A,B)(  is  2— dimensional  and  spanned  by  (x0,0,B0)  and 
(y0,A0,0)  corresponding  to  the  solutions  with  (5m,dJ)  =  (0,J)  and 
(m,0)  respectively. 

For  a=o,  putting  A=0  means  dm=0  so  c  and  u  vanish, 
as  well  as  ^  by  (5.37).  Putting  B=0  gives  [3=0  .  For  o  in 

rsj 

some  (small)  neighborhood  of  <j  in  the  operator 

r^J  ,  .  r^j  r^j  .  ro.  . 

L(<r):P — >Q  is  injective  since  L(cj)  and  L(<r)  then  satisfy  the 
hypothesis  of  theorem  3.14. 

Thus  there  are  unique  x0,y0^P  such  that 


L(cr)x0=q((j)(0>l)T 

(5.60) 

L(cj)y0=q(cr)(hO)T 

(5.61) 

In  fact  x0=(0,0,[3  0,0)  where  (30=  0(1x1  2)  at  infinity  while  y0G:P\Pe 

where  P^S^  1©Mg_11©Xg_u©Xg_1)0  is  a  subspace  of  faster  fall  off. 
Suppose  a  is  in  a  neighborhood  of  o  such  that 


. 


' 
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I!l(c7)-L(ct)II<s1  and  l!q(cT)-q(a),iPiS_3i(5+3<£2  (5.62) 

for  some  small  ev  e2>0  .  Then  by  theorem  3.14  kerpL(a)=|0( 
and  there  are  unique  solutions  x,y3P  of 

L(cr)x  =  q(a)(0,l)T  (5.63) 

and 

L(tf)y=q(cr)(l,0)T  (5.64) 

r^i  r~sj 

Using  the  fact  that  L(cf):P — >Q  is  injective  and  has  closed  range 
and  so  has  a  bounded  inverse  we  can  show,  exactly  as  was  done 
in  the  static  case,  that  x  is  arbitrarily  close  to  x0  and  y  is 
arbitrarily  close  to  y0  for  small  enough  and  e2.  Therefore 

y  cannot  lie  in  P  and  x  =  (c,u ,/?,<f)  must  have  /S=0(|x|  ~)  at 
infinity  and  thus  corresponds  to  a  solution  with  <5J^0. 

By  the  remark  below  (5.34)  we  then  have 

Theorem  5.3:  If  c:[0,l] — Sp(-p)  is  a  C1  curve  of 
solutions  having  all  the  same  (small)  angular 
momentum  J  then  c  is  constant  if  the  slice  tf  is 
contained  in  a  small  enough  neighborhood  of  <j. 

Proof:  Exactly  as  in  the  static  case. 

Note  that  since  Nirenberg  and  Walker’s  theorem  3.14  is 
proven  by  contradiction  this  method  cannot  be  used  to 
investigate  how  large  J  can  become  before  there  is  a 
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bifurcation  of  the  solution  space. 

5.5  The  problem  of  surjectivity 

Whether  there  is  a  physical  (non— singular)  source  for  the 
Kerr  solution  is  a  long  outstanding  question.  (See  Krasinski 
(1978)  for  a  recent  review  of  investigations  of  this  question.) 

More  particularly,  whether  a  rigidly  rotating  perfect  fluid  could 
be  the  source  is  unknown.  Indeed,  it  has  not  even  been  shown 
that  there  exist  rigidly  rotating  perfect  fluid  solutions  of 
Einstein’s  equations.  We  would  like  to  be  able  to  answer  this 

rsj 

latter  question  by  showing  that  £(cr)  is  surjective,  for  then  we 

would  have  shown,  by  theorem  3.2,  that  £  (0)  is  a  submanifold 
of  dimension  one  in  a  neighborhood  of  a.  (We  would  then  also 

rsj 

have  that  £  was  linearization  stable  at  a  and  the  uniqueness 
theorems  for  both  the  static  and  stationary  cases  would  be 
stronger.) 

r*j 

Unfortunately,  neither  £'{<j)  nor  £'(&)  are  surjective.  We 

rsj 

will  show  this  only  for  £'(a)  since  the  argument  for  £'(a)  is 
exactly  analogous  and  can  be  obtained  by  dropping  the  tildes 
and  using  equation  (4.13)  instead  of  equations  (5.26-5.29)  to 
define  the  map 


2(a) 


:S 


p 

S  1,1 


>2 


(c,/?)1 — KiK7°div7c,  VUJdivyc,  ^d(div^),  ^-7?Jd(div7c)) 


so  that 


. 

•>'*y 

' 
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rvj 

«^(o‘)(c)u,/S,^,A,B):  =  L(o')(c,u,^,(')— q(a)(A,B)T=c^,'(cr)(c,u,/S,<f,A,B)+^(cr)(c,/S) 


for  (c,u,p,^, A,B)£Tgt^)g_ltl  .  For  0  =  <j  we  have  the  splitting 

D  rs-'  r^' 

s-i,r^/'a+^(i^  and  our  above  anal37sis  has  shown  that 

t^(cr)iT£7Cp=^(CJ)lssp_u\J7©x^_u\/7  while  ^(cr)|T^=2'(cr)restirctT^  .  Thus, 
if  £'{cr)  is  surjective  so  is  33 (or)  .  Now  ker  S(crH(0,0,j3(B),0,0,B)( 


is  one  dimensional  so  any 


Xj— (i£^7,d£^UJ<£^a+d^,i£^?7,O,O)G!To.0o.  with 


£  nonzero  is  clearly  not  in  ker  §?(<j)  and  S(n)x1=y1^0  .  If  2'(u) 


is  surjective  there  exists  x2=(£,u,/3,<f, A,B)£T  ^  such  that 

«®(a)x2=y1  so  S(cr)(x1-x2)=0  .  But  then  £^y-l=0  with  div7f=0 
which  by  the  Berger-Ebin  type  decomposition  implies  i— 0,  £=0 

which  contradicts  our  assumption. 

So  what  can  be  done?  The  basic  problem  is  clearly  that 


ker  £'{u)  is  too  large  to  allow  surjectivity  onto  the  entire  space 

rsj  rw 

J2,  or  from  another  point  of  view,  3  is  too  large.  But  it  is  by 


no  means  clear  how  to  pick  a  submanifold  of  3  of  the  right 
size.  In  fact  has  too  a  fa^  rate  to  even  be  able  to 

r*j 

apply  Cantor’s  isomorphism  theorem  3.7  to  show  that  S(cr)  and 

thus  £'{a)  have  closed  range.  However,  this  may  not  be  a 
fundamental  block  as  McCowen  (1979)  has  shown  that  the 
Laplacian  on  flat  space  En  has  a  finite  dimensional  kernel  and 
closed  range  as  long  as  <5^n-2-n/p+k  for  some  kCN  so  it  may 
be  possible  to  obtain  a  similar  generalization  to  elliptic 
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operators  with  variable  coefficients.  (Cantor’s  theorem  3.7  is  a 
generalization  of  the  above  result  for  0^(5<n— 2  — n/p  .)  McCowen  is 
presently  working  on  this  type  of  generalization.  Having  a 
closed  range,  however,  does  not  guarantee  that  one  can  find  an 

appropriate  submanifold  of  d2  onto  which  i^'(cr)  will  be 
surjective  but  it  does  suggest  there  is  hope. 

If  such  a  submanifold  exists  it  is  clear  that  there  must  be 

a  splitting  of  the  tangent  space  T02=2  (recall  2(cr)=0  )  which 
could  possibly  be  obtained  from  determining  the  kernel  of  the 

rsj  * 

adjoint  of  ker  £'{(j)  ,  since  we  could  then  write,  for  any 

yGc2,  y=^'(cj)x+yg  with  y2dker  where  x  might  then  be 

r^j  *  rsj  ^  r^j 

uniquely  determined  since  then  £'{a)  y=£'(u)  i£'(cr)x  and  the 
operator  on  the  right  hand  side  is  formally  elliptic  and  self 
adjoint  and  so  could  be  expected  to  be  invertible.  However,  it  is 

a  4th  order  differential  operator  acting  on  Mg_33  tensor  spaces 
with  s^3  so  neither  the  differentiability  nor  the  fall  off  are 
suitable  for  applying  Cantor’s  isomorphism  theorem  (and  in  fact 

PP  PQ  V 

one  is  mapping  into  some  distribution  space  for  small  s  ). 

Only  for  large  s  do  we  even  have  the  use  of  some  of  the 
decomposition  theorems,  so  there  seems  little  hope  of 
determining  an  appropriate  submanifold  in  this  fashion. 

Another  possible  avenue  is  to  make  an  initial  coordinate 

restriction  such  as  restricting  to  the  set  of  metrics  which  are 
such  that  the  Euclidean  coordinates  are  harmonic  and  then 
hoping  to  find  an  appropriate  image  space  on  which  one  could 
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show  surjectivity.  For  metrics  which  are  close  to  Euclidean  this 

set  is  a  subset  of  5^^  but  it  is  not  known  if  this  will  be  true 
far  away  from  flat  space  and  it  is  not  known  if  this  will  lead  to 
any  simplification  in  trying  to  obtain  a  suitable  Banach  space 

onto  which  £'(&)  will  be  surjective. 

In  short,  the  problems  of  existence  and  obtaining  stronger 
uniqueness  theorems  are  still  open  and  require  much  work. 
However,  the  uniqueness  theorems  we  have  shown  are  of  the 
same  rigor  as  Carter’s  stationary  black  hole  uniqueness  result 

and,  since  we  have  shown  is  not  surjective,  improvements 

on  them  will  take  some  additional  insight  into  the  appropriate 
spaces  to  use  for  modelling  stationary  equilibrium  relativistic 
fluids. 
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APPENDIX  1:  A  LEMMA  OF  AVEZ 


The  lemma  we  prove  here  is  used  to  study  a  differential 
operator  on  the  compact  topological  2  — spheres  U=const.  that 
arises  in  the  investigation  of  the  linearized  (2+1)  — dimensional 
form  of  Einstein’s  equations  (c.f.  appendix  2).  Since  all 
discontinuities  in  derivatives  of  the  tensor  fields  are  going  to 
arise  only  across  p=const.  surfaces  which  coincide  in  the 
spherical  case  with  the  U  =  const.  surfaces  we  can  assume  that 

CO 

the  U=const.  surfaces  are  C  manifolds  and  that  the  tensor 

CO 

fields  on  them  are  also  C  .  The  latter  condition  could  be 
relaxed  in  the  proof  we  give  below. 

Let  (N,g)  be  a  2-dimensional  C  Riemannian  manifold,  A1 

the  fibre  of  1  — forms  and  S2  (S2)  the  fiber  of  symmetric 
covariant  2  — tensors  (with  null  trace)  on  N.  The  following  lemma 
was  given  by  Avez  (1977)  (c.f.  Klinzle  (1971)). 

Lemma  Al.l:  If  K<EC°°(S2(N))  satisfies  VBKBA=0  and 

KA=C=const.  then  there  exists  a  harmonic  1-form  0 

such  that  KAB=-^(C-02)gAB+0A0B  where  02=0A0A. 

This  lemma  is  equivalent  to  the  following. 

Lemma  A1.2:  KGC°0(S2(N)),  ^BKBA=0  implies  that  there 
exists  a  harmonic  1-form  0  such  that 

K ab=  ^ a^b  "2”^  Sab-  (Al-l) 
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Proof:  (c.f.  Kiinzle  (1971))  There  exists  a  function  f  and  a 
1— form  0  such  that  KAB=f  gAB+0A0B-  Clearly  f  =  ^cp2.  Now 

V  KAB=20AVj-A(/>Bjh(VA^A)<5!)B=O.  Contracting  with  0B  implies  VA0A=O 

A 

so  0A  is  coclosed.  Then  0  Vj-A0Bj=O.  But  locally  there  exists  a 

1  2 

coordinate  system  such  that  0  =0,  0  ^0.  Hence  V^201j=O  or 

d0  =  O  so  0  is  closed  and  thus  harmonic.  ■ 

Kiinzle  applied  lemma  Al.l  to  the  case  N  =  S  where  0=0 
because  there  are  no  harmonic  1  — forms  on  the  sphere.  However, 
the  proof  given  is  not  valid  because  0,  defined  by  (Al.l),  is  not 
differentiable  at  points  xC!N  where  K  (or  0  )  vanish. 

The  study  of  the  differential  operator 

D:C”(s"(N))-*D"(A1(N)):KA^ViKiB 

and  the  more  restricted  lemma  (which  is  still  sufficient  for  our 
purposes)  which  we  give  below  was  first  presented  by  Kiinzle  and 
Savage  (1978). 

Lemma  A1.3:  If  N  is  connected  and  is  not  conformally 

equivalent  to  a  complete  manifold  of  constant 

negative  curvature,  the  harmonic  1— form  of  lemma 
00 

A1.2  is  C  . 

Proof:  A  connected  Riemannian  manifold  of  dimension  2  is 
conformally  equivalent  to  a  complete  manifold  of  constant 
curvature  (Wolf  (1977)  p.83).  It  is  thus  sufficient  to  prove  the 
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lemma  on  the  surfaces  of  constant  non  — negative  curvature  and 
to  show  that  the  kernel  of  D  is  invariant  under  conformal 
transformations. 

Let  f:(N,g) — >(N,g)  be  a  diffeomorphism  such  that 
f  (g)=X  g,  xec  (N,a),  A(x)t^0  for  all  x£N.  If  KGC°°(Sg(N))  one 

finds  that  K=f^(K)eC°°(Sg(N))  and  ( VBKBA)(f(x))  =  X_2(x)( VEKBA)(x) 
using  the  fact  that  dim  N  =  2  and  L  Ac=TAc+2<5AB3c)log  X 
~6b cS  dDloS  Suppose  that  (N,g)  has  a  constant  curvature 

r*j 

R^O.  By  lemma  A1.2  we  can  write 

nsj  <  r^j  P 

^ab~^>a^>b_2" 0  Sab  (A1.2) 


where  0  is  a  harmonic,  differentiable  1-form  on 
NHxGN|K(x)*0}.  Also,  0  =  0  on  N\N'  and,  in  particular,  on  the 

r^j 

boundary  3N'. 

On  N'  one  finds  that  A(0~)  =  2(V A0BV  A0B+ JrR0~)AO  since 

r\j  rsj  rsj 

R^O.  If  5N'  is  not  empty,  then  <p=  0  on  N.  The  decomposition 
(A1.2)  then  always  holds  with  0  differentiable  since  if  K(x)#0 

for  one  point  xON,  K(x)^0  for  all  xON  .  B 

This  lemma  is  sufficient  for  our  purposes  since  we  will  be 
dealing  with  topological  2— spheres  but  let  us  make  a  few  general 

2 

remarks  about  the  operator  D.  If  N'A'S  ,  d:=dim  ker  D  =  0  while  if 
N^'T2=S1XS1,  d  =  2.  (On  the  flat  torus,  ds8—  dx2-i-dy2,  KGker  D  iff 
KAS=c°nst  .) 


■ 
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The  method  of  the  preceding  proof  cannot  be  applied  to 
orientable  surfaces  which  are  conformally  equivalent  to  surfaces 
of  constant  negative  curvature.  However,  it  is  easily  verified 
that  D  is  elliptic  which  suggests  that  lemma  A1.3  remains  valid 
for  all  2— surfaces.  The  kernel  of  D  could  be  investigated  as 

follows.  Using  the  scalar  products  (K,L):=y jKABLABvolg, 

K,LUC  (Sg)  and  (0,^):  =  jf^'Vuvolg,  0,fGC°?(A1)  the  adjoint 
operator  is  found  to  be 

D*:C°0(A1)-^COO(Sg):fA^-2V(A^B)— VC^cgAB. 

*  k  A 

Thus  ker  D  consists  of  the  1-forms  Y7  f°r  which  Y7  dA  is  the 
generator  of  the  group  of  conformal  transformations  (the  group 
of  infinitesimal  coordinate  changes  which  leaves  the  conformal 

metric  gAB=g  g AB,  which  is  independent  of  arbitrary  scale 

changes  in  g,  invariant).  By  the  Atiyah-Singer  index  theorem 

* 

(c.f.  Palais  (1965))  the  analytic  index  i=dim  ker  D  —  dim  ker  D  is  a 
topological  invariant  which  can  be  calculated  from  the 
topological  properties.  One  finds  that  i=— 6  for  the  sphere  (the 
group  of  conformal  transformations  is  SL(2,CC))  and  i=0  for  the 
torus  but  unfortunately,  for  other  compact  surfaces  neither 
ker  D  nor  the  conformal  group  nor  the  topological  index  are 
easy  to  determine. 


i 


APPENDIX  2:  SOLUTION  OF  THE  LINEARIZED,  (2+l)-DIMENSI0NAL 
STATIC  FIELD  EQUATIONS  ON  THE  SPHERICAL  BACKGROUND 


The  following  result  was  first  given  by  Kunzle  (1971)  with 
the  slight  exceptions  that  he  assumed  0^3p^p  and  that  PA=71A 
was  not  necessarily  zero  while  below  we  assume  that  O^p^p  and, 
for  a  slight  gain  in  simplicity  with  no  loss  in  generality  (c.f. 

§1.4),  we  use  only  coordinates  for  which  y1A=0.  From  §1.4  it  is 
evident  that  we  can  regard  a  static,  perfect  fluid  spacetime  with 
a  fixed  equation  of  state,  fixed  central  potential  Uc  and  fixed 
central  pressure  pc  as  characterized  by  the  set  $=f(7,fl,W,p)£ 

where  all  these  tensors  are  functions  of  UH(Uc,0)  and  where  we 
will  drop  the  bars  for  the  rest  of  this  appendix  since  we  will 

always  be  dealing  with  tensors  in  the  (2+1)— dimensional 
formalism. 

In  order  to  solve  the  linearized  field  equations  on  the 
spherical  background,  consider  a  1-parameter  family  $( A)  where 
S(0)  is  the  spherically  symmetric  solution  and  follow  the 
linearization  procedure  as  described  in  §3.1.  We  will  factor  out 
coordinate  transformations  in  the  same  manner  as  was  done  for 

the  three  dimensional  case.  Specifically,  we  write  ciy  AB=:CAB, 

<5W=w.  On  the  compact  Riemannian  manifold  SC=U  \c)  there  is  a 
Berger-Ebin  (1969)  decomposition 

^ab=^ab^^ab’  ^B0A=O  (A2.1) 
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where  0  is  unique  and  0  is  unique  up  to  a  Killing  vector  field 
of  7.  Since  the  Lie  derivatives  of  all  the  equations  for  & 
simply  give  the  linearized  equations  for  variations  of  the  form 

for  f£<5,  and  since  =  0  (we  are  on  the  spherical 

background)  we  compute  the  variations  of  (1.59—1.62,  1.32) 
assuming  that 


^7aB~ ^AB’  ^b^a  0-  (A2.2) 

Since  the  central  potential  and  pressure  are  fixed,  <5p  is 
given  by  (4.10)  with  k  =  G  When  Sc  is  a  Euclidean  sphere  (so 
aAfi=0  ),  equation  (1.61)  together  with  lemma  A1.3  implies  that 

Qab=yq7ab-  (Recall  also  that  R=2WD(log  r)  and  R=2r  2  .)  The 
linearized  Einstein  equations  (1.59-1.62)  are  now  found  to  be 


(MW  Vn)w+WDw  +  ^W2D0=O  (A2.3) 

QW  23Aw+9AD0=O  (A2.4) 

(A+-^R)0+-^WfiD0=-2W_1(R  +  2p)w  (A2.5) 

A(W_1w)-WQD(W“1w)-4pW_1w-YW2DD0-^-MD0=O.  (A2.6) 


Using  (A2.3)  and  its  D  derivative  in  (A2.6)  one  can 


calculate  that 


* 
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Aw+W  DDw+(3M  2w"D(log  r))Dw+(2DM  — 2M  W  2+8MD(log  r)— 4p)w=0. 

(A2.7) 

The  three-dimensional  Laplacian  A  has  in  the  spherically 
symmetric  case  the  form 

Af=W2DDf  +  MDf+ Af  (A2.8) 

for  any  function  f  on  E.  After  some  calculations  (A2.7)  can  be 
seen  to  be  equivalent  to 


Aw=Fw  (A2.9) 

where  w=m(U)r  V,  m(U)=2  j^Mr^U  (  m(U)  agrees  with  the 
gravitational  mass  m  in  vacuo)  and 

F(U)=-DM+2MSW~2-6MD(log  r)+3p+W8.  (A2.10) 

Now,  (47AB)(U)eM^^  since  the  discontinuities  in 
derivatives  of  y AB  occur  only  in  the  normal  derivatives  D yAB. 

Thus  ^(U)^Mg_1 5+1  and  (A2.3)  then  implies  that  w(U)CMg_1 6+1  so 

w£Mg_12CC.  If  F(U)  were  everywhere  positive  it  would  follow 
immediately  from  (A2.9)  that  w  vanishes.  However  this  is  not 
quite  true  due  to  the  only  negative  term  — 6MD(log  r),  which 
makes  F  negative  near  the  center. 

To  make  a  closer  analysis,  let 
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w=EK£(U)Y,(^)  (A2.ll) 

where  is  any  normalized  linear  combination  of  spherical 
harmonics  of  order  i  so  they  satisfy  A*Yg+£ (£+l)Yg=0  and 

where  A  is  the  Laplacian  restricted  to  the  unit 
sphere.  Then  A=r  2A  and  (A2.9)  reduces  to 

wVDDK£+Mr2DKy=[r2F  +  f(f  +  l)]K£,  £=0,1,2,-.  (A2.12) 

We  first  investigate  this  equation  near  the  center.  Since 
w  is  C  ,  w  vanishes  at  the  center  and  in  fact  must  vanish  at 

least  like  0(z2)  where  z=(l/6)Mc(U— Uc)1//2.  Since  all  our  functions 
are  regular  at  the  origin  we  can  expand  in  terms  of  at  least  a 
few  powers  of  z.  Comparing  the  first  coefficients  of  the  powers 

t  t  +  l 

of  z  in  (A2.12)  shows  that  K^(z)=aogz  +0(z  )  near  the  center  so 

that  K0=K1=0. 

The  remaining  multipoles  can  be  shown  to  vanish  by 
showing  that 

Fp=F+£(£+l)r~2>0  for  all  2  (A2.13) 

everywhere  on  S.  For  then  (A2.12)  implies  AK^=2K^AK^+2W  (DK^) 
=2F^K2+2W2(DK^)2^0  so  that  K2  can  not  have  a  maximum  on  S. 
Since  lim  K£=0  this  implies  that  K£=0  on  all  of  S  so  w=0. 

That  0AB  then  vanishes  can  be  seen  as  follows.  Equations  (A2.3) 


1  I 
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and  (A2.5)  with  w  — 0  imply  that  A0+r  <f>  =  0  or,  restricted  to  any 

sphere  Sc,  A  0  +  0=0.  Since  the  first  nonzero  eigenvalue  of  A  on 
a  2— sphere  is  —2  there  are  no  regular  solutions  of  this  equation 

so  it  follows  that  0=0.  Since  V  0AB=O,  lemma  A1.3  implies  that 
0ab=O- 

It  thus  remains  only  to  show  (A2.13).  Let 

_  j  g  ,  g 

u=mr  ,  v  =  pr  ,  x=^Mr  u  .  (A2.14) 

Then  O^p^p  implies  0^|v^xu  and  (A2.13)  is  implied  by 

Fr2^18x2+15ux+u2— 18xVl+u2+v4— 6.  (A2.15) 

This  in  turn  is  implied  by 

(324x4-108x2+36)+x2(486ux-324v) 

+(30u3x  — 63x2u2+54uxj)+u4+180ux+12u2^0.  (A2.16) 

2  2  2 

The  first  term  is  just.  (I8x  —6)  +108x  ^0  while  the  second  is 
486u  (ux-v2/3)^0  and  the  third  is  greater  than 
ux(V54x- V/30u)2+18u2x2^0,  which  proves  the  relation  (A2.13). 


APPENDIX  3:  LINEARIZATION  OF  THE  (2+1)— DIMENSION AL 
STATIONARY  FIELD  EQUATIONS  ON  THE  SPHERICAL  BACKGROUND 

The  equations  which  we  present  here  were  first  given  by 
Kunzle  and  Savage(l980a).  They  are  mainly  of  interest  because 
of  their  similarity  to  the  linearized  3  — dimensional  equations  (on 

the  spherical  background)  involving  dh1  and  the  linearized 
static  equations  in  appendix  2. 

From  §1.4  it  is  seen  that  stationary,  rigid,  perfect  fluid 
spacetimes  with  a  fixed  equation  of  state,  fixed  central  potential 
Uc  and  fixed  central  pressure  pc  can  be  characterized  by  the 

set  5H(7,aW,hC,h,p,0,T)f  where  all  these  tensors  are  functions 
of  UC(Uc,0)  and  where  we  drop  the  bars  from  now  on  when 

there  is  no  danger  of  confusion.  We  write  6 W=w,  dh°=k, 

A  A 

dhA=*;A,  66  ,  dT=r  and  as  in  appendix  2  we  use  the 

Berger— Ebin  decomposition  to  write 

(57ab=^>ab+^+-7aB’  ^b0a=O  (A3.1) 

where  cfr  is  unique  and  is  unique  up  to  a  Killing  vector  field 
of  7.  Again  we  can  eliminate  the  coordinate  freedom 
corresponding  to  Lie  derivatives  of  the  field  equations  with 
respect  to  the  vector  field  f  and  compute  the  variations  of 
(1.48-1.58)  assuming  that 
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^  y  AB~  0  AB’  (A3. 2) 

A  Riemannian  2— sphere  S  has  no  closed  two-dimensional 
subgroup  of  its  isometry  group  I,  and  is  isometric  to  the 
Euclidean  2-sphere  if  dim  1  =  3  (c.f.  Kobayashi  1972  p.47). 
Therefore,  in  the  case  dim  I<3,  we  have  that  inside  matter 
A  A 

—c8  for  some  cEtR.  But  all  the  quantities  are  invariant 

A 

under  8  so  all  variations  will  then  vanish  under  a  coordinate 

A 

transformation.  Outside  matter  8  and  r  are  no  longer 

defined,  but  we  have  assumed  axisymmetry,  so  all  variations  will 
vanish  under  a  coordinate  transformation.  If  S  is  isometric  to 

the  Euclidean  2-sphere  (for  all  U  ),  all  quantities  become 

A  A 

functions  of  U  only,  6  =h  =0  and  all  variations  will  vanish 
under  a  coordinate  transformation. 

A  A  0 

On  a  spherical  background  we  have  6  =h  =h  =0  and  as  in 
appendix  2,  0ab=T°7ab'  Linearization  of  (1.57)  gives  Dt+t=0 

=3at  whence  r=rceU°  U  which  vanishes  since  we  are  restricting 

to  models  with  the  same  surface  temperature  Tb.  Note  also  that 

_ ^ 

the  integration  constant  c  appearing  in  dp=(p+p)(T  r-c)  must 

then  be  zero  (recall  pEMg_3n  but  <5pEMg_2n)  since  we  keep  the 
central  pressure  fixed  under  the  variation.  Using  these  facts  we 
now  find  for  the  linearized  Einstein  equations  (1.48-1.53) 


(MW  1+D)w  +  WDw+^  W2D0=O, 


(A3. 3) 


. 
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QW  2Oaw+3aD0=O,  (A3. 4) 

(A+^R)0+^WRD0=-2W_1(R+2£)w,  (A3.5) 

A(W“1w)-WOD(W“1w)-4pW“1w-^W2DD0-^MD0=O,  (A3. 6) 

D/cA+4/cA-W~23Ak-2W“1(p+p)e_6UTceUc£ABT9B=0,  (A3. 7) 

^[a^b]^,  (A3. 8) 

where  0  =  7  0AB.  Linearization  of  the  remaining  equilibrium 
conditions  (1.55)  yields 

V(A^b)=0  and  Di?A=0,  (A3. 9) 

when  p^O,  and  similarly,  linearization  of  Vrh  =0  gives 

Dk+(20W-1-MW~2)k+VA/cA=0.  (A3.10) 


Equations  (A3. 3  — A3. 6)  are  just  those  obtained  in  appendix  2  for 
the  static  case  while  (A3. 7  — A3. 10)  are  solved  in  chapter  5 


